OER CPA 
May 4761 
RATIONAL MECHANICS 


and 


ANALYSIS 


Edited by 
Cy TRUESDELL 


Volume 7, Number 2 


SPRINGER-VERLAG 
BERLIN-GOTTINGEN- HEIDELBERG 
(Postwerlagsort Berlin - 28. 2.1961) 


< 


Mechanicam vero duplicem Veteres constituerunt: Rationalem quae per De- 
monstrationes accurate procedit, & Practicam. Ad practicam spectant Artes omnes 
Manuales, a quibus utique Mechanica nomen mutuata est. Cum autem Artifices 
parum accurate operari soleant, fit ut Mechanica omms a Geometria ita distinguatur, 
ut quicquid accuratum sit ad Geometriam referatur, quicquid minus accuratum ad 
Mechanicam. Attamen errores non sunt Artis sed Artificum. Qui minus accurate 
operatur, imperfectior est Mechanicus, & si quis accuratissime operari posset, hic 


foret Mechanicus omnium perfectissimus. NEWTON 


La généralité que j’embrasse, au lieu d’éblouir nos lumieres, nous découvrira 
plutét les véritables loix de la Nature dans tout leur éclat, & on y trouvera des raisons 


encore plus fortes, d’en admirer la beauté & la simplicité. EULER 


.. ut proinde his paucis consideratis tota haec materia redacta sit ad puram 


Geometriam, quod in physicis & mechanicis unice desideratum. LEIBNIZ 


The ARCHIVE FOR RATIONAL MECHANICS AND ANALYSIS nourishes 
the discipline of mechanics as a deductive, mathematical science in the classical 
tradition and promotes pure analysis, particularly in contexts of application. Its 
purpose is to give rapid and full publication to researches of exceptional moment, 
- depth, and permanence. 


Each memoir must meet a standard of rigor set by the best work in its field. 
Contributions must consist largely in original research; on occasion, an expository 
paper may be invited. 


English, French, German, Italian, and Latin are the languages of the Archive. 
Authors are urged to write clearly and well, avoiding an excessively condensed or 
crabbed style. 


Manuscripts intended for the Archive should be submitted to an appropriate 
member of the Editorial Board. 


The ARCHIVE FOR RATIONAL MECHANICS AND ANALYSIS appears in 
numbers struck off as the material reaches the press; five numbers constitute a 
volume. Subscriptions may be entered through any agent. The price is DM 96.— 
per volume. 


Notice is hereby given that for all articles published exclusive rights in all 
languages and countries rest with Springer-Verlag. Without express permission of 
Springer-Verlag, no reproduction of any kind is allowed. 


For each paper 75 offprints are provided free of charge. 


‘os’ Aplin wpbhist: ve. 


Nl ped tet 


The Geometry of the Orbits of Artificial Satellites 


RAIMOND A. STRUBLE 


Communicated by E. LEIMANIS 


Abstract 


A novel approach to the study of the orbits of artificial satellites is presented. 
Emphasis is placed upon the basic geometry and other aspects of satellite motion 
which are of first importance to satellite engineering. The motion of the orbital 
plane as a rigid body is introduced and a non-elliptical orbit motion in this plane 
is defined. The plane orbit so defined possesses the very desirable feature of represent- 
ing a succession of satellite positions and hence reveals the true motion of the satellite. 
An analytical treatment yields a completely general second order theory of earth satellite 
motion which is suitable for engineering purposes. In the latter development, particu- 
lar attention is paid to the apsidal motion of the orbit and the concomitant resonance 
effects at the critical orbit inclination. The basic nonlinear features of the apsidal 
motion, which have not been recognized in earlier theories, are incorporated in the 
analytical development so as to produce a theory valid at all angles of inclination 
of the orbit. 


I. Introduction 

The orbit of an artificial earth satellite is a twisted space curve, wound about 
the earth in a complicated wave. The waves in this space curve are removed by 
introducing a precessing plane, called the orbital plane, upon which the orbit may 
be mapped as a plane curve. The latter plane curve is often referred to as the 
orbit itself though, in general, it is not a closed circuit even for bound orbits. 
Thus an additional artifice is required before one may properly speak of the 
motion of the orbit. Tradition calls for the introduction of an osculating ellipse 
which varies in size and shape throughout the satellite motion. Unfortunately 
(for engineering purposes), the osculating ellipse does not represent a succession 
of satellite positions and so reveals neither where the satellite has been nor 
where it is going. The motion (if any) of the osculating ellipse is interpreted 
as the motion of the satellite orbit though, of course, the two are somewhat 
unrelated. 

Here, a novel approach to the study of the geometry of satellite orbits is 
presented. The motion of the orbital plane as a rigid body is introduced, and a 
non-elliptical satellite motion in this plane is defined. The plane orbit so defined 
does represent a succession of satellite positions and hence reveals the true motion 
of the satellite. 

A second order theory of the motion of earth satellites is obtained. This 
theory is based on the gravitational potential of an oblate spheroid suggested 
by JEFFERIES [/] which contains, in addition to the fundamental, the second 
and fourth spherical harmonics. The second harmonic is considered to be small 
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of first order (in the oblateness parameter J), and the fourth harmonic is con- 
sidered to be small of second order. Expansions of the solutions to second order 
are obtained without the usual limitations as to eccentricity or orbit inclination. 
However, such generality calls for a number of innovations in the standard per- 
turbational techniques and concepts of celestial mechanics*. It behooves the 
reader, therefore, to take careful note of the details of the development and 
particularly of the definitions of the quantities introduced here. For ease of 
expression, considerable effort has been made to retain or parallel the traditional 
orbital concepts of astronomy, but this in turn may lead to confusion unless 
the definitions, as given here, are carefully adhered to. 


II. Orbital Geometry 


As indicated in Fig. 1, we refer the satellite motion to an inertial system, 
xyz. The z-axis is called the polar axis and the x y-plane is called the equatorial 
plane. Though our primary concern 
is with conservative motion of arti- 
ficial earth satellites, the geometrical 
Velocity vector treatment is quite general. 

At each instant, the orientation 
of the orbital plane is determined 
by the position vector of the satellite 
and the velocity vector of the satel- 
lite. The “‘line of nodes”’ is defined 
as the line of intersection of the 
y-os Orbital plane and the equatorial 
plane. Thus the argument of the 
line of nodes, 2, varies continually 
and corresponds to a “‘node’’ only 

as the satellite passes through the 

ee Meridian equatorial plane. Similarly, the 
Fig. 1. Orbital geometry angle of incidence 7 is defined as the 
instantaneous angle formed by the 

orbital plane and the equatorial plane. The angles Q and 7, of course, characterize 
the orientation of the orbital plane. The central angle f in the orbital plane 
locates the satellite relative to the instantaneous line of nodes. Thus the equations 


3-0x1s (polar axis) 


Line of nodes 


cot? =tanisin(p +2), (1) 
cos 3? = sini sin, (2) 
cos B = sin? cos (mp + Q), (3) 


which express the orbital coordinates in terms of the usual spherical coordinates, 
are treated as identities. 

We choose to think of the orbital plane itself as a rigid body with the origin 
as a fixed point. Then the motion of the orbital plane (as a rigid body) is an 


* We deal here with a troublesome degenerate system, discussed recently by 
BRILLOUIN [2], in which the variables do not separate. 
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instantaneous rotation about some axis through the origin. It is not difficult 
to see that this axis of rotation necessarily lies along the position vector 7 and 
merely reflects the rotation (if any) of the velocity vector about 7. This rotation 
is in turn responsible for the changes in 2 and 7. In fact, if the rotation y is 
positive as shown in Fig. 1 (counter clockwise as viewed from the origin), then 
from the geometry alone we derive the relations 


dQ sinBp dy 


Tio sie (Ue (4) 
and 

iis 2 

= cos p AY - ; (5) 


III. Equations of Motion 


The equations of motion of a satellite in spherical coordinates are well known. 
For an arbitrary potential U they may be expressed: 


SR cos al! 
rie sin? }- af) = ap” (6) 
dil pdt \E =, 25 ew _ 
al? ah a sind cos 0 (P| = Fae (7) 
OI ea O Ne ON et etn | Gann, CU 
ee Ue avg eee ee oe (8) 
Here we limit our considerations to a potential of the form*, 
eee eee iL: 9) a8 ( re 2) 
U pea hs (5 = cos? #) + wa (00S a : cos? } + —_}. (9) 
Because U is independent of gy, (6) pa one integral of the system 
sine go“? —p, (10) 


a 


where # is a constant (the angular momentum about the polar axis). Using 
this integral and equations (1) —(5), one may show ** that (7) reduces very simply 
to the first-order equation 

dy COSC mG) 

ie ita ek 11 

dt psin® ad — ae 


For the potential (9), (11) becomes 


e a Eee eens oo : (3sint — 7sin?isin?6)|}sinB. (12) 
The major variation in (12) is due to the factor sinf which is positive in the 
northern hemisphere and negative in the southern hemisphere. This in turn 
causes 7 to oscillate and 2 to increase steadily (at a variable rate, however). 


In fact, (4), (5) and (12) yield the equations 


dQ J g R* cos? 2 2D¢g RS cos? t rae eine 4 — cos 28 7 

dt ig or 7Tpyr® (3 — 7sin® 7 sin® B) | ( ), (13) 
aa til Tg R! cos?72 sin 2 2D g R®cos?i Sh) 7sin37 sin2B | sin 28 14 
dt ; pr 7pY* B : ) ( 


*x Here g, R, J, and D are constants. See [1]. 
xk For details, see [3]. 
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Q di 


~ and me a at each node (i.e., as the satellite 


Thus we see that both @ 
passes through the me ees plane) while @ a ’ also vanishes at maximum latitude. 


Clearly, 7 dbs is positive in the southern hemisphere and negative in the northern 
salle tend) but = is always non-negative. These facts could have been 
t 
inferred directly from the physics and the geometry. 
The radial equation (8) may be expressed in the form 
ay p*sects 00 
dt? 7 or 


or, with U given by (9), in the form 


Bie PESCON a. BREN, 37s OR cent ott? hi 
dt? Vad ie r* (15) 


6 . . . . = 
a on [sint i sintB — ° sin?7 sin? B + ei . 


Finally using (1)—(5) and (10), one may obtain the following equation for the 
central angle fp: 
oF 1+ a costisin?B +4 oe gh co s47(3 — 7sin7 sin?) sin? |. (16) 
Equations (14), (15), and (16), together with appropriate initial values, 
determine the exact motion* of the satellite in the orbital plane. In particular, 
it should be observed that the argument of the line of nodes [7.e., equation (13) | 
is not involved. However, to locate the orbital plane itself, equation (13) must 
be integrated (by a direct quadrature) using the solution of the system (14), 
(15), and (16). It is possible to reduce the fundamental system to third-order 
using 6 as the independent variable. The resolution of the third-order system 
is then followed by two direct quadratures, one to obtain ¢ and one to obtain Q. 
The derivation of the satellite equations (14), (15), and (16) from (9) is, of 
course, perfectly ngorous. To integrate the system, one could revert, at this 
point, to the traditional methods of celestial mechanics by introducing an osculat- 
ing ellipse. One would seek a solution for 7 in the form 
a(1—e?) 
1+e cos (B—o) 


ioe (17) 


wherein a, e and w, as well as f, are variables. Together with 7 and Q, they 
are referred to as the orbit elements: a is the semi-major axis, e the eccentricity 
and w the argument of perigee. If (17) is used directly in (14), (45), and (16), 
one has an extra degree of freedom: five variables, a, e, a, B andi and a fourth- 
order system. However, the requirement of ‘“‘osculation’? removes the extra 
degree of freedom. It amounts to nothing more than imposing the condition 
dy _ ae(1—e?) sin (B—o) 
dp [1-++e cos(B—o) |? 
upon the four variables, a, e, m and B. In the Keplerian case**, a, e, and w are 
constants, and (17) is merely the polar equation (in 7 and f) of an ellipse with 


(18) 


* Corresponding to the potential (9). 
ak Y= D=0. 
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the origin at one focus. Since in that case, 7 and 2 are also constants, the orbit 
is a fixed ellipse lying in a fixed plane. Clearly (18) is an identity in the Keplerian 
case, and imposing it in the general case provides for the osculatory interpretation 
of (17), z.e., the instantaneous elliptic motion which would prevail if all per- 
turbations were suddenly removed. The problem then reduces to the study 
(rigorous or otherwise) of the variations of the orbit elements. As indicated in 
the Introduction, the ellipse (17) does not represent a succession of satellite 
positions, and it is equally disturbing to find that the variations of the orbital 
elements do not represent the true motion of the orbit *. 

Here we are primarily interested in the orbit geometry and so introduce a 
new independent variable B, which is (geometrically) much like the central angle 
f but which preserves some of the mathematical simplicity of the system which 
would be lost if 8 were used as the independent variable. Equations (14), (45), 
and (16) become 


o = — 2 cost sini Jer u+ 2J*?D, g Rois = = sin?é sin®)] sin 2B, (19) 
G 
Nea Aneedae repens a 1S gare in27 sin2 
Bae Pea il oy cord 7 |“ &—_ y?({ — 3 sin?2 sin? 8) + 
k ad (; d : A 
4 ; a et . 2) 
fy ant (35 sint7 sin? 6 — 30sin?7 sin? + 3), 
7 ae seact f me 2i= u cos!z sin? B + 
eoncat (21) 
+ ALE costs u® (3 sin? 8 — 7sin?7 sin’ f)|, 
where 1 
Ua, (22) 
Vv 
dt k 
—— 2 
é dp pw ae 
an 
[2D = De (24) 


Actually equation (23) is the definition of B, and k is chosen so that # and B 
exhibit the same secular behavior. The latter is accomplished by choosing k 
so that the right-hand member in (21) has a mean value (with respect to f) 
equal to one**. 

We seek a solution of the system (19), (20), and (21) with w in the form 


u = — [1+ ecos(B —w) —Jo+ Jw], (25) 


wherein 7%), ¢, w, v and w are variables. This form should be compared with 
that in (17). Again, one has more degrees of freedom in the solution than the 
system (19), (20), and (21) warrants, but the extra degrees of freedom are “used 
up” by an appropriate separation of the terms that occur in (20). It will be 
possible to obtain a solution of the form (25) which is correct to second order 
in J and such that 7), e, and w exhibit only long-period *** perturbations. The 


* This is particularly evident for small eccentricities. 
*x* We do not know k a priori. It is determined as the solution develops. 
*xk Approximate period of the rotation of perigee. 
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latter, of course, also exhibits in general a secular perturbation. Thus all short- 
period* perturbations will appear in v and w. 

When e is of nominal size (e>>/J), the first two terms in (25) dominate, and 
the orbit is essentially an ellipse rotating slowly in the orbital plane (the apsidal 
rotation). In our solution, the mean radius 7), the mean eccentricity e, and 
mean argument of perigee « experience only long-period oscillations. The tra- 
ditional osculating values of these quantities in (17) would in addition reflect 
the short-period oscillations ** which here have been relegated to the additive 
terms Jv and /?w in (25). This form, in which the long and short-period per- 
turbations are separately displayed, serves to clarify the total description of the 
satellite motion. 

When e is small (of order J) the total orbit defined by (25) bears little re- 
semblance to an ellipse. In fact, Jv is then very nearly given by the single term 


4 7(2)’sin®i cos (28), (26) 


which becomes dominant as e->0. Since the latter is of double frequency, more 
than one perigee (z.e., local minima of the radial distance ***) may occur in a 
single revolution of the satellite. In addition, (26) does not share in the “‘apsidal 
motion’’ associated with the second term in (25), and so the orbit shape is 
continually subjected to major variations. It is clear that traditional orbital 
concepts based on the elements of an osculating ellipse become illusory for these 
nearly circular orbits. None the less, throughout this paper we shall refer to 
the variable e in (25) as the eccentricity of the orbit and to the variable m in 
(25) as the argument of perigee. In all cases, apsidal motion will refer to the 
variations in @, and we shall think of w as determining the “‘line of apsides’”’. 

It turns out that the apsidal motion is intricately coupled with the long- 
period oscillations of the angle of incidence of the orbital plane. This fundamental 
nonlinearity of the system has not been treated in earlier theories ****. Once the 
essential coupling is recognized, however, it is a rather straightforward matter 
to obtain the true characteristics of the apsidal motion. The latter are embodied 
in a modified pendulum equation. The familiar, nearly steady monotonic rotation 
of the line of apsides turns out to be the analog of a pendulum swinging com- 
pletely about its pinion, 7.e., a pendulum going “over the top’’, while the hitherto 
undisclosed apsidal motion for orbits near the critical incidence angle, 63.4°, 
is the analog of an “‘oscillating’’ pendulum. 

The quadratures called for in (23) to compute the time ¢ and in (13) to compute 
the argument of the line of nodes Q are not considered in this paper. Rather it 
is suggested that these quantities be introduced into the solution via the method 
of averaging developed in [7] and [8]. The latter is designed specifically for 
accurate predictions of the secular characteristics of the motion and for this 


purpose is probably superior to any of the perturbational theories which have 
been proffered. 


* Approximate period of the rotation of the satellite in its orbit. 
SS SOO (PA. 
ae The radial distance y may exhibit as many as thyee minima per revolution. 
‘* See, for example, the second order theories in [4], [6] and [6]. First order 
theories reflect none of the really interesting nonlinear features of the motion, 
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IV. The Second Order Theory of Earth Satellite Motion 
For /=0, the system (19), (20), (21) possesses the general solution 


C= 

k==COS ty 

B=B (27) 
og 6 naar 

eer COS? Zp 


Wa [1+ ecos(B — w)], 
0 
where 7), e, and w are arbitrary constants*. If only terms through the first 
order in J are retained, the system admits the solution 


1 
4 


: 2 
k = cos iol! — es COS? 25 
Yo 
uae 


1 =, 41 


7(F) sin 2%, |e cos (B+ w) + cos 2p 5 cos (38 —0)| 


B=B+2I(7 J 4ecostiy sin (B w) + 2e(4 — 2cos® 7) sin (B + @) + 


+ (14 — 3.cos?4,) sin 2B + sell — 4cos?%,) sin (38 — w) 


Aye cos? zy 4 = J (= (1 + _ (2 — 3sin®%9) (28) 


v9 pau Di V5 Me 
W = pa [1 +e cos (B — w) — J] 
0 
eee (= sin*ég| (2+ =| cos 28 + ecos (3B — w) + ae cos (4B — 2m) + 
0 3) 


2 = 
, (=) e? (2 — 3sin?7,) cos(28 — 2m) 
Yo 


d R \2 : 
arma! (x) ET ise a 


where now 7, ¢, and the initial value of m are arbitrary constants. 


If the first order solution is used in the right-hand member of (19), a 
may be determined to second order in J. One obtains 
is eae ti(2 J'sin 2i,[e sin (B +o) + 2sin 28 + esin(3B —@)] — 
dB iNT (29) 
— 5 (Fy sin 2i4|— 14 +15sin?7) + 2D, (6 — 7sin?%9)| sin 2w + Jy, 
a 


where y=y (B) represents numerous short-period trigonometrical terms which are 
easily integrated. The explicit second order term in (29) is a long-period term 


: . dw . 
which, however, must be given special treatment. For when aE is zero **, 
: 


* One constant of integration is determined by the condition f =. 


**x See the last equation in (28). The derivative a is always small of order /, 


but for i) =63.4° actually appears to vanish identically. 
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this term would apparently lead to secular variations in the incidence angle 7! 
Let us now introduce variations in the mean incidence angle 7) which satisfy the 
equation 


dt _ = 72(2)' sin 2i9|—14 + 15sintiy + +> Dy (6 — 7sin®) sin2w. (30) 
0 


apté‘ BS; 
Then all long-period variations in the incidence angle are absorbed (notationally 
for the moment) into 7)*. Since a is of order J?, the use of a constant 1, 


in first or second order terms is valid and, in particular, our zero and first order 
solutions (27) and (28) remain intact. Integration of (29) leads to 


t= + ; 7(2)'sin 2i9|e cos (B + @) + Cos OB 7 cos (38 — 0)| + 
amet ys (31) 
: Me ye 
Tees es Sin 219%, 
where y,=»,(8) is given in the Appendix and consists entirely of short-period 
terms. 


Using (31) and the first order solution (28) in (21), both k and 7 may 
be determined to second order in J. However, now k must exhibit long-period 
variations similar to 7, and, in fact, to second order in J, satisfies the equation 


dk d (COS 19) 


=27r(4) sin igsin 2ig|— 14 +15 sin*ig + 1° D, (6—7sin*ig)] sin 20. 

0 

Also k must absorb all the potentially troublesome, long-period terms which 
occur in (24). Specifically then we obtain 


k = cos ig — I * V'cos? to tm J? (=y cos? io ; (33) 
"0 Yo 
where 7, is given in the Appendix, and upon integration of (21), 
p=p+ £ (2) [4e cosy sin (6 — w) + 2e(1 — 2cos®%,) sin (B + w) + 
0 


ms _ 4 
+ (14 — 3.cos? 49) sin 26 + 3 e(1 — 4cos*%) sin (38 — w) an 


R\4 
+(ns 
Yo 
where »,=9,(B) is given in the Appendix and consists entirely of short-period 

terms. 


Having obtained 7 and k to second order, the integration of the radial equation 
(20) to second order may be considered. We assume a solution for ~ in the form, 


u=—-[1+¢cosB —o) — Ju+ Jrw), (35) 


where v is given in (28) and w will afford the second order corrections. In addition 
it will be necessary to re-define 7) and w and to permit long-period variations 


. us The arbitrary constant, formerly identified with ij, may now be identified 
with the initial value of 7). 
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in the mean eccentricity e. To second order in J, 7, will satisfy 


a7 p? d(sec? %9) 


Oya eRe ap .. 
=— aan, sin?iy| — 14+ 15sin?7 + 7 P16 — 7sin*t9) sin 2@. 


Guided by the extant first order solution, we shall seek a second order solution 


in which a is not greater than first order in J while = is not greater 


than third order. Similarly, we shall find that 4 is not greater than second 
de 


order while iP is not greater than third order. We note, in passing, then 
that each of “ a a and io is also of second order while each of 
diy d2h dey, 

—_.°, —— and —— is not greater than third order*. On the other hand, 
Ge CHEF d p 


d?w 


the additive term /?w in (35) and its derivative /?— iP is of second order and 


must be retained. 
f (31)—(36) are used in equation (20) and all terms through second order 
in : retained, then essentially six types of terms appear: (a) the second order 


2 
terms = and J*w; (b) explicit second order, short-period terms in harmonics 


appreciably different than the fundamental (7.e., terms of at least double fre- 
quency in B); ; (c) explicit second order, long-period terms; (d) non-harmonic 
terms; (c) explicit fundamental, harmonic terms (i.e. , cos (6 —w) and sin (6 —«) 
terms) ; (f) explicit nearly fundamental, harmonic terms (1.¢., cos (B+), cos (B—30) 
and cos(B+3) terms). 

The first three types (a), (b), and (c) are combined and determine w. The 
resulting w-equation is of the form 


—~ +w=%, (37) 


where ?3=%3 (B) contains only long-period terms (which are harmless in the 
second order equation (37)) and non-resonant harmonic terms. Hence the 
solution for w to zero order in J is written down immediately and is given in 
the Appendix. 

The non-harmonic terms (d) in (20) determine 7, according to 


1 _ gR? cos? ty ! J (FY (1 ! 5) (2—3sin®ig) + 5 (=) m. (38) 


2 4 
Yo p 27%) \ % 2 Vo Vo 


where 7 is given in the Appendix. We note that this is consistent with (36). 


The fundamental (e) and nearly fundamental harmonic terms (f) are classified 
either as cos(f —q) or sin (8 —w) terms. For example, cos (B+), expressed in 


the form cos 2m cos(B —w) —sin 2 sin(B —w), contributes both a cos (B —w) 
term and a sin($ —w) term. When all of the fundamental and nearly fundamental 


* See (30), (32) and (36). 
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harmonic terms are thus combined, these two classifications lead to the two 
equations 
de ae (= y e(ng sin 2m + 4 sin 4a) (39) 
dp 2 Yo 


and 


an = sI(F =) (5.cos*74—= 1) + re (n; +g cos 2m +n,cos4@), (40) 
(a 0 
where the quantities 73, j4, 7; and 7 are given in the Appendix. The first of 
these may be integrated by direct quadrature in terms of the integral of the 
second. Hence we are primarily concerned with the integration of (40). 
When i, is appreciably away from the critical value 63.4°, the first term on 
the right in (40) is truly of first order in J and hence dominates. The integration 
of (40), as well as other long-period variations*, to order /? is then straight- 
forward. One obtains the following long-period perturbations of the motion: 


o=o+|FI(F J (5.c0s%ig — 1) +15 J*(5 |B + 


z (7s Sin 20 + — sta sin 4a} 
BREN ag: (41) 

< jatareatst liga sin 21, 

Oa ao Tatler: 


14+ 15sin24, + as D6 7sin?i,)| cos 20 


5cos?7,—1 
R 2 
IG) ae 1 
Gre wv COST.) = 
0 Sco a (5 73 @ =— I M4 COS 400) ‘ 


where W , %o, and é) are integration constants. However, when 7) is near the 
critical incidence angle, the quantity (5cos?z)—1) is small and the first term 
on the right in (40) may not dominate. In such a case, the expansions in (41) 
can no longer be justified. In fact, a solution of (40) then evidently becomes 
sensitive to small variations in the mean inclination 7), and the integrations of 
(30) and (40) cannot be performed independently. This pair of equations is a 
two-variable version of the pendulum equation 


Pa 
a — J® A, sin 2+ J*/A,sin 4o = 0, (42) 


where 


i= ( Saat enti) 
A |— 14+ 45sin?2,+ > D, (6 — 7sin*i9)| — 216 (5 cos? ty — nf (43) 
A= 40, (= =| (5 Cos* 4 — 1), 


which is obtained by differentiating (40), using (30) and (36), and retaining only 
the dominant terms**. The long-period perturbations of the motion of an orbit 


near the critical incidence angle are governed by this modified pendulum equation 
and wall be discussed in the following section. 


SDE equations (30) and (39). 
** Due to the differentiation, the orders of all terms increase by at least one. 
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V. The Apsidal Motion of Orbits Near the Critical Incidence Angle 


As indicated above, the motion of the line of apsides of an orbit near the 
critical incidence angle is manifest in the modified pendulum equation (42). 
The solutions of (42) are qualitatively similar to those of the well known simple 
pendulum equation. (See [9].) Indeed, for 7, sufficiently near the critical value 
63.4°, Ay becomes negligible, and (42) reduces to the simple pendulum equation. 
Essentially all that is required here is an interpretation of the ‘‘pendulum motions” 
in terms of the present application. The phase-space solutions of (42) are the 
family of curves 


(3 
dp 
where A; is an arbitrary constant and are illustrated in Fig. 2 for 2,>0*. The 
appropriate solution curves are determined, say, by specifying the initial values 


of w and (the latter through equation (40) and the initial value of the 


ya A, CoS 2H — 5 F°dacos 40 = J* Aa, (44) 


mean inclination 7)), which in turn determine the constant A; in (44). 


The closed trajectories, interior to the heavy loops** in Fig. 2, represent 
the “oscillatory” or libratory motions of a pendulum. Here they also depict 
“oscillatory’” motions of the 
line of apsides about points 
of maximum latitude in either = - 
hemisphere. It is interesting to 
note that when A, is negative, 
these “‘oscillatory’’ motions 
center about the nodes in the 
equatorial plane***. On the 
other hand, the sign of A, in 
(43) is determined by the nume- 
rical value of D,, which for 
the earth is approximately 4 
(see [Z]), and turns out to 
be positive by virtue of the 


LW _gyig 
ap 


D,-term in (43). If the fourth a 2s 
harmonic in the earth’s poten- 

tial (9) were absent, somewhat =< 

smaller or negative, then A Fig, 2. Apsidal motion (arrows indicate increasing f) 


would actually be negative. 

As it is, the absolute value of 4, is but a half of what it might be if there were 
no cancellation of the effects of the second and fourth harmonics in the earth’s 
potential. As indicated in Fig. 2, the “oscillatory” motions of the line of apsides 
correspond to variations in w of amplitudes less than a/2 and to variations in 


a of amplitudes up to J? A}. The latter figure leads to variations in the 


* For earth satellites, J,>0 at the critical incidence angle and /,=0. 
*x*x The loops, of course, are separatricies [10]. 
xxx For 1,<0, the phase-portrait is similar to that in Fig. 2 except for a shift 
of origin. 
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mean eccentricity e and the mean inclination 1 of surprisingly large amplitudes. 
In fact from (30) and (42) * we have 
diy ine Ap RR ae aw 
4 a 


dp 48 J dp’ 
where es 
Were sin 2i,|— 14 + 15sin?%) + = D,(6 — 7sin* 1%) |, 
Al 
and hence 
ee. ee (ZY. 1 dw 
1G NO aaa eee 
Using the phase-space solution (44) *, this may be written 
lt Petich aan xy i (45) 
i / cos 2M)?, 
0 = ing F 8 (—) (ta — Ax 008 200) 


where there are two changes of sign per apsidal oscillation. The integration of 
(39) is only slightly more complicated. Using (42)* and (44)* we write 


oak Aas) a (rye ‘2 Fee a 
dp 2 
1 3 alk, 4¢e d*w (= y Beis 1 ido 5 
a HE J dp J a 7 4 ee ap) 
1 , (R\4 a As =y- e ad? sin 20 (a 
=|5 ae =| 22 ee 48 d Bp? (G) dp 


Again using (44) *, the last term may be ere in the rn 
+ na(3-) € J* sin 2m (Ag — A, cos 2m)? Ad 
70 dp 


2 


and the integration is now straightforward. There results 
7 nq A 1 
e=e+Je -) ia isa) (Ae A, COS 2m)? 
= a te = (A, — A, cos 2) 8. 


Thus both e and 7) experience ani perturbations of order /*! 


The non-closed trajectories in Fig. 2 represent the “rotary’’ or spinning 
motions of a pendulum. Here they also depict “‘rotary’’ motions of the line of 


apsides. The nearly level curves of Fig. 2, appreciably away from the w-axis, 


correspond to nearly constant values for - , 1.e., constant rates of rotation. 


These solutions merge with the solutions for orbits inclined away from the 
critical value of 7). Here one must retain the 4,-term in (44). 

The transition from “oscillatory” to “rotary’’ apsidal motion occurs along 
the heavy loops in Fig. 2. Since the motions represented by the trajectories on 
either side of the loops and the loops themselves are dissimilar, it follows that 
the solutions are not continuously dependent upon the initial conditions. It is 
clear that they are not continuously dependent upon the perturbation parameter / 
either. Though this behavior might appear to be somewhat disconcerting, it 
should not be considered at all unusual, for the loss of continuity (a fortiori, 
analyticity) under certain critical conditions is to be expected (see [2]) and is 

mOVVith A, —=\0) 
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closely allied to the inherent coupling which may exist between different variables 
in a system. 

To complete the integration of (44), one separates variables and obtains B 
as a function of w: 


Q 1 dw 
rh ny aa | 
‘ fe ate Vas —, cos 2H+4$/, cos 4m 
With A,=0, this integral is hyper-elliptic and thus not tabulated. However, 


for 7) near the critical value, A, is small, and it is sufficient to consider the simple 
pendulum case, A,=0. An elementary change of variable results in 


1/2—ay 


a) (47) 


3 (A +49)* +V1—2 sin? 0 
where eee 
j 
ae 8 
r AeA (48) 


For ki<1, either the plus sign or the minus sign in (47) prevails throughout, 
and (47) then corresponds to “‘rotary” apsidal motion. In particular, the period 


(in £) of one apsidal rotation * is 
ZI 


ee (49) 
m J* (Ay A)* 
satellite revolutions per apsidal revolution, where 
7/2 
de 


IC et 
("4) V1— 2 sin? 9 


0 
is the complete elliptic integral of the first kind of modulus hy. 
For ki>1, (47) may be transformed to 


Oe ep ey (50) 
2? J? 7 J 4V1—kesin? 6 
where, say, a 
siny =| ede J sin (= 0) 
, soe 2 
an 
yo At (51) 
Die 


The integral in (50) is again a Legendre elliptic integral of the first kind and 


corresponds to “‘oscillatory’”’ apsidal motion. In particular, the period (in 4) 
of one apsidal oscillation about a point of maximum latitude is 
2? K(k) 


fe : 
me eft al 


, (52) 


satellite revolutions per apsidal oscillation. 


This work was sponsored by the Office of Ordnance Research, U.S.Army. The 
writer wishes to acknowledge the valuable assistance of Messrs. W. F. CAMPBELL, 
J. E. Fretcuer and S. M. Yronouris in the preparation of this paper. 


* i.e., a change of 27 in w. 
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Sulla determinazione della densita neutronica critica 
nella teoria matematica della pila atomica a fissione 


ANTONIO PIGNEDOLI 


Memoria presentata da DARIO GRAFFI 


There is a large bibliography on the mathematical theory of the retardation 
and diffusion of neutrons within the moderating medium of an atomic fission pile. 

The methods used for such study must allow for the differences in reactivity 
that take place in the reactors, chiefly while starting but also because of the lowering 
and raising of the control bars, and must provide a quantitative model for the behavior 
of the fast neutrons produced in nuclear fission and slowed in the moderating medium 
until they fall within the range of “‘thermal’’ energies. The neutrons endowed with 
such energies are trapped, in turn, by the uranium nuclei so as to produce new nuclear 
fissions, and so on, by the well known “‘chain reaction’’. 

In this paper, in order to determine the “‘neutron density”’ relative to the ‘‘critical 
condition”’ for moderating media of various geometrical forms, we shall study the 
general equation of diffusion, particularly in the case of a multiplying medium as 
occurring in the atomic fission pile. 

Study of neutron diffusion in a moderator requires an appropriate solution of 
the Maxwell-Boltzmann integro-differential equation of the ‘‘transport’’ theory or 
an appropriate solution of the partial differential equation of diffusion, derived from 
the “‘phenomenological”’ theory. Both of these mathematical theories seek to determine 
analytically the way the fast neutrons produced by nuclear fission behave in the 
circumstances described above. 

In practice, for the study of neutron diffusion in relation to the operating of 
nuclear fission reactors, the phenomenological theory yields a completely satisfactory 
model. Within this theory, we consider determination of the neutron density at 
“critical conditions’”’ for the nuclear fission reactor, studying the problem from the 
general point of view and also, in particular, for spherical and cylindrical moderators 
with elliptical or circular boundary. The problem is studied also from the point of 
view of the “age theory’, introducing the solution of the differential equation of the 
“age theory’’ as the source function in the differential equation of thermal neutron 
diffusion. 

In the first part of this paper we give references to earlier works on the subject. 
The reader should refer also to the bibliography printed in the volumes of the review 
“Energia nucleare’’ edited by Crise at Milan and to the one published in the Pro- 
ceedings of the two international meetings held at Geneva on the use of nuclear energy. 


1. Vasta é la letteratura concernente la teoria matematica del rallentamento 
e della diffusione dei neutroni in un mezzo moderatore, teoria interessante, 
quindi, il funzionamento di una pila atomica basata sul fenomeno di fissione 
nucleare. Ne diamo in calce un ampio cenno bibliografico }. 


1 Cfr. Fermi, E.: Ricerca scientifica 1936. — AGENo, M.: Nuovo Cimento 1943. — 
AMALDI, E., D. BocciARELLi, B. FERRETTI e G.C. TRaBAccHI: Ricerca scientifica 
1941. — Wick, G. C.: Rendiconti della R. Accademia dei Lincei 1936. — Wick, G. C.: 
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I metodi usati per tale studio debbono tenere conto, naturalmente, delle 
variazioni di reattivita che si verificano nelle unita reagenti nucleari stesse 
soprattutto in fase di avviamento o corrispondentemente alle manovre, co- 
mandate dall’esterno e consistenti nell’affondamento o nell’estrazione delle sbarre 
dicontrollo. I metodi in questione debbono, in definitiva, concretare uno schema 
quantitativo soddisfacente delle vicissitudini subite dai neutroni veloci emessi 
per fissione nucleare, per giungere, attraverso rallentamento in un mezzo mode- 
ratore, alle energie cosiddette ,,termiche™. 

In corrispondenza di queste ultime, i neutroni vengono, come ben sl sa, a 
loro volta catturati da nuclei d’uranio per dare luogo a nuove fissioni, e cosi 
via, secondo il meccanismo della ,,reazione a catena”. 


In questo lavoro, con lo scopo di potere assegnare poi la ,,densita neutronica 
relativa alla condizione critica‘’ per pile dotate di mezzi moderatori di varie 
forme geometriche, prenderemo in esame l’equazione del fenomeno di diffusione 
in generale, indi con particolare riguardo al caso di un mezzo moltiplicante, 
di forma sferica o cilindrica (a sezione ellittica o circolare), per passare infine 
allo studio del problema in questione in una pit generale impostazione che 
proviene dalla ,,age theory”. 

I metodi usati per lo studio del fenomeno di diffusione dei neutroni in un 
mezzo moderatore si basano o sulla ricerca, sotto opportune condizioni, di soluzioni 
della equazione integro-differenziale del tipo di MAXWELL-BOLTZMANN della 
teoria del ,,trasporto“ oppure sulla ricerca, con opportune condizioni supplemen- 
tari, di soluzioni della equazione differenziale alle derivate parziali di diffusione, 
discendente dalla applicazione della cosiddetta teoria ,,fenomenologica™. 

Le teorie matematiche in questione puntano tutte a stabilire, col maggiore 
rigore possibile, quale sia l’aspetto analitico delle vicissitudini subite dai neutron 
veloci emessi per scissione nucleare per giungere, attraverso rallentamento e 
diffusione per urto in un mezzo moderatore, all’energia termica, alla quale i 
neutroni vengono in parte ulteriormente diffusi ed in parte catturati dai nuclei 
d’uranio per dare luogo a nuove scissioni. Agli effetti pratici, per lo studio dei 
fenomeni diffusivi dei neutroni in relazione al funzionamento dei reattori nucleari 
a fissione, l’impostazione derivante dalla teoria fenomenologica costituisce uno 


Zeitschrift ftir Physik 1943. — Ariry, N.: Dansk. Vid. Selks. Medd. 1938. — 
AMALDI, E., e G.C. Wick: Appunti di Fisica nucleare. Roma: Pioda 1946. — 
MarsHak, R. E.: Reviews of modern Physics 1947. — Fermi, E.: Science 1947. — 
VERDE, M.: Nuovo Cimento 1947. — PiGNEepott, A.: Atti del Seminario matematico 
e fisico della Universita di Modena 1947—48. — Satvertt, C.: Nuovo Cimento 1949. — 
GALLong, G., e C. Satvetti: Nuovo Cimento 1950. — PrcNEport, A.: Rendiconti del 
Seminario matematico di Torino 1950—51. — Rendiconti del Seminario matem. e 
fis. di Milano 1951. — Annali di matematica pura ed applicata 1951. — Rendiconti 
della Accademia delle Scienze dell’Istituto di Bologna 1954. — Journal of Rational 
Mechanics and Analysis 1955. — Rendiconti del Seminario matem. dell’ Universita 
di Padova 1956. — Grasstone, S., e M. C. Eptunp: The elements of nuclear reactor 
theory. New York: Van Nostrand 1957. — Murray, R. L.: Nuclear reactor physics. 
Prentice Hall, U.S.A., 1957. — Mazzoent, F.: Ingegneria nucleare. Milano: Hoepli 
1956. — Prrsico, E.: Gli atomi e la loro energia. Bologna: Zanichelli 1959. — 
I lettore terra presente anche la bibliografia esistente nei volumi della Rivista 
, Energia nucleare“‘ edita dal C1sr di Milano e quella pubblicata nei ,, Prodeedings‘ 
delle due conferenze internazionali di Ginevra sulla utilizzazione dell’energia nucleare. 
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schema di lavoro del tutto soddisfacente. Secondo tale impostazione, ci occupe- 
remo qui della determinazione della densita dei neutroni in condizioni ,,critiche‘‘ 
per il reattore nucleare a fissione, considerando il problema sotto il suo aspetto 
generale e per moderatori di forma sferica e cilindrica (a sezione ellittica 0, in 
particolare, circolare). 


2. L’equazione generale che regge il fenomeno di diffusione dei neutroni in 
un mezzo moderatore pel quale la probabilita di diffusione sia maggiore della 
probabilita di assorbimento si scrive, come é noto, operando il bilancio del numero 
di neutroni per unita di volume, precisamente notando che la variazione della 
densita o dei neutroni nella unita di tempo é data, per unita di volume, dal 
numero N, di neutroni prodotti, diminuito del numero N, di neutroni sfuggiti 
al contorno, diminuito ancora del numero WN, di neutroni assorbiti. Tale equazione 
é, dunque: 


CQ = om 
(1) 2¢ _N—N_N, 


In termini della funzione incognita ,,densita dei neutroni 0“ dipendente dal 
posto e dal tempo, l’equazione differenziale (1) diventa: 


(2) Alec +D-A,o0—o,v0, (4, =operatore laplaciano) 


dove S é la funzione, dipendente dal posto e dal tempo, ,,delle sorgenti neutro- 
niche“, a, € la ,,sezione d’urto di assorbimento™, v é il modulo della velocita dei 
neutroni termici diffusi (che si suppone costante), /, il ,,ibero cammino medio 
di diffusione“ e dove, infine, il coefficiente costante D vale D=vil,/3 ?. 

Nel caso particolare del regime stazionario di neutroni termici, per cui é 


SO: a = 0, la (2) si riduce alla: 


1 l, ee eee 
(3) A,o 2 es = 0, (L? =m ae ar (Cale ati) 


30, va, 6 


2 Essendo N il numero degli urti verificantisi nel mezzo, indicati con s; (¢= 1, 2, 
., N) gli spostamenti vettoriali di cui si compone il cammino r di un neutrone, sara: 


N N N 
PPP = Ly By BPS a Dis ae 2 OF: 
= g—1 =1 Ee) 


Prendendo la media di ambo i membri, si ha: 


CO 
Se Ss Z 
(77) medio=N -° (s?)medio=N - [= exp| ] Jas =2N ij. 
: 5 s! 


La quantita v /,/3 si indichera nel seguito con D e si chiama , coefficiente di diffusione“‘ 
Si definiscono anche una ,,lunghezza di rallentamento‘’ L, ed una ,,lunghezza di 
migrazione’‘ M date rispettivamente da: 


5 ae 
Here: (7?) vel. medio » M? = % (7 7”)tot. medio, 
essendo 


(77) tot. medio = (v7) term. medio “> (7) vel. medio + 
con 


5 g AG 2 23 = [2 
(7?) term. medio = 2 Neem. Ike , (7?) vel. medio = 2 Neel. Ik. 


2 
(7? \tot. medio = 2 (Merm. ae Nvel.) Bice 
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dove la quantita L, dimensionalmente una lunghezza, viene appunto chiamata 
,lunghezza di diffusione. Alla (2) va associata una condizione al contorno. 
In generale, tenuto conto del fatto che, attraverso il contorno del moderatore, 
si verifica una certa diffusione di neutroni, l’equazione generale di diffusione 
(2) dovra soddisfare sul contorno (£) del mezzo moderatore all’equazione: 


(4) As +ho=0, (n = normale esterna, / = cost. nota), 


che esprime, come si sa, la classica condizione al contorno di tipo diffusivo. 
(Non opereremo in questo lavoro considerando, come si fa spesso, il cosi detto 
,contorno estrapolato, che viene sostituito al contorno vero e sul quale con- 
torno estrapolato si impone una condizione di annullamento della funzione in- 
cognita @ *.) . 


3. Ponendo, anzitutto: 


o(x, y, 2,4) =exp(— a, v2) R(x, y, z, 2), 


0) S(x, y, 2, t) =— Dexp(— a, v2) Q(x, y, 2,2), 


la (2) si trasforma in una equazione del calore con ,,funzione delle sorgenti 
termiche Q(x, y, z, t), cloé nell’equazione: 


ANS ONE an Waco aes ; 
(6) AcR (x, y; 4, t) D et : Q(x, Ws 4, t) ; 


ed a questa va associata la condizione al contorno: 


oh +hR=0, su (£) per ogni ?. 


Ponendo quindi: 


Q(x, y, 2,t) =f, (x, v, 2) coswt+ f.(x, y, z) sinwt, 


(7) 
R(x, y, 2, t) = U(x, y, 2) coswt + U,(x, y, 2) sinwt, 


per le funzioni U, e U, si ottiene il sistema di equazioni differenziali alle derivate 
parziali: 


(8) A,U,— > =h, A,U,+-| U=he, 
il che ¢ come dire che, per la funzione U=U,+7U,, ponendof=f,+7f2, (i= ¥—1), 


vale l’equazione differenziale alle derivate parziali, relativa alla regione @ oc- 
cupata dal moderatore: 


(9) A, U(x, y,2) +AU(x,y,2)=f(%,%,2),  (A=t@/D). 


A questa va associata la condizione al contorno 


(10) fy LOO, suey: 
on 


° Cfr. per es. GLASsTONE, S., e M. Eptunp: Nuclear Reactor Theory, gia cit 
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1] problema in questione si riduce alle equazioni integrali. Precisamente, per la 
funzione U, vale con chiaro simbolismo, nel campo complesso, l’equazione inte- 
grale: 


(141) «-U(P’) =< [or P)U(P)d@— qn f P) f(P) dp, 
‘6 6 


in cul G(P’, P) é la funzione di GREEN, ed ogni soluzione della quale soddisfa 
all’equazione differenziale (9) ed alla condizione ai limiti (10). Nel campo reale 
la (11) si scinde nel sistema di equazioni integrali: 


’ ; 1 , 
U=— 225 [ CP. P) UP) aS — aL [G(R P)A(P) dG, 
€ C 


(12) ; 
U,= 2, f GP, P)U,(P) 46 — 4 { GP", P) f,(P) 4p. 
E 6 


Se A=7w/D non é un autovalore del nucleo della (11), essa soddisfa, nella regione 
occupata dal moderatore, con la condizione sopra detta sul contorno (£), al 
teorema di unicita, cloé ammette in @ una ed una sola soluzione continua soddis- 
facente alla data condizione ai limiti. 


Si consideri allora il problema omogeneo: 
Aa Ua Vet) He AO %,.9, 2) = OF — Ie: 


(13) aU 
ai WU 30, sure): 


Se A, sono gli autovalori reali e positivi del parametro A e se U, ne sono le cor- 
rispondenti autofunzioni costituenti un sistema ortogonale, tali dunque che é 


(14) fi, Ud C@=054 per k= 7; 
¢ 
si suppone / sviluppabile in serie di autofunzioni U, della forma: 
(15) f(x, ¥, 2) = a, U,(*, 9,2), % ee y, 2) U,(%, y, 2) 4%, 
k 
e ponendo poi 
(16) Olt 9nt) = 2A Us (44 YA), 


si ricavano i valori delle costanti A,, che sono forniti dalla: 
Ap 

A,= — 7. 
2) he A=); 
Occorre dunque determinare gli autovalori e le autofunzioni del problema omo- 
geneo considerato (13). Esistono, come si sa, infiniti autovalori del parametro 
A tutti positivi e disposti in successione crescente con unico punto di accumu- 
lazione -+ coe. 


4. Veniamo ora, pit in particolare, al problema della pila atomica di fissione. 
Lo studio di tale problema si riduce, nello schema del nostro lavoro, allo studio, 


4 Vedi Porncar#&, H.: Théorie analytique de la propagation de la chaleur, cap. XV. 
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con la solita condizione al contorno, della equazione differenziale di citioeeue 
per un mezzo moltiplicante, pel quale si fa ordinariamente, come é lecito, 1 ipotesi 
semplificatrice consistente nel supporre che i neutroni prodotti nel processo di 
fissione siano esclusivamente neutroni termici?. 

Se si suppone, dunque, che, per ogni neutrone termico catturato, siano pro- 
dotti K,, neutroni termici, la funzione delle sorgenti neutroniche S diventa: 


(18) =O, U0; 
e l’equazione differenziale di diffusione (2) si trasforma nella seguente: 
(19) £8 = D- Aye + 0¢0(Keo— 1)0=0, 


per la quale deve valere, al solito, la condizione diffusiva al contorno: 
oer hoa 0, st (2). 
Essendo o, v/D =1/L?, la (19) puo scriversi: 


Ke 1 ee eee 
. 0+ Ae Om Sa 


e per tener conto, sia pure in approssimazione, del fatto che i neutroni di fissione 
nascono come neutroni veloci, sostituiremo, come si fa, Z con la lunghezza di 
,migrazione‘‘ M, considerando dunque, per il fenomeno diffusivo in parola, 
Vequazione differenziale: 


24) Avot Kt =, 22, (Kx = Ket) 


dove ty=1/o,v prende il nome di ,,vita media di un neutrone in un sistema 
infinito‘‘. Ricordiamo ora che, per un reattore nucleare a fissione, si chiama 
,coefficiente di moltiplicazione effettivo’’ AK, il rapporto: 


(22) fore numero tot. dineutr. generati 
a * numero tot. di neutr. assorb. + numero tot. dineutr. perduti ° 


Un reattore si dice poi ,,critico‘’ od anche ,,di dimensioni critiche‘‘ quando é 
K,=1. Per K,<1 l’unita reagente tende a ,,spegnersi‘‘; per K,> 1, invece, 
il numero di neutroni presenti nella unita reagente e la potenza vanno crescendo, 
tanto che, se le energie di fissione prodotte sono molto alte, il reattore si avvia 
a diventare esplosivo. Lo studio delle condizioni critiche di un reattore, di 
per se della massima importanza, esige in pratica grande accuratezza. Basta 
invero per K, una differenza da wno di poche unita dell’ordine di 10~? perché 
il regime del reattore nucleare varii con grande rapidita. Ora é importante 
tenere presente che, in un reattore nucleare di dimensioni critiche, la densita 
neutronica diventa una grandezza di tipo conservativo, nel senso che sussiste, 
in condizioni critiche, equilibrio fra i neutroni che vengono generati e quelli 
che sfuggono al contorno o che vengono catturati. La funzione di distribuzione 
della densita neutronica g diventa dunque indipendente dal tempo. Inoltre, 
come si dimostra, non si ha pit a che fare con le armoniche pit alte della funzione 


» Cfr. Litter, D. J., e J. F. Rarrre: I reattori nucleari, gia cit. 
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di distribuzione della densita neutronica stessa. Tale distribuzione, indipendente 


dal tempo, si consegue determinando la soluzione prima armonica della equazione 
differenziale: 


(23) Ay Ory He = 0. AEG Sa K oe 


che si ottiene dalla (21) facendovi (stazionarieta della oye = = 0, con la condizione 
al contorno: 


(24) “8 + he= 4, su. (HE); 


5. Daremo qui la soluzione del problema della determinazione della densita 
critica nei casi di moderatore a forma sferica e cilindrica circolare e nell’ipotesi 
generale di dipendenza della 9 da tutte le variabili di posizione, ma, soprattutto, 
assegneremo la soluzione per il caso della pila cilindrica a sezione ellittica. 

Nel caso in cui il mezzo moderatore occupi una regione sferica di raggio a, 
esprimendo l’operatore laplaciano in coordinate polari 7, #, y ed operando suc- 
cessivamente la classica separazione di variabili del tipo: 


(25) o(7, ¥, 9) =A(r)- (9,9), 


a (23) si spezza nelle due equazioni differenziali: 


Lim 52 1 eY ~ 
ae = foe (sin oe a ante O92 LGW == ©. 
i EER): + K? jo= (a Crcostante) . 


R dr 


La prima di queste ammette la soluzione generale: 


2n-+1 
(27) y”) — > m, P™(cosy,), 

s=1 
dove le m, sono 2” +1 costantie si ha cos y,=cos ? cos #, + sin # sin #, cos (p—9,), 
essendo #, e gy, le coordinate polari di 2”-+-1 punti appartenenti alla superficie 
sferica di raggio unitario e P™) il simbolo dei polinomi sferici di LEGENDRE di 
prima specie. La condizione al contorno da imporre sulla superficie sferica, cioé 


si riduce alla: 
aR Qe = 
(29) iets RA =O Wiper 1 —a- 


La seconda delle (26) si puo dunque scrivere: 


(0) OF) + Ky R= n(n +1) 2 


Questa ammette la soluzione (regolare nell’intorno dell’origine) : 


(31) Ba (2G ep (14 — #2)" cos(¢r K) dt, 
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la quale é scritta a meno di una costante addittiva inessenziale. La condizione 
ai limiti (29) diventa, dunque: 
1 
a” K” i (1— #)"|" cos(taK) —#K sin(taK)| dt + 
0 1 
4 hark” [ (1 — #2) cos(ta kK) dt=0O, 


0 


ed, escludendo I’autovalore banale K =0 e ponendo tak =f, fornisce la seguente 
equazione trascendente determinatrice degli autovalori del parametro K: 


ak 
hs 2 \n 2 
(32) | (1 oz) [(n + ah) cos — Bsinf] dB =0. 
0 
Questa, per »=1, si riduce alla seguente semplice equazione trascendente: 
, , ah—2 
(33) tg(aK)=ak ah—2+a? k2° 


Con cio si é in grado di scrivere le autofunzioni della equazione differenziale 
diffusiva (23) con la voluta condizione al contorno e, in particolare, di assegnare 
la densita neutronica critica. 


6. Il caso, particolarmente interessante dal punto di vista fisico-matematico, 
del moderatore a forma di cilindro ellittico, costituisce problema che si puod 
risolvere, prescindendo dalle funzioni di MATHIEUv, con l’ausilio delle ,,funzioni 
trascendenti del cilindro epicicloidale‘‘ di AGOSTINELLI ®. 

Sia (e) la superficie laterale del cilindro; supponiamo poi che la base inferiore 
del cilindro stesso giaccia sul piano coordinato Oxy e che la base superiore si 
trovi alla quota z=/. Cercando di soddisfare alla (23) ponendo: 


(34) @(%, y, 2) = u(x, y)-Z(2), 
si ottiene il sistema differenziale: 
A,u(x, y) + pu =0 (4. Disgcaeie 


(35) amen 
sd 4 @Z(2)=0 (K*=p2+¢9); 


al quale vanno associate le condizioni ai limiti: 


(36) M+ hu=0 su (e), 
i dZ z / 
(30’) a ea per z= 0; £2. P20 per 210. 


Riferiamo i punti del campo piano ad un sistema di coordinate ellittiche (é, 7) 
ed indichiamo con c¢ la semidistanza focale delle ellissi ed iperboli omofocali 
E=cost, 7=cost (OSESE,, OSS 2a con & valore di & sul contorno (e)). La 


®§ Vedi AGOSTINELLI,C.: Atti della Accademia delle Scienze di Torino 86 
(1951—52) e suoi lavori precedenti ivi citati. 
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trasformata in coordinate ellittiche della prima delle (35) é la seguente equazione 
differenziale : 
eu a 


(37) aD 4 - ; p* c? (cosh 2€ — cos 2n) u=0. 


Quest’ultima, ponendo r=exp é, si trasforma nella equazione differenziale che 
si ottiene dalla 


1 01 1 ¢ —2(n n 
(38) a4 5+ a tp? n2a5 [1 +0207 2"*) 2407-49 cos(n +1) n]u=0, 


. . . C 
con m intero diverso da zero, facendovin=1, «=—1, ape € che ammette 
soluzioni dei quattro tipi seguenti: 


(39) V= eeu ie UE p) = Soos( (2s 7) ) Sean F 2m, AS p); 


(40) V= > coma Uom+1(§, 7, Pp) = = 21008 (2s 1) )9 X Coma F 2m-+1,s (So) 
(41) V*¥—) ch, Un (én, ) = 7a (287) 2; Cm Gom,s(E, B), 
m=1 Sz m=1 


(42) Vie a 2) Eimst Ud im +1( (S 1); p) = oy sin ( 2S st 1) \n >: emit Gom+1, es p) 


m= 


con le Com, Com+i> lm» Clm+1 costanti arbitrarie e dove si ha: 
foe) 
(43) U2m(&,% P) = 2 Fans, p) cos (2s%), (m= 0,1, 2,..-), 


(44) toms €s% #) = 5 Fair el6. 8) cos (2s + 1)%, (m = 0, 1, 2,...), 


(45) wba (E, 1, P) =D, Gon o(€. 6) sin (257), AA. 2e), 
(46) tna (6.14 ®) = 3, Gam: BRE prem (2s tin. = 0,120) 


essendosi posto (J,, funzione di BEssEL di prima specie di ordine m): 

(47) FamolE, 2) = In (32 6#) In(BD Ce); 

(48) Fam, s(E, ?) = Ings (EP 6C*) In—s(BP CC *) + Im—s(BP Ce) InaclBb Oe"), 
49) Gam cl6sP) = InleP oC) Imas(BP CS) — In—a(BP°*) Inte (BP CE), 
50) Foyt (6.0) Seen PCO) I BPO ) + Ince P CO) Imgeri(aP ce *), 
54) Goya (GP) = Inge GPC) Ip gE PCE) — Ina s (PCO) IntsriaP ce). 
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Sulla superficie laterale del cilindro ellittico considerato, cioé per =&, deve 
essere: 


av eve +e Me ADV ea a 
(52) S-+hV=0, +h4V=0, 5 +hV¥*=0, 4 +hV*=0. 


Le condizioni (52) si esplicitano come segue: 


(53) Seos| (2s) ieee s(Go,P PIA Oem, Be 2m ls a (£5, 2) —= OF 


m= 


(54) © cos( (25-4) ine sie al s(&, p +h cos ( (28%) 2 y Cam 2m, s(€o, 2) == 


m= m= 


(55) Ssin( (25%) ) Dn Gh (Eas A) +h Ysin( (2s) ale ie 2m, {eps p)= 0; 


m= m=0 


> sin (2s +. 1) 1] 3 Caras Goyiatlens p) =i 
(56) s—0 m=0 


AP h> sin (2s = 1) Ds C3m41 Gom-s,s (So p) = 0, 


s=0 m=0 


dove gli apici indicano derivazione rispetto alla variabile €. Poiché i coefficienti 
rispettivi di cos(2s7), cos(2s+1)y, sin(2sy), sin(2s+1)n, debbono essere nulli 
per ogni valore dell’indice s, risulteranno, per la determinazione degli autovalori 
del parametro #, i quattro seguenti sistemi di infinite equazioni lineari in infinite 
incognite: 

a) [Fo,0 (£0, Pye hFy 9 (Eo, PD) | + cg [Fo (EqsP yak hE, o(&o, p)|+ 
TT Com Bp o(&o, p) sh h Pom, o(&o, p)] a a 0, 
Co [Fo.1 (Eo, p) thy 1(§o, P)] + ¢, [Fo (9,2) +h ales: Pi 


(57) Peal Seley 2m, sa ae 


ARP) A +e CR B+ AR, oA + 
Walre See ae pS rls Sie 0) 


[Fo (6o, 2) +4F, 9 (Eo, £)] + ¢3[Fro (Eo, P) +hFs 9 (Eq, P)] + 
“Ht Com 43 [F3m41,0(E, D) + AF ym+1,0(o, ~)] +--+: =0, 
[Fis (Eo, 2) +E, 1 (Eo, B)] + cg [Fe 1 (Ey, 2) + hF3 1 (&>, £)] + 
(58) es piece) eg oa pe give 


CLF (Eo, 6) +h, ae + cs [F% , (Eo, 2) eee )| + 
a ae Com+1 [se omaieae (ean p) Sima S yevrecwe (Shr p)| teem 
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60 [Go,1 (Eo, 2) +4 Gor (Eo, )] +8 [Gs,1 (Ey, b) HAGs 1 (Ey, B)] + 
Bore ag Cam cop Pe! Ga (eon Poe 0. 
C0 [Go,2 (Eo, b) +h Go, (&,p)| +cF [G22 (Fo, d) +hG, o(&, p)] + 
(59) -++ + 63m [Go m,2 (Eo, P) + Gamo (Eo, P)] +---=0, 


C8 [Go,s(Eo, 2) +4 Go, « (Eo, B)] + cf [G2 (Eo, 2) + AGe, « (Eo, B)] + 
“=> + 63 m[Gam,s (Ep, P) thGam, s(Eo,P)] + --- = 0, 


ct [Gio (Eo, p) + AG, 9 (So, p)] + cs [G3,0(€, 2) +hGz3 o(&, Pia 


fe als C3 m-+1,0Gam+1,0 (Eo, pb) + h Gom+1,0(§0; p)| +---=0, 
co [G11 (Eo, b) +4 G11 (Eo. 6)] + oF [G3,1 (E>, B) +4Gs,1 (Eo, B)] + 
(60) vee Chm41 [Gomt1,1(€, p) +h Gom+1,1 (Eo; p) | +.--=0, 


of [Gis (Go. b) + 4G, 5 (Eo, P)) + oF [s,s (Eo, B) + Gs, .(E0,P)] + 
ag Be rary erskn Cae p) af h Gopi seo? p)| “ie = OY 


Applicando a tali sistemi di equazioni lineari il ,,principio delle ridotte“ di 
FourIER, gli autovalori del parametro # si determineranno per mezzo delle 
equazioni trascendenti che si ottengono uguagliando a zero i determinanti di 
ordine finito dei coefficienti delle incognite c; scaturiti dalla applicazione del 
principio suddetto. 

Supponendo poi ricavati tali autovalori, si considereranno le quattro cor- 
rispondenti successioni di autofunzioni Ve Va. eC Vi che costituiscono un 
sistema completo di funzioni di € e di 7 soddisfacenti a due a due alla condizione 
di ortogonalita. 

Riferiamoci ora alla seconda delle (35). Essa fornisce: 


Z(z) =C,cosgz+C,singz, 


con C, e C, costanti arbitrarie. Imponendo le condizioni ai limiti (36’), dopo 
aver posto, come é lecito fare, C,=1, si ottiene il sistema: 


(61) q—-hQ=0, GQ [h—gqtgql])+q+htgqi=o. 
Eliminando C,, si ottiene per il parametro g, l’equazione trascendente: 
(62) (h? — gq?) tg gl + 2hq=0, 


la quale ammette infinite radici; determinate tali radici, cioé gli infiniti valori 
di qg, si otterranno gli infiniti valori corrispondenti di C, dalla relazione C;=q/h. 
Si puo dunque ritenere risolto anche per il caso del moderatore a forma di 
cilindro ellittico, il problema della determinazione degli autovalori e delle auto- 
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funzioni dell’equazione diffusiva in questione con le volute condizioni ai limiti, 
hornet ere es ee 

e, in particolare, il problema della determinazione della ,,densita critica’. 

’ 


7. Evidentemente, nel caso semplice del cilindro a sezione circolare di raggio 
7, e di altezza /, riferendosi alle coordinate cilindriche 7, y, z e cercando della 
equazione differenziale (23) soluzioni della forma @ (7, p, 2)=ul(r, p) - Z(z), soddis- 
facenti alle condizioni ai limiti 


<0 +he=0, per 775; “2 _he=0 per “71; 


(63) 
<8 +ho=0, per 21, 


ci si ridurra allo studio del sistema differenziale: 


nu 14. du , 1 OFu oe 
(64) Ast) a PO) = sa aah gaage Mipeaagt a! lace 
ee AO Cd dart 


La seconda delle (64) dara luogo, con le sopra scritte condizioni ai limiti, ai 
risultati gia visti; la prima fornira autovalori del parametro # che sono le infinite 
radici della equazione trascendente (J,, funzione di BEssEL di prima specie di 
ordine 7) : 


(65) SEL hb) =O, per sf = Fy 


Dopo di che si scriveranno le autofunzioni per il caso considerato e si potra 
assegnare, in particolare, la densita neutronica critica. 


8. Va tenuto presente che lo studio matematico del funzionamento di un 
reattore atomico a fissione in regime transitorio si presenta assai pil complesso. 
Cid avviene soprattutto se non si mette esplicitamente in conto il contributo 
dei ,,neutroni ritardati“ cioé di quei neutroni che non vengono emessi istantanea- 
mente all’atto della fissione nucleare ma che vengono emessi con un certo ritardo 
per disintegrazione da parte di nuclei formatisi per fissione e che, prima di essere 
emessi, passano una certa fase della loro ,,vita‘‘ come ,,neutroni latenti‘‘. Invero, 
nella ,,vita‘‘ di un neutrone in una pila atomica di fissione, si delineano nettamente 
due fasi. La prima fase riguarda il tempo intercorrente fra l’emissione dei neutroni 
per scissione (neutroni veloci) ed il loro giungere alla soglia della energia termica 
(e cio vale tanto per neutroni istantanei quanto per neutroni ritardati). La 
seconda fase riguarda il tempo trascorso dai neutroni come termici fino alla 
loro cattura per dare luogo a nuove fissioni (€ il periodo di ,,vita‘‘ pit lungo, 
ma non il solo). Occorre dunque tenere conto del tempo finito di rallentamento 
dei neutroni veloci, tempo finito di rallentamento che perderebbe, invece, im- 
portanza qualora si considerasse, oltre al contributo dei neutroni istantanei, 
anche quello dei neutroni ritardati. Di tale tempo di rallentamento si tiene 
conto, come noto, scrivendo una equazione differenziale del rallentamento con 
diffusione dei neutroni veloci (‘age theory”), cioé una equazione che regge il 
comportamento dei neutroni nella pila e che vale fino ad energie immediatamente 
superiori al pit alto dei livelli di risonanza dell’uranio.? 


? Cir. Marsuak, R. E.: Reviews of modern Physics 19, 185, 1947. 
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Precisamente, considerato un mezzo moderatore ideale, in cui non avvengano 
fenomeni di cattura di neutroni, costituito da nuclei atomici di massa M e nel 
quale siano contenute sorgenti di neutroni, siano, rispettivamente: 


I. x(x, y, 2, #, t) la funzione ,,densita di rallentamento“ dei neutroni di energia 
E nel punto P(x, y, z), cioé il numero di neutroni che per unita di volume e di 
tempo raggiungono la age“ t definita come segue: 


E 
66 ea laeht aks Pa gsi i 
(66) r yeaa—ne | eee, 


rea 


dove # é l’angolo di ,,scattering nel sistema dell’osservatore, cos # ne é i 
valor medio del coseno, /(£) é il libero cammino medio di scattering dei neutroni 
di energia E ed Ey é l’energia dei neutroni veloci di fissione emessi (monocinetici), 
& é, infine, la perdita logaritmica media di energia di un neutrone per urto contro 
un nucleo del moderatore; 


II. S(x, y, z,¢, t) la funzione di distribuzione delle ,,sorgenti‘’ neutroniche 
riferita anch’essa all’unita di volume e di tempo; 


III. Hr) un coefficiente diffusivo definito dalla formula: 


(67) G(r) =v1(t)/3 (1 — cos #). 

Allora l’equazione di rallentamento con diffusione ed assenza di cattura dei 
neutroni veloci nel mezzo moderatore od equazione della ,,age theory“ é la 
seguente: 


O4(% 247) , 1 axlx, Wty " =A,y(x, y,2,t, t) + Sy(x, y, 2,4, 7), 
Ot " “G(t) al ( 
(68) , st HH 
OS pO 
a” 8x2 Tay? * azt 


Essa esprime, dunque, le vicissitudini dei neutroni aventi energia compresa fra 
la soglia di fissione dell’uranio e la soglia della zona termica. 

Le vicissitudini dei neutroni termici sono invece espresse dalla equazione di 
diffusione dei neutroni termici: 


(69) — LS S(x,y, 2,1) +D-A,0(%, y, zt) —o,v0(%, y,24), 
di cui ci siamo gia occupati e per la quale é noto il significato dei simboli. Come 
ho fatto vedere gia in un mio precedente lavoro’, una impostazione assai generale 
della teoria matematica del comportamento dei neutroni in una pila atomica 
di fissione si pud fare, allora, introducendo come ,,funzione delle sorgenti di 
neutroni termici’ nella (69) la soluzione, soddisfacente alle condizioni supple- 
mentari imposte dal problema, della equazione differenziale della ,,age theory”, 
e studiando poi la (69) con tale funzione delle sorgenti presa con t=T , valore 
di t corrispondente alla soglia termica e con le condizioni supplementari pure 
imposte dal problema. Si opereranno alcune semplificazioni. Anzitutto, essendo 


8 Cfr, PrGNEDOLI, A.: Rendiconti del Seminario matematico della Universita di 
Padova 1956. 
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relativamente breve il tempo di rallentamento dei neutroni veloci nel moderatore, 
si sostituira al coefficiente A(t) il suo valore medio &,, fra gli istanti corrispondenti 
all’emissione ed al raggiungimento della soglia termica. Inoltre alla funzione 
incognita y (x, y, 2, t, t) siimporra di soddisfare ad una certa condizione ,,iniziale“ 
rispetto alla ,,age‘‘ t, condizione consistente nell’ammettere che la distribuzione 
x(x, y, 2, t, 0) sia rappresentata da una funzione nota. Tale funzione, con ipotesi 
soddisfacente, potra essere assunta uguale a zero, poiché la densita di neutroni 
che, in ogni istante generico ¢ ed in un punto qualsiasi del moderatore, di co- 
ordinate cartesiane x, y, z, hanno ,,age‘‘ nulla appare ragionevole pensarla uguale 


a Zero. 


Dovra poi valere la condizione ai limiti: 
(70) oh +hy=0 (sul contorno del moderatore), 


dove é chiaro ed usuale il significato dei simboli. Operando sulla equazione 
differenziale (68) una trasformazione di Laplace semplice ed unilatera su inter- 
vallo infinito rispetto al parametro ,,age“ t, avremo, essendo k un numero 


positivo: 
us +00 
Spe ORME, WE Gat, 7) 1 -é =f 
gee DPE OSE os a fe te Nin es Pe dt 
i OT Dp ot a : 
Gi)aauae : 
+00 Pied 
=A, f e*" x(x, Vea ts adr + f OE SEV aelak) 2 Cs 
6 0 


cioé, essendo x(x, y, z,t, p) ed S,(x, y, 2,4, p) le rispettive L-trasformate di y e 
di S,; avremo l’equazione differenziale alle derivate parziali: 


s ‘ 1 0x (x, v, 4, 2, D) = 
(72) Any (*, y; A) t, p) Dn ot Dx (x, y; ay t, p)=—- Si (%, y; 2, b, p) 


alla quale andra associata la condizione al contorno: 
ox = 
ms +hy =0. 


Considereremo pure assegnata per y una condizione iniziale rispetto al tempo #, 
il che equivale ad assegnare una tale condizione per la funzione incognita 4. 
La (72), che é, sostanzialmente, la (2), la penseremo risolta per la regione occupata 
dal mezzo moderatore considerato e con le volute condizioni supplementari. 
Otterremo, quindi, la funzione y mediante L-antitrasformazione. 

Dopo cid, studieremo, come abbiamo detto, Vequazione di diffusione dei 
neutroni termici (69) ed il problema che ci eravamo proposto potra essere con- 
siderato risolto. 


Istituto matematico ,,S. Pincherle‘ 
dell’ Universita di Bologna 
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Stability in the Large of Systems of Two Equations 


IH; MUFTI 


Communicated by E. LEIMANIS 


1. Introduction 
Let us consider a system of linear equations 


dx, 


= hy % +My %2 


(1.1) 
dX, =5 hea xy a= No» Xo. 


It is well known that for the trivial solution of (1.1) to be asymptotically stable 
in the large the Routh-Hurwitz conditions 


yy lgg< 0,  Iyyhyg—Iyghgy > 0 (1.2) 


are necessary and sufficient. 

Here, we consider the problem of stability of the system (1.1) when h,, (7,7 = 
1, 2) are not necessarily all constants. 

In Section 2, we shall take /,,, 4, as continuous functions of x,, *, and in 
Section 3, #,, and h,, will be taken as continuous functions of x, and x, respec- 
tively. 

Finally, in Section 4, we correct the proofs of the three theorems proved by 
B. A. ErsHov [2] regarding the stability in the large of the system 


Es where o=¢,% — d, y. 


2. The stability in the large of <*. =h(y)x+ay, , =h.(x)x+by 


Let us consider the system 


(2.1) 


£Y —hpy(x)ax+by. 


We assume that the right-hand sides of the system (2.1) satisfy the conditions 
guaranteeing the existence and uniqueness of every solution. 


2a. Let us assume that 


Die K) =O 2 (OR S204 = OMiand) Ay(y)i=a0hs yylaso: 
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For the sake of definiteness we take a<0, hg(x)>0, »=+-0. We represent the 
right-hand sides of (2.1) on the (x, y) plane. The ordinates of the graphs of 


the curves are given by 


y= Fy (oi) x20 according as x [0 
a 3 (2.2) 
AG 


4 = = x20 according as x20. 

This shows that the curves /, (y) x-+-ay=0 and h,(x) x+by=0 he in the second, 
fourth and first, third quadrants respectively and consequently the origin is the 
only point of equilibrium. It is not difficult to see that 


ioe i(y)x+ay<0_ to the right of the curve h,(y)x+ay=0 
eS) xt+ay>o0 to the left of the curve h,(y)x+ay=0 
ou =h,(x)x+by<0 to the left of the curve h,(x)x+by=0 
cE ho(x)x+by>0_ to the right of the curve h,(x) x + by =0. 


The function y(t) attains its maxi- 
Y mum on the curve /,.(x)x+by=0 for 
y>O0O and its minimum for y<0; x (é) is 
minimized on the curve /,(y)x+ay=0 
for y>0O and maximized for y<0. The 
direction of motion and the curves 
obtained by putting the right-hand sides 
of (2.1) equal to zero are represented 
in Fig. 14. The curves and the co-ordinate 
axes divide the plane into eight regions. 
Using polar coordinates we see that 
the signs of 7 and @ in these regions are 


(14,5) * may be =0, p>0 
2,6) *<0, @ may be =O 
4 = 
(3,7) # may be =0, g>0 
(4,8) *<0, @ may be =0. 


Fig. 1 


The trivial solution is easily shown to be asymptotically stable in the sense 
of Lyapunov. In fact, let a V-function be defined by 


2V=2f «ho(x)dx —ay?. 
0 


Then : 
V = xhg(x) (Iy(y) x + ay) — ay(xh,(x) + by) 
= x*hy(y) hg(x) —aby?< 0, y += 0 
=0 possibly for y=0. 
Obviously, V is a positive-definite function. The derivative V is negative for 
y=+-0 and possibly zero for y=0. Since y=0 does not contain a positive half- 
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trajectory of the system (2.1) except *«=0=y, the trivial solution is asymp- 
totically stable according to Lyapunoy. Since b+, (y)<0, there are no periodic 
solutions according to the criterion of Bendixson. This is obvious from the 
fact that a Lyapunov function for the system is constructed. For the straight 
line L(0, cc) appearing in Theorem 1.3 (ERuGIN [3]) we take the positive half 
y-axis. We now show that the motions with bounded polar angles are bounded. 
We consider first the regions (2, 4, 6, 8). In these regions any motion with bounded 
polar angle is bounded, since 7<0. Next we consider the regions (1, 3, 5, 7). 
Here we show that there are no motions with bounded polar angles, 7.e., any 
motion starting or entering these regions must leave them after a sufficient 
time. For this we write 


> =—[— sin g (iy (y) x + ay) +008 o (p(x) x + by). 


Let us suppose that the motion started in the region (1) does not cross the 
curve h,(x)x+by=0, for, since x is decreasing and y is increasing in (1), y 
becomes infinitely large, and consequently ¢> «> 0, and hence with the increase 
of time the motion leaves the region (1). This contradicts our assumption. The 
contradiction shows that the motion must leave the region (1) and enter the 
region (2). The same reasoning holds for the region (5). Next we consider the 
rate of change of the quantity (bx—ay) along the trajectories of (2.1): 


d 
+ (bx - ay) =b(hy(y) x +ay) —a(hy(x) x+by) os 
= x(bh,(y) —ah,(x)) 20 according as x20. ‘ 
Consider the straight line 
bx —ay=A, (2.4) 


The straight line (2.4) intersects the curve xh,(y)+ay=0 for all A. Let A 
be negative; then the straight line (2.4) intersects the curve h,(y)*+ay=0 in 
the fourth quadrant. We consider the region bounded by the negative half 
y-axis, the straight line 6x —ay=A and the curve /,(y)«+-ay=0. The motion 
entering this region must cross the curve h,(y)x+ay=0, since it cannot cross 
the straight line bx —ay=A because of (2.3) and cannot go to the origin, since 
~> 0 in this region. A similar argument holds in region (3). The above analysis 
shows that all motions with bounded polar angles are bounded, and we have 
the following theorem: 


Theorem 2.1. If ahy(x)<0 for «+0, b<0 and h,(y)<0 for y+0, then the 
trivial solution of (2.1) 1s asymptotically stable in the large. 


2b. We shall now discuss the case when ah,(x)>0 and the conditions 
b+h(y)<0, bh(y) —ah.(x)>0, ee ee a0) (2.5) 


are satisfied. From (2.5) it follows that if ah,(«)>0, we must necessarily have 
b<0, h,(y)<0 for y=-0. We assume a> 0, h(x) >0, x==0. From (2.2) we have 


hy (vy h 
Mie = Spee i eee wa 
: (2.6) 
ee ae according as x20. 


Q* 
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Since y, and yy are positive or negative according as x is positive or negative, 
the two curves 
x% h(x) + by=0 
xh(y) +ay=0 
are situated in the first and third quadrants and are such that the curve 
«h,(y)+ay=0 is above the curve xh2(x)+by=0 in the first quadrant and 
below it in the third quadrant (because of (2.6)). We further note that 


(2,7) 


ae =h,(y)x+ay>0 tothe left of the curve %/,(y)+ay=0 
oe =h,(y)x+ay<0 totherightofthecurve xh,(y)+ay=0 
a =h,(x)x+by<0 to the left of the curve xh,.(x)+by=0 
@ =h,(x)x+by>0_ tothe right ofthecurve xh,(x)+by=0. 


The function x(f) is maximized on 
the curve x/,(y)+ay=0 for y>0O and 
minimized for y<0; y(t) is maximized 
on the curve *h,(x)+by=0 for y>0 
and minimized for y<0. The curves 
and the co-ordinate axes divide the (x, y) 
plane in to eight regions. The direction 
of motion and the curves are shown in 
Fig. 2. As before we introduce polar 
co-ordinates and notice that the signs of 
7 and @ in different regions are as follows: 

(1,5) * may be £0, ¢<0 

(2,6) *<0, P may be ZO 

(3,7) # may be £0, p>0 

(4,8) *<0, @ may be =O. 

In each of the regions (4) and (8) there will be at least one integral curve 
going to the origin, the only point of equilibrium. This follows from Theorem 2.4 
of EruaGin’s work [3]. Other motions started in the regions (4) and (8) either 
go to the origin or enter the regions (3,5) or (1,7) (since 7<0, @ may be =O 
in the regions (4,8)). We now show that motions entering the regions (1), (3), 
(5) and (7) must leave these regions and enter the regions (2) or (6). To show 
this we consider the rate of change of the quantity (bx—ay), 1.e., 


ss (6% —ay)= x (b hy(y) — a hy (x)) 20 according as x20. (2.8) 


Fig. 2 


The straight line bx —ay=A (A negative) has positive intercepts with the 
axes of co-ordinates and hence intersects the curves (2.7) in the first quadrant. 
Consider the region bounded by the straight line bx—ay=A (negative), the 
curve x/,(y)-+-ay=0 and the positive half y-axis. The motion M(é) entering 
this region cannot cross the line b«—ay=A (A negative) because of (2.8) and 
cannot go to the origin, since ¢<0, and hence must leave this region after a 
sufficient time. A similar argument holds for the regions (3), (5) and (7). The 
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motions after entering the regions (2) or (6) tend to the origin as t-+co. Thus 
we have proved the following theorem: 


Theorem 2.2. If ahy(x)>0 for x=+:0, then the trivial solution of (2.1) is 
asymptotically stable in the large under conditions (2.5). 


2c. We now assume that either b>0 or h,(y) >0 for yO and the conditions 
(2.5) are satisfied. For the sake of definiteness we assume that b>0. Then 
from (2.5) it follows that /,(y)<0 and ah,(x)<0. We suppose that a<0, 
hy(x)>0O for x=-0. The ordinates of the curves in (2.7) are 


y,=—%x LAA, according as x SO 
a 
hy (x) ee) 
ae? a 20 according as x SO. 


We compute the difference between the ordinates 

hy (y h(x 

| glee + % - ee 

x Daly) —4 hy (*) 

ab 

From (2.9) and (2.10) it follows that the two curves lie in the second and fourth 

quadrants, and the curve xh,(y)+av=0 lies above the curve xh,(x)+by=0 

in the fourth quadrant and below it in the second quadrant. It is not difficult 
to verify that 


(2.10) 


20 according as x20. 


=h,(y)xtay>0 _ to the left of the curve h,(y)x+ay=0 
ee =h,(y)x+ay<0 to the nght of the curve h,(y)*+ay=0 
A =h,({%)x-+ by=0 to the right of the curve h,(x) x + by=0 
2 =h,(x)x+by<0 tothe left of the curve A,(x)x+bdy=0. 


The function x(¢) is minimized on 
the curve x/,(y)+ay=0 for y>0 and 
maximized for y<0, and y(¢) 1s maxi- 
mized on the curve xh,(x)+by=0 for 
y>0 and minimized for y<0. The cur- 
ves (2.7) and the co-ordinate axes divide 
the plane into eight regions. The curves 
and the direction of motion are shown 
in Fig.3. The signs of 7 and @ in 
different regions are as follows: 


1,5) * may be 20, p>0 
2,6) *>0, p may be =0 
7) + may be =0) p>0 


) 


( 
( 
(3,7) 

(4,8) *<0, g may be =O. Fig. 3 


Let us compute @ in the regions (2) and (6) and see whether it is positive or not. 
Let c, be the greatest lower bound of the function /,(y). c¢ is finite, since 
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otherwise the condition bh, (y) — ah, (x)>0 for x0, y +0 Is violated. We then 
ee Shy) <5. 


From the inequality b/,(y) — ah, (x) >0 we have bc, — ah,(x)20. The equality 
sign is admitted if h,(y) does not attain its greatest lower bound; otherwise 
strict inequality holds. Let us suppose that c, is attained; then 


and Sih (2.41) 


If we let Ay(y)=c,+o,(y), then 0S a (y)< —(6+¢,) from (2.11). For he (x) 
we can take 


hy (4) = se + a(x) where a(x) >0. 
Now 
rp~=—xsing+ycos@p 
= —sing(xh,(y) +ay) + cosy (hg(x) x+ by) 
= — sin p[x (cy + (y)) + ay] + 08 | oa + a(x)) x+y]. 
Then 


p= — asin? y +" cost y + (6 — c,) sin g cos m — 
— o%(¥) sin y cos y + a(x) cos? p 


= — acos? @ [tan? y =| tan p ae | (2.12) 
— a (y) sin y cos gm + ay (x) cos? p (y + multiple of 2/2) 
= — acos? p(tan yg — 2) (tan g + - — a, (vy) sin pcos @ + a(x) cos? p. 


a 


The expression (tan p — =} (tan p +2) can change its sign only while 


passing through the values 2 and — 4. The value se > a , since 6+c¢,<0 
and a<0. Since ae — a, the straight line with slope R lies above the straight 
line with slope — <1 for x>0 and below for x<0. The straight lines y = = x 
and y= — 7 x lie in the second and fourth quadrants. We show that they do 
not lie in the regions (2) and (6). For this we have only to show that the straight 
line y= — 1x lies below the curve xh,(x)+-by=0 in the second quadrant 
and above in the fourth quadrant. In fact, this is so because 


byt pect 
a b 


deere Da ( = hy (x)) SO according as x SO. 
Thus we see from (2.12) that @ keeps the same sign throughout the regions (2) 
and (6) and which is easily seen to be positive. 

According to conditions (2.5) the origin is the only point of equilibrium, and 
there are no periodic motions. Let a motion M(é) start in region (1). The motion 
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M(t) must enter the region (2) after a sufficient time, since otherwise y becomes 
infinite, and then from Le =h,(y)x+ay it follows that the motion M(t) cannot 
remain in the first quadrant. We now show that it cannot remain in the region 
(2) for all time. For this we consider the straight lines 
bx—ay=A and ay+cx=0. 

A ¢,A 
b+c, ’ a(b+¢,) 
in the second or fourth quadrant according as A is positive or negative. We 
consider the region bounded by the straight line bx —ay=A (positive), the 


Their point of intersection is given by | ) This point lies 


positive half y-axis and the straight line ay+c,x=0. Since a (bx—ay)<0 


in the second quadrant and g>0 in (2), it follows that the motion enters the 
region (3). Here in the region (3), since g>O and the curve xh,(y)+ay=0 
lies above the straight line bx —ay=0O and y<0, the motion enters the region 
(4) after a sufficient time. The motion after entering the region (4) either tends 
to the origin or enters the region (5), since 7< 0 in region (4). Similar reasoning 
holds for the regions (5), (6), (7) and (8). The above analysis shows that any 
motion with bounded polar angle is bounded. For the straight line L(0, «) 
appearing in Theorem 1.3 (ERUGIN [3}) we can take the positive half y-axis. 
Now if we assume that the trivial solution is asymptotically stable in the sense 
of Lyapunov, then we have proved the following theorem: 


Theorem 2.3. If either b>0 or hy (vy) >0, y +0 and conditions (2.5) are satisfied, 
then the trivial solution of system (2.1) ts asymptotically stable in the large provided 
it 1s asymptotically stable in the small (i.e., according to Lyapunov). 

The requirement that the trivial solution be stable in the small can be 
realized, if we take, for example, 


ha (x) = 9% + a(x) = 9% +m + a9(x), 


where m > 0 and «3(x) does not contain any constant term and h, (y) =c,+% (y), 
where «, (y) does not contain any constant term. The equations of first approxima- 
tion can be written as 


dx 
aa aie ar Ey 
Cb AM aioe b 
dt ( a +m) fen 
The characteristic equation is 
A— ty a 
SPN yat PRESEN ys 
a 


or 


2 (b+c4)A4—am=0. 
The roots of this equation have negative real parts, since 
b+eaq<0 and —am>0, 


which ensures the stability in the small and which in turn ensures the stability 
in the large. 
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The origin is a node or a focus according as A= (b+ ¢)?+4amZ20. In case 
A~0, the directions along which the motions tend to the origin are given by 
au® + (cy — b)u pee m)=0. 


aa 


This is obtained by putting g=0 in the expression 


bc 
a 


+m) x+ by) 


—yvygp=sing(4*+7ay) — cos ¢(( 


and by writing uw = tan g, where tan y = = . The two directions are called critical 


directions. A critical direction is called singular if it satisfies the equation 
¢;+au=0; 


otherwise it is an ordinary critical direction. Along the ordinary critical direction 
enters either an infinite number of trajectories or only one integral curve. The 
question whether a finite number of trajectories or only one enters along a 
particular ordinary direction can be decided by using FROMMER’s criterion [4]. 


Along similar lines it is easy to prove the following theorem: 


Theorem 2.4. If b=0, then the trivial solution of (2.1) 1s asymptotically stable 
in the large under conditions (2.5). 


“4 = bax +Iy(x)y 


3. The stability in the large of Ss = rh (y)+ay, : 


Let us consider the system 


d 
Tp = ey) tay 
(3.1) 
Ug ae eM e 7e 
Oe + hy (x) ¥- 


We assume that the right-hand sides of (3.1) satisfy conditions which guarantee 
the existence and uniqueness of every solution. 


We first prove the following theorem: 

Theorem 3.1. 1f ab<0 and h,(y)<0 when y +0, and if hy(x)<0 when x0, 
the strict inequality being satisfied in at least in one of these conditions, then the 
trivial solution of (3.1) ts asymptotically stable in the large. 


Proof. We assume that a<0, b>0. A Lyapunov function for (3.1) under 
the conditions of the theorem is 


9 


2V (x, vy) =b x? — a y?. 


Its total time derivative in view of (3.1) is 


V=bdx(xiy(y) +ay) —ay(b% + hy(x) y) 
= b x*h(y)—ay*h.(x)\<0 for x+0, y+0 
=0 possibly on x=0 or y=O. 


Obviously, V(x, y) is an infinitely large positive definite function and x—0 
or y=0 does not contain any positive half trajectory of (3.1) except the origin. 
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Hence all the conditions of Theorem 4 (BARBASHIN [/]) are satisfied, which 
proves the above theorem. 


Next we consider the system (3.1) under conditions 


hy) +ha(x)<0, bly) ho(*%)—ab>0 for «+0, y=+0. (3.2) 


The following theorem can be proved in the same way as Theorem 2.2. 

Theorem 3.2. If ab>0, then the trivial solution of (3.1) 1s asymptotically 
stable in the large under conditions (3.2). 

The proof of Theorem 3.3 goes along similar lines as the proof of Theorem 2.3. 

Theorem 3.3. If either h,(y)>0 for y +0 or hy(x)>0 for x +0 and conditions 


(3.2) are satisfied, then the trivial solution of system (3.1) ts asymptotically stable 
in the large provided it is asymptotically stable in the small. 


4. The stability in the large of the system << =f (xy); 
dy =e 


is [(e) using qualitative methods 


Let us consider the system 


ax 


=F (x,y) 

dt 

Ls (4.1) 
Fe = I (0) »-where ¢ = ¢, % —d,' 


and c,, d, are positive constants. We shall assume that I(x, y) is a continuous 
function, having first order partial derivatives with respect to x and y, for all 
values of x, y. Further we assume that 


see aie AO. 0) tO (4.2) 
: RO, LOLS Sry ode 
We shall discuss three cases (1) - = 0, (i) ==), Amt) a = 0. 


The continuous function f(c), appearing in the right-hand side of the second 
equation of (4.1), is subjected to the following conditions: 


Gf) =0 fora ==0 (4.3) 
{(O)j=0 (4.4) 
mie = 0pon alka. (4.5) 
(xe) 
In the condition (4.5) a is not taken to be identically zero. The equations (4.1) 
can be written as 
dx OF OF 
dt ie (Se) y=0 an ay 4 15 a x (%, y) 
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where eal 
ok oF 
1%) =Flx 9) —(F) _, o* — (By lec, yao” 


w(x, vy) =f(o) an sO ey” a ae: 


For simplicity we write 


Seas lew 2) lott 


The system (4.6) can then be written in the form 


Aen ai 


00 age y=0 


a —Nax—by+y4(%,¥) 
(4.7) 
& =cx—dy+y(x,9). 
The three cases thus correspond to the values 1, 0, —1 of N in (4.7) respectively. 
Gase Li = <0. From equations (4.7) we have for N= 1 
x 
f _ _ax—by+z(x,9) 
dt 
ty (4.8) 
Wp ot 4Y Tle Y)- 
The equations of the first approximation are 
ve =—ax—by 
| (4.9) 
albest cx—dy 
dt : 
The characteristic equation of (4.9) is 
A —0b 
+a =, 
C At+d 
or (4.10) 


AM+(atdAtad+bc=0. 
The roots of this equation have negative real parts, since 


at+td>0O and ad+bc>0O. 


Since the roots of (4.10) have negative real parts, it follows that the trivial 
solution x =0=y of (4.8) is asymptotically stable according to Lyapunov. The 
absence of periodic solutions is easily seen by using the criterion of Bendixson. 
In fact 


Ba Op ine Shes ea 
ax * dy On 1 G0 
in view of condition (4.5) and the fact that a = O70: 
x 
We represent the curves F(x, y)=0 and /(c)=0 on the (x, y) plane. By 
virtue of conditions (4.3) and (4.4) f(o)=0 represents the straight line o=0. 
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Since c, and d, are positive constants, o=0 is situated in the first and third 
quadrants. The curve F(x, y)=0 passes through the origin since F(0, 0)=0. 
It is situated in the second and fourth 
quadrants since the slope of F(x, y)=0, 
ee. 


id hoed EGE is negative 
dx —OF/éy”’ shee ts 


These considerations show that (0, 0) 
is the only point of equilibrium of (4.8). 
We have shown so far that (i) the origin 
is the only point of equilibrium, (ii) it 
is asymptotically stable in the sense of 
Lyapunov, (iii) there exist no periodic 
solutions. We shall now show that there 
exists a straight line L(0, co) which is 
intersected by the motions in one direc- 
tion only, and all motions with bounded 
polar angles are bounded. For this purpose we examine the directions of motions 
(see Fig. 4). It is easy to see that 


Fig. 4 


ae = F(x, y)>0 for points lying to the left of the curve F(x, y) =0 


a =F (x, y)<0 for points lying to the right of the curve F(x, y) =0 
a =f (oc) 0 for points lying below the straight ine o=0O 


oe (Oi) <i) for points lying above the straight line o=0. 


The function y (¢) attains its maximum for y> 0 andits minimum for y< 0 on the 
straight line f (oc) =0; x(t) attains its maximum for «> 0 and its minimum for *< 0 
on the curve F(x, y)=0. The straight line f(c)=0, the curve F(x, y)=0 and 
the coordinate axes divide the (x, y) plane into eight regions. We introduce 
polar coordinates 

L==7 COS, == 7S Oe 
Then 
7=xcosy+ysiny 
and 
rp=—xsing+ycos@q; 
the signs of 7 and @ in different regions are as follows: 

(1,5) 7 may be =0, ¢>0 

(2,6) *<0, p may be =0 

(3,7) 7 may be =0, ¢>0 

(4,8) 7<0, @ may be =O. 


Ersuov [2] argued that since y>0 in the regions (1), (3), (5) and (7) there 
cannot be any motion with bounded polar angle w in these regions, and any 
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motion falling into these regions or starting there has to get out of them after 
intersecting either the straight line /(¢)=0 or the curve F(x, y)=0. To us this 
reasoning is doubtful, since p> 0 is not sufficient to guarantee that there cannot 
be any motion with bounded polar angle in these regions. However the above 
assertion would remain true if we could show that ¢>e>0 in these regions. In 
fact, if @>e, then integrating we have p—q@>e(t—t), whence it follows 
that as ¢ increases, p increases, and hence there will be an instant of time when 
the motion leaves these regions. ; 
Let a motion M (i) after intersecting the negative half y-axis enter the region 
(7). We show that motion M(/) cannot remain in this region, and consequently 
it will enter the region (8). To show this we write 
gas [— xsing +¥ cos ] 
‘ (4.11) 
== ‘ (— F(x, y) sing + f(a) cosy) > 0 in region (7). 
If M(t) does not cross the curve I(x, y)=0, then x(#) becomes infinite since 
x(t) and y(é) are increasing and p>0. For sufficiently large x, with the increase 
of 
éo 
therefore from (4.11) it follows that g can always be taken greater than ¢«>0, 
and hence after a sufficient time the motion must leave the region (7), which 
contradicts our assumption. After entering the region (8) it either goes to the 
origin as ¢—>-++ co or intersects the x-axis and enter the region (1). This follows 
from the fact in this region *<0 and @ may be =0. The motion cannot remain 
in the region (1) and must enter the region (2). To see this, consider the region 
bounded by the straight lines c=0, x=A and the x-axis. Since in the region 
(1) x decreases, the motion cannot intersect the line x =A. It cannot go to 
the origin, since p> 0, and therefore it must necessarily go out of the region (1) 
and enter the region (2). Here, 7.e., in region (2), since *<0 the motion either 
goes to the point of equilibrium or enters the region (3). Similar reasoning 
holds for the rest of the regions. All this shows that any motion with bounded 
polar angle is bounded. For the straight line L(0, «) we can take the positive 
half y-axis. Thus all the conditions of Theorem 1.3 (ERUGIN [3]) are satisfied, 
and we have the following theorem: 


of x, o increases. Since it is assumed that = 0, f(a) is non-decreasing, and 


Theorem 4.1. Let <0. Then under conditions (4.2)—(4.5) the trivial 
solution of (4.1) is asynrptotically stable in the large. 
OF : : : 
Case II, Frees! The equations (4.7) in this case take the form 
dx 
Fe OD 
a (4.12) 
dp = o%— ay + p(x, 9). 
The equations of first approximation are 
ax 
Lie * oy 
Lier ae (4.13) 
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The characteristic equation of (4.13) is given by 


A —b Br 
ey ea, 51 

or 
A+di+bc=0. 


The real parts of the roots of this equation are negative, since d>0 and bc>0, 
whence follows the asymptotic stability of the trivial solution of system (4.12) 
in the sense of Lyapunov. 

The curve F(y)=0 represents the straight line y=0, and f(a) =0 represents 
the straight line c=0. The origin is the only common point of f(¢)=0, and 
F(y)=0. It is easy to see that 


—- =F(y)>0_ for points lying below the x-axis 
=F (y)<0_ for points lying above the x-axis 
oe =f(c)>0 _ for points below the straight line /(c) = 0 
oe =f(c)<0 for points above the straight line {(¢) =0. 


The function y(¢) is maximized for 
y>0O and minimized for y<0O for the 
points on f(c)=0, and x (¢) is maximized 
for *>0O and minimized for x<0 for 
the points on the w-axis. The curves 
F(y) =0, f(¢)=0 and the axes of coordi- 
nates divide the plane into six regions. 
The direction of motion is represented 
in Fig.5. The signs of 7 and @ in different 
regions are given as below: 


(1,4) 7 may be =0, p>0 
(2,5) 7<0, g may be =0 


(3,6) * may be £0, >0. 


In the region (6) it is easy to see that a =>1> 0, whence it follows that 


there exists an instant ¢ when the motion intersects the x-axis. As in Case I 
it can be shown that the motion M(é) intersects the straight line /(¢)=0 and 
enters the region (2) after a sufficient time. In region (2), after a sufficient time, 
the motion either goes to the origin or enters the region (3). Similar arguments 
can be applied for the regions (3), (4) and (5). Thus we have shown that (i) (0, 0) 
is the only point of equilibrium, (ii) it is asymptotically stable in the small, 
(iii) any motion with bounded polar angle is bounded. To show that there exist 
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no periodic motions, we use the criterion of Bendixson 
Ge + oat dy, - oe (not identically equal to zero). 
00 


For the straight line L(0, cc) we can take the positive half x-axis. Thus all 
the conditions of Theorem 1.3 (ERUGIN [3]) are satisfied, and we have the follow- 


ing theorem: 


Theorem 4.2. Let oF _. Then under conditions (4.2) —(4.5) the trovial 


Ox 
solution of system (4.1) is asymptotically stable im the large. 
Caselil. <-> .0. The equations (4.7) in this case take the form 
x 
ax —by + 4(%9) ; 
d (4.14) 
ees s . 

a OF dy+y(x, ). 

The characteristic equation of the first approximation is 
Mere, Sea al: ‘ 
c At+d 

oe (4.15) 


M+A(d—a)+bcec—ad=0. 
The roots of (4.15) will have negative real parts under the following conditions: 
d—a>0O, (4.16) 
be—da>0O. (4.17) 
If the conditions (4.16) and (4.17) are satisfied, then the trivial solution 
x=O=y of (4.14) is asymptotically stable in the small. We now represent the 


curves f(x, y)=0 and the straight line /(c) =0 on the (x, y) plane. Since 7 > 0; 


the curve F(x, y)=0 is situated in the first and third quadrants and so is the 
straight line given by f(¢)=0. Since we are interested in having a unique point 
of equilibrium, we shall have to impose an extra condition: 

Condition 1. The curve F(x, y)=0 is situated between the x-axis and the 
straight line /(o)=0. 

Obviously, Condition 1 implies the condition (4.17). It is not difficult to 
see that 


ii a F(x, y) > 0 for the points lying to the right of the curve F(x, y) =0 


= = F(x, y)<0 for the points lying to the left of the curve F(x, y) =0 
d ; 

= = f(o)>0 for the points below the straight line fia) = 
dy 


oF Fes f(o) <0 for the points above the straight line HiG\== 103 
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The function y(¢) is maximized for y>0 and minimized for y<0 for the 
points on the straight line f(¢)=0, and x(¢) attains its maximum for x>0 and 
its minimum for x<0 for the points lying on the curve F(x, y)=0; the direction 
of motion is shown in Fig. 6. The curve F(x, y)=0, the straight line #(c)=0 
and the axes of co-ordinates divide the 
plane (x, y) into eight regions. As before 
we introduce polar co-ordinates 


X%=rcosp, Yy=rsing. 
Then 
7=—xcosy+ysing 


1p=—xsing+ycos@. 
The signs of 7 and @ in different regions 
are given as 
(1,5) *>0, @ may be £0 
(2,6) 7 may be £0, p>0 
(3,7) 7<0, » may be == (0 
(4,8) * may be =0, g>0. 


Fig. 6 


In the regions (1) and (5) we see that 7>0 and @ may be —0. First of all 
we need a condition that makes @ positive in these regions. This is necessary 
in order to ensure that there is no motion with bounded polar angle in the regions 
(1,5). Our purpose is served if for example we assume that 

WAG Neo), 1X, V) Ow in the. region. (4) 


and ; (4.18) 
U(x, V) 0, wary) = 0” inthe’ region (5). 


Next we require that there be no periodic motions, and for that we must have 
OF of 
we ell eS 
Ox a 00 =0 (49) 

according to the Bendixson criterion. 

Let a motion M(é) enter the region (8) after intersecting the negative half 
y-axis. It is easy to see as in Case 2 that the motion M(#) crosses the x-axis 
and enters the region (1). We now show that it cannot remain in the region (1) 
and must enter the region (2). This is done as follows. We consider the rate 
of change of (cx —ay) along the motion M (8), 1.e., 

a 
dt 


(cx —ay)=acx—bcy+cy(x, y)—acx+ady—ay(x, y) 
= y(ad — be) +cyx(x,y) —ay(x, y) <0 
in the region (1) in view of (4.17) and (4.18). 
We consider the straight lines cx —ay=A (A positive) and cx —dy=0O and 


find their point of intersection foes : 


A —— : 
Ta ee ): It lies in the first quadrant, 
since A, c, d, d—a are all positive quantities. Since F(x, y)=0 lies always 
below the straight line /(o)=0, the straight line cx—ay=A _ intersects the 
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curve F(x, y)=0 for all 4. Thus the motion M(é) entering the region (1) must 
cross the curve F(x, y)=0, because it cannot go to the origin since p>0 and 
cannot cross the straight line cx —ay=A due to (4.20). After intersecting the 
curve F(x, y)=0, it cannot even remain in the region (2) because of the same 
reason. Thus the motion enters the region (3), where it either tends to the origin 
with the increase of time or goes out of this region and enters the region (4). 
Similar argument holds for the regions (4), (5), (6) and (7). This shows that 
there cannot be any motion with bounded polar angle in the regions (1), (2), 
(4), (5), (6) and (8). The motion with bounded polar angle can only occur in 
regions (3) and (7) where it is bounded, since 7< 0. Thus any motion with bounded 
polar angle is bounded. For the straight line L(0, cc) we can take the positive 
half x-axis. Hence all the requirements of Theorem 1.3 (ERUGIN [3]) are satisfied, 
and we have the following theorem: 


Theorem 4.3. Let - >0. Then under conditions (4.2)—(4.5), (4.16), (4.18), 


(4.19) and Condition 1, the trivial solution x =0=y of system (4.1) 1s asymptotically 
stable in the large. 
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Feblerabschatzungen aus den Defekten 
bei Randwertaufgaben gewohnlicher und partieller 
Differentialgleichungen 2. Ordnung 


Hoa GLoLsteHN 


Vorgelegt von L. COLLATZ 


In einer friiher erschienenen Arbeit hat J. ScHRODER [6] Abschadtzungen fiir 
die Lésung w der Differentialgleichung M[w]=/(x, w) mit gewissen Rand- oder 
Anfangsbedingungen angegeben, wobei M einen linearen homogenen Differential- 
ausdruck darstellt. Die Beweise wurden dabei mit Hilfe der Existenz einer nicht- 
negativen Greenschen Funktion durchgefiihrt. In dieser Arbeit sollen allgemeiner 
fiir in w und den Ableitungen 1. Ordnung nichtlineare Differentialgleichungen 
Fehlerabschatzungen hergeleitet werden. Die Voraussetzungen sind hierbei andere 
als bei SCHRODER, speziell wird z.B. bei den Beweisen nicht die Greensche Funk- 
tion herangezogen. 


§ 1. Bezeichnungen und Voraussetzungen 
Fir die Differentialgleichung 


an Tu=—4;,uU,,+F(%,u4,u4,)=r(x) in B, 
Ae 

U,u=u—a ot =}, auf mit «20 
sollen Abschatzungen fiir «—v, wobei v(x) eine beliebige Naherungslésung fiir 
u(x) bedeutet, aus den Defekten auf J’ und in B angegeben werden. »% repra- 
sentiere die m Verdnderlichen (*,, %,,..., %,), B sei ein beschranktes offenes 
Gebiet des n-dimensionalen (x)-Raumes, J’ dessen Rand. Die Indizes bei Funk- 


tionen w(x), v(x) usw. sollen partielle Ableitungen nach der betreffenden x-Ko- 
ordinate bedeuten, also 4;= » Uj R= oe . Bei den a;,(x) dienen die ange- 
hangten Indizes dagegen nur zur Unterscheidung der Funktionen @,,, a9, .--, @y- 
In (1.4) wurde von der Summenkonvention Gebrauch gemacht, d.h. es wurde 
iiber zwei in einem Produkt auftretende gleiche Indizes von 1 bis m summiert. 
Die Funktionen a,,(x) seien in B stetig, die symmetrische Matrix a;,=4a,,; sei 
dort positiv definit. @u/do ist die Ableitung auf dem Rande J’ in Richtung der 
Konormalen, die unter den Voraussetzungen tiber die a;, stets in das Innere 
von B weist (s. Coriarz [J], S. 23). S sei ein hinsichtlich « und (w,) konvexer, 
abgeschlossener Bereich des (%, u, u;)-Raumes mit (x) ¢ B+". 
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v(x) €K soll heiBen: v besitze in B+ I'stetige partielle Ableitungen 1. Ordnung, 
in B stetige partielle Ableitungen 2. Ordnung einschlieBlich und der Punkt 
(x, v(x), v;(x)) liege in S. 

Die Funktion F(x, uv, w;) besitze in S stetige partielle Ableitungen nach w 
und den w,. ! E 

Es werde vorausgesetzt, daB die Randwertaufgabe (1.1) eine Lésung w(x) CA 
besitze. Auf Existenzfragen wird also nicht eingegangen. Die Eindeutigkeit der 


Lésung von (1.1) wird dagegen nicht verlangt. 
Zu einer beliebigen Naherung v(x)€K existiere eine (2mal stetig differen- 


zierbare) Funktion g(x) mit den Eigenschaften 
peo und Lg =— aj, 94+ Fy(x,0,0) 9+ F,(%2,0)-9>0 in B, 
Ure 0 aut: 


I 


(1.2) 


Fiir eine derart festgewahlte Funktion (x) erfiille die Funktion 


(1.3) G(x, 2)=F(x,v+@2, v,+9;2)—F(x, v, v;) -E, (%, 2, 0;) - pj 2—F,,(%, v, 04) - pz 


die Bedingung 
(1.4) NG) zs G (42) = Naya it a ae 


Voraussetzung in (1.3) ist selbstverstandlich, daB mit (x, v, v;) auch (x,v+@z, 
v,+ 9,2) in S liegen soll. 


§ 2. Herleitung einer Abschatzungsformel 
Fiir eine Naherungslésung v(x) € K der Lésung w(x) von (1.1) bilden wir eine 


7) 


Fehlerfunktion z= “—" mit einer festgewahlten Funktion g(x), die den Bedin- 


We 
gungen (1.2) geniigt, und definieren den linearen Differentialoperator 


Mit (¥u—v),=p2,+9,;2 und (4@—v);,= 9 2;,+ (My 2%; +Y; %) + Vj, 2 erhalten wir 
nach einer kee Rechnung 


~ 


Lz=Tu—Tv—([F(4,0+92,0,+ 9,2 


1 


| 
T 
= Fi, (wiv, ¥,) pee =F ,@, 9) 2])) mB! 


Uj 


(222) Y 2;) F(x, OE v,) ee 


Auf dem Rande /’ wird 
Oz 
(2.3) U,p-2—ag <2 =U, (u—2). 


Fiir z(x) sollen jetzt in B+J” Schranken hergeleitet werden. Sei 7= Mae das 
in einem Punkte 2)= (x9) angenommen werden mége. 
Fir cl folgt aus dem Taylorschen Satz (02/0c)pS0, also wegen ap=O 
auf J” aus (2.3) 
U,,(u—v) U,,(u—v) 
= << xe es ———— = 
: s(“* U,@ ), = Max U.P aa. 
Fir (cB besitzt z in P, ein differenzierbares Maximum mit z, j(%)=0 und 
(4;4 2)n)p,S0 (s. Cortatz [1], S. 374). Aus den Darstellungen (2.1) und (2.2) 
und mit der in (1.3) eingefithrten Funktion G(x, z) folgt unter Beriicksichtigung 
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der Bedingung (1.4) somit 


iS hey Gia 2) S 
<< ye 2 
z= ( LG i ( oj Sot Mi? 


mit den Abkiirzungen 


Lady 
0» = sup ee 


und M, = sup Fa 20. 


2 geniigt also der Bedingung 
(2.4) g(2) = M, 22-74 6,20. 


1. Fall. Ist M,=0, so folgt aus (2.4) unmittelbar 7<6,. 
2. Fall. Ist M,>0, so kénnen die in Fig. 1a oder 1b fiir 6,.>0 bzw. 6,.<0 
dargestellten Méglichkeiten eintreten. 
g Y 


4 06,/7,>7 406,/4,<7 


Fig. 1a. (d2>0) Fig. 1b. (6.<0) 


Im Falle 6.<0 erhalten wir stets eine Schranke fiir Z, die negativ ausfallt: 


TS a (1 yi—46, M,). Ist dagegen 6,;>0, so kénnen wir 2 nur dann ab- 
1 


schatzen, wenn die durch g(2) dargestellte Kurve die z-Achse schneidet. Die 

Bedingung hierfiir lautet 46,M,<1, was stets zu erreichen sein diirfte, wenn 

die Naherung v(x) ,,nahe genug“ bei w(x) hegt. Dann wird wie oben 

C= ar (i! Vi 40) ,M,). In diesen Schranken ist der Grenzfall 6,=0 mit 
1 

2<0 enthalten. 


In ganz entsprechender Weise kann Min z abgeschatzt werden. Die vorlie- 
+r 


genden Ergebnisse stellen wir zusammen im 


Satz 1. Sei u(x) die Lésung der Randwertaufgabe (1.1), v(x) CK ene Naherung 
fiir u(x) und p(x) eine Funktion, die den Bedingungen (1.2) geniigt. Erfiillt die 
durch (1.3) definierte Funktion G(x, z) die Ungleichungen (1.4), so galt mat den 
A bkiirzungen 


d, sa Min U,(u— : Oy as inf ‘le — = IPG : lV Be — sup — aw N,(¥). = 0, 
(2.5) if ne L* = @P 
U,,(u—v) Tu—Tv : ss) 
Ss {2 a25 zee i <a 
dy Max Ute 0» SUD eras ta M, inf — 7 =30., 


in B-+I" die Fehlerabschatzung: 
im Falle M,.M,<0, 46,M,<1 und 46,M,<1 
: 1 (ee SE pe a 
x) Min |dy; sy (1 1—40,M,) 
p ( ) i 2M, ( V 1 ) 


(2.6) 
<u—v<sq(x)- Max 


1 
dy ap (VAs 1—40,M) 
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und im Falle M,=M,=0 
(2.7) (x) - Min [d,; 6,] Su —vS p(x) - Max[d,; 6,]. 

Gilt U,(w—v) <0 auf I’ und Tu — TvS0 in B, so folgt aus dem obigen Satz 
der in [3] angefiihrte Monotoniesatz, d.h. w—vS0 in B+’, hier allerdings 
unter etwas anderen Voraussetzungen iiber die Funktion G(x, z) als friher. Fur 
eine auf irgendeine Art (z.B. durch Kollokation, Ritzsches Verfahren, Fehler- 
quadratmethode usw.) bestimmte Naherung v (2) werden jedoch im allgemeinen 
die Monotonievoraussetzungen auf J’ und in B nicht erfiillt sein. In diesem 
Falle kénnen die im Satz1 angegebenen Fehlerabschatzungen herangezogen 
werden. In welcher Weise dieses geschehen kann, soll im folgenden ausgefiihrt 
werden. 

§ 3. Untersuchung einiger Sonderfalle 

Bei der Aufstellung einer Naherungslosung v(x) fiir die Lésung u(x) von 
(1.1) wird man von v die Erfiillung gewisser Bedingungen verlangen. Bei partiellen 
(linearen) Differentialgleichungen ist es zuweilen méglich, v(x) so zu wahlen, daB 
die Differentialgleichung exakt erfiillt ist. LaBt sich v(x) nicht so wahlen, z.B. 
bei gewohnlichen und partiellen nichtlinearen (in den meisten Fallen auch bereits 
bei den linearen) Differentialgleichungen, so wird man von v(x) méglichst die 
Erfiillung der Randbedingungen verlangen. 


a) v(x) erfiillt die Differentialglerchung exakt 
Ist Tu=Tv in B, so wird 6,=6,=0 und wir erhalten den 
Satz 2. Evfiillt v(x) die Differentialgleichung (1.1) exakt, nimmt der Defekt 


U,(u—v) auf I’den Wert Null an und existiert eine den Bedingungen (1.2) gentigende 
Funktion (x), so gut in B+ die Fehlerabschatzung 


- Uy(u—v U, (u—v 
(3.1) (2) - Min #7") < u(x) — 0(x) Sox) - Max) 


Daber ist p(x) eine fret wahlbare Funktion mit p>0 und L* y=0 in B, U, p>0 
auf I’. — Ist U,(w—v) auf I 2.B. nur positiv, so ist die linke Schranke in (3.1) 
durch Null zu ersetzen. 


Der Beweis ergibt sich sofort aus Satz 1, es mu8 nur noch gezeigt werden, 
daB L*py=0 in B zulassig ist, wahrend frither L*g>0 in B gefordert wurde. 
Dazu brauchen wir nur $=p-+eg mit e>0 zu setzen, dann gilt @>0 und 
L*~>0 in B, U,~>0 auf I, dh. @ erfiillt die Voraussetzungen des Satzes 1. 
Fiir e—0 folgt daraus (3.1). 


b) v(x) erfiillt die Randbedingungen exakt 
Satz 3. Ist v(x) eine die Randbedingungen exakt erfiillende Naherung fiir u (a 
so gelten die Abschiitzungen (2.6) bzw. (2.7) mit d,=d,=0 und einer den Bedin- 
gungen p>0 und L*p=p>0 in B, U, y= auf I geniigenden Funktion (x). 


Zum Beweise setzen wir G=g-+e mit e>0 und ¢ so klein, daB E* es 
L*op+e - F(x, v, v;)>0 in B wird, was wegen F, beschrankt in S stets moglich 
ist. Dann erfiillt p mit U,~=e>0 auf I" die Voraussetzungen des Satzes 1 und 
mit ¢—0 ergibt sich der Satz 3. 
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c) Bemerkungen zur praktischen Durchfiihrung der Abschdtzung 
Wir wollen uns im folgenden auf den Fall 


(3.2) F(x, u, u;) = — b;(x)-u;+ f(x, u) 

beschranken (bei den 6; dienen die angehangten Indizes nur zur Unterscheidung 
der Funktionen 0,, b,;...,6,, sie bedeuten hier keine partiellen Ableitungen). 
Dann wird 


G(x, z) =f(x,0 + pz) — f(x, 2) — f,(%, 0) + pz. 
Setzen wir von fin S die Existenz und Beschranktheit der partiellen Ableitung /, , 


voraus, so wird G(x, z)=$f,,,(%, ) - (py z)2, wobei 3 eine Zwischenstelle zwischen 
v und v+ @z bedeutet, und wir kénnen abschiatzen: 


eee oN FONG é Gn oe 
(3.3) Ni (x) =->-- Min [0; Min RA Ca < - Max [0; Max ale 
Ist _ 1, (x, uw) 2c(x), so kénnen wir g(x) >0 in B aus . 
L*¥ gp ZL** p= — 4;,9;,— 5,9; + op >0 
(bzw. 20 im Fallea) in B, U,y>0 (bzw. U,g 20 im Falle b) 


7 


(3.4) 


auf J” bestimmen. Dann brauchen wir nicht fiir jede Naherungsfunktion v (x) 
erneut L*q@ zu berechnen, sondern kénnen fiir alle v(x) mit dem festen L** 
rechnen. Allerdings tritt hierdurch im allgemeinen eine Vergréberung der Ab- 
schatzung ein. 

Bei gewohnlichen Differentialgleichungen wird man meistens von Naherungs- 
lo6sungen v(x) ausgehen, die die Randbedingungen bereits erfiillen. Ist die Diffe- 
rentialgleichung im Intervall a<x<b gegeben, so fiihrt haufig der Ansatz 
y (x) =(* —a) (b— x) zum Ziel. Aber auch andere Funktionen sind méglich und 
k6nnen unter Umstanden bessere Schranken liefern. — Bei partiellen Differential- 
gleichungen hangt die Wahl von @ sehr stark vom Bereich B+ TI’ ab, allgemein- 
giiltige Regeln lassen sich hier kaum angeben. Liegt B+ J" z.B. ganz in einem 
Quader |x;— %);|Sa; (¢=1, 2,..., ”), so liegt ein Ansatz p=1—A,(x;— %o,)? 
mit passend gewahlten Konstanten A, nahe. 

Mit einer festgewahlten Funktion g(x) lassen sich die im Satz 1 angegebenen 
Konstanten d, 6 und M bestimmen. Wahrend man d und 6 moglichst genau 
berechnen oder abschatzen sollte, geniigt fiir M (und damit auch fiir N) haufig 
eine grdébere Abschatzung. Ist namlich 46M<1, so wird die Schranke 
oar (1— V1—46 M) ~ 6(1+6M) verhaltnismaBig unempfindlich gegentiber Ande- 
rungen von M. 

Es sei noch bemerkt, daB sich die Abschatzungen (2.6) bzw. (2.7) auf beliebige 
Iterationsverfahren (wobei besonders das Newtonsche Verfahren erwahnt sei) 
anwenden lassen. In diesen Fallen kénnen die Gebietsdefekte 6 in eine fiir 
numerische Rechnungen zweckmaBigere Form gebracht werden, doch soll darauf 
hier nicht eingegangen werden. 


§ 4. Beispiele 
Beispiel 1. Nichtlineare gewohnliche Differentialgleichung 


(44) Tu=—u"’—(1+%) (1+u')u=2 fiir 0<x<1, u(0)=u(1)=0 ( a): 
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Wir fragen nach nichtnegativen Lésungen u(x) und konstruieren zunachst 
eine Naherungslésung v(x), die die Randbedingungen bereits erfiillt. Fiir 
v(x) =x(1—x)(@+bx)+c- sin zx mit Konstanten a, b, c wird 


Tu—Tv=2(1—a+)) —6bx 
—g?c-sinax+ (1+ x)(1+’)v. 


Die Forderung Tu — Tv=0 fir +=0; 331 
liefert die Werte a,=1, b,=0, c,=0,0395 
und fiir Tu—Tv etwa die in Fig. 2 ge- 
strichelt dargestellte Kurve. Daraufhin 
wurden a, b,c erneut aus der Forderung 
Tu— Tv=— 0,09 fir +=0, =O fir se 
— 0,09 fiir x1 errechnet. Auf diese Weise 
wurden 


v(x) = «(1 — x) (1,015 — 0,03 x) + 0,043 88 sin x x 


Ju-v 


Fig. 2 


und 
Tu —Tv=— 0,09 + 0,18 x — 0,433 08sina x + (1+ x) (1+ ')v 


(in Fig. 2 die ausgezogene Kurve) erhalten, deren Werte in der Tabelle in Spalte 2 
bzw. 4 stehen. Mit y(x)=x(1—~) wird 
L*yp=2— (14+ x) [v(1— 2x) + (1+) x(1—x)] (Spalte 5). 


In Spalte 6 der Tabelle stehen schlieBlich die Werte ite Te 


und Minimum fiir die Fehlerabschatzung ben6tigt werden. 


, deren Maximum 


* 


Tabelle 1 
oe v (x) v(x) Tu—Tov L* | co ae 
L*g 
O O 1,1529 = 009005 4)9 <2 —0,045 
0,1 0,1046 0,9380 0,0172 1,716 0,010 
0,2 0,1873 Owed O;O763i ie e536 0,050 
0,3 0,2468 0,4771 0,0875 1,475 0,059 
0,4 0,2825 | 0,2360 0,0589 1,506 0,039 
0,5 0,2939 = O0075 — 0,0044 O23" = 0,008 
0,6 0,2810 — 0,2492 —0,0563 1,802 | —0,031 
0,7 0,2442 —0,4849 —0;1005 1,982 =O O54 
0,8 0,1844 =O 7/1OS: —0,1056 2,116 —0,050 
0,9 | 0,1015 OOo 2, —0,0488 2135 a OLO22 
1 O 12295) 0,0900 2 0,045 
Mit Fy=—(1+x)u und F,=—(1+.%)(1+w’) wird nach (1.3) G(x, 2)= 


—(1+%) py’ 2=— x(1—**) (12x) 22 und hiermit kénnen wir abschatzen: 


M,<0,085 und M,>—0,075. Mit —0,054< a <0,062 erhalten wir aus 
(2.6) die Fehlerabschatzung: v 


— 0,0542% (41 — x) S u(x) — v(x) = 0,0624x (4 — x), 
z.B. im Punkte *=0,5° — 0,0136< (0,5) — 0,2939<0,0156. 
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Da die Schranken bei x=0 und x=1 den Wert Null annehmen, la8t sich 
auch sofort in diesen Punkten eine Abschatzung fiir die Ableitung w’ angeben: 


— 0,0542 S u’(0) — v’(0) $ 0,0624, —0,0624< w'(1) — v'(1) <0,0542 
oder 
| «’(0) — 1,1570| S$ 0,0583 und |w’(1) + 1,1270] S 0,0583. 


Beisprel 2. Nichtlineare partielle Differentialgleichung in zwei Veranderlichen 
ek. X= V5 
Tu=—Au—w=2—% ftir 0<r<i, 


4.2 
of) (Brehm IB gala ae ets ae 


Gesucht sind Lésungen aus S: 0Sr<1, OSuSX1. Der Ansatz v(x, y)= 
(1—r?) (a+b x*%+ cy?) mit den Konstanten a, 6, c liefert 


Tu—Tv=2—4a+2(b+c) — (1+146 4 2c) x? — (264 14c) y? + v?. 


a, b, c werden aus der Forderung Tu — Tv=0 in den Punkten (0, 0), (1, 0), (0, 1) 
VAN 
1 


= (8 —1/41,5) = 0,5194, b=1—*4 __ 09,0898, c= +b=—0,0065 


12, 12 


= 


bestimmt. Damit wird 


2 
Tu—Tv=—a(i—?r*) E (1 — 7?) (1 a y) | =.0 
und 
Die Tie a(t 7)72 4 — 


hur p—1—/* ergibt sich nach (1.2) L*p=4—20(¢-—77) = 4— 24a (1—7*)? und 
somit 


Tu—Tv a Mate Plas or) 2+(27))| 
Loe = 4—2a(1—7?)? 


lV 


= 70,0220) tur, OS 71 


Weiter ist /,,—=—2, so daB wir mit M,>0, M,=0 und 6,>—0,0220, 6,=0 
nach (2.6) und (2.7) fiir O0<v<1 die Abschatzung erhalten: 


— 0,0220(1 — 7?) S u(x, y) — v(x, y) SO 


mit v(x, y) = (1 — 7) (0,5194 — 0,0898 x2 — 0,0065 y2). 


Zum Beispiel wird 
| ~(0, 0) — 0,5084|<0,0110 und |u(},3) — 0,3657| <0,0083. 


Da die Schranken fiir « die Randbedingungen erfiillen, lieBen sich auf I’ 
auch sofort Schranken fiir die partiellen Ableitungen von w angeben. 
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Uber die Approximation analytischer Funktionen 


in einem veellen Intervall 


G. MEINARDUS 


Vorgelegt von L. COLLATZ 


Die Tschebyscheffsche Approximationsaufgabe, eine gegebene reelle stetige 
Funktion /(x) in einem reellen Intervall [a, b] durch stetige Funktionen w, (x) 
(v=0, 1, ..., #) so zu approximieren, daB der Ausdruck 


eens) I(x) 2% es 
durch geeignete Wahl der reellen Zahlen «, méglichst klein wird, gewinnt in 
neuerer Zeit standig an Bedeutung. Die Berechnung der «, ist meist recht 
schwierig. Setzt man von den Funktionen f(x) und w,(x) nur die Stetigkeit 
voraus, so ist man auf iterative Verfahren zur Bestimmung der «, angewiesen, 
die einen erheblichen Rechenaufwand mit sich bringen. 

Haufig betrachtet man den Fall der Approximation von f(x) durch Polynome, 
deh. man setzt w, (x)=? fiir v=0,4,..., 2. Dies geschieht z.B. oft bei der 
Aufstellung von Unterprogrammen fiir elektronische Rechenmaschinen. In vielen 
Fallen ist nun f(x) eine analytische Funktion, genauer: Es gibt eine Funktion 
f(z), 2=x-+1y, die in einem Gebiet G, welches das reelle Intervall [a,b] im 
Innern enthalt, holomorph ist, und die auf diesem Intervall mit der gegebenen 
Funktion /(x) iibereinstimmt. S. N. BERNSTEIN [7] hat gezeigt, daB die Giite 
der Approximation einer solchen Funktion /(x) durch Polynome wesentlich von 
der Lage der Singularitaten von f(z) abhangt. Fiihrt man die Uberlegungen von 
BERNSTEIN weiter, was hier geschehen soll, so gelangt man zu einfachen Kon- 
struktionsverfahren, die das Polynom bester Approximation mit hinreichender 
Genauigkeit lefern. 

In der vorliegenden Arbeit sollen derartige Verfahren entwickelt werden. Wir 
geben am Anfang einen Uberblick iiber die Grundbegriffe der Theorie in einer 
fiir unsere Zwecke geeigneten Form. Sodann werden Abschatzungsformeln ge- 
wonnen und Konstruktionsverfahren zur Gewinnung eines mdéglichst guten 
Approximationspolynoms aufgestellt, wobei bestimmte Klassen analytischer 
Funktionen im Vordergrund stehen. Zum AbschluB werden noch einige Be- 
merkungen zur numerischen Durchfiihrung der Berechnungsmethoden gemacht. 
Der numerische Aufwand setzt sich zusammen aus der Nullstellenbestimmung 
eines Polynoms mit lauter reellen und einfachen Nullstellen sowie aus der Auf- 
lésung eines linearen Gleichungssystems. Sind die Singularitaten von /(z) weit 
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genug vom Approximationsintervall entfernt, so lassen sich die gesuchten Koef- 
fizienten besonders einfach bestimmen. 

Die hier geschilderten Verfahren verursachen sehr viel weniger Rechenarbeit 
als die bisher bekannten Methoden von MURNAGHAN und WRENCH [6] oder etwa 
die Simplex-Methode des Linear Programming. Man kann sie auch auffassen 
als besonders giinstige Ausgangsnaherungen fiir die oben erwahnten Iterations- 
verfahren. 

§ 1. Grundbegriffe 

Wir betrachten ein Intervall [a,b] der x-Achse und auf ihm stetige Funk- 
tionen f(x). Der Einfachheit halber wahlen wir stets das Intervall [—1, +1], 
auf welches sich jedes Intervall [a,b] mit a6 zuriickfiihren 14Bt. Die Norm 
von f sei erklart durch 


[fl =, Max | f(a). 


#€(—1, +1) 
Wir interessieren uns fiir die Approximation einer Funktion /(x) durch Polynome 
von héchstens m-tem Grad. Ist P,(x) ein Polynom mit Grad P,(x)<n, und gilt 


If—F.lS|f— Qa 


fiir jedes Polynom Q,,(x) mit Grad Q,,(«) <n, so nennen wir P, (x) ein Minimal- 
polynom fiir die Funktion /(x). Bekanntlich gibt es genau ein solches Polynom. 
Wir setzen dann 
Offenbar gilt 

E,,(f +8) SE, (f) + £,(g). (1) 


Es sei L(f) ein lineares stetiges Funktional auf dem Raum C [—1, +1] 
der in —15%<+1 stetigen Funktionen und 


E WSAe tl (80 (tessiOeActentehs (2) 
Dann ist sicher 
|Z(/))| SZ, (A). (3) 
Man kann sich Funktionale mit den Eigenschaften (2) sofort verschaffen, indem 
man von einer Menge {x,} von Punkten x, (w=0,1,...,2+4) mit 
Ss MoS ty SES Hees Sa (4) 
ausgeht und dann fiir w=0,1,...,7-+1 Zahlen i,, aus dem linearen Gleichungs- 
system et 
Ae eae (v= 0,4, 255%) 
u=0 
n+1 (5) 
2(—1)K A, = 4 
u=0 


bestimmt. Hier ist stets A, 4-0 und es gilt 
sgndA,ii=—sgnd, fir w=0,1,...,n. 


Das gewiinschte Funktional lautet dann 
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Nach TSCHEBYSCHEFF gibt es zu jeder Funktion /(x) eine sogenannte Alter- 
nante, d.h. eine Menge {x,,} mit der Anordnung (4), so daB, mit P, (x) als Minimal- 
polynom, die folgenden Beziehungen gelten: 


|#(%,) — F,(%,)| =£,,(f) fiir alle w (7) 
und ; 
f (X41) il 1a (%,41) = —(f (x, = (a) fiir L Sa 0, 1, vee MN. (8) 
Bildet man mit einer Alternante das Funktional L(f) nach der Vorschrift (5), so ist 
IL(\| = E,(/)- 


Man nennt dann L auch ein maximales Funktional beziiglich f. 


Wir ordnen einem Funktional zu einer Menge {x,,} ein Polynom 


Qn(%) = Ditty x” 


zu, indem wir die Koeffizienten «, dem linearen Gleichungssystem 


Dey (SR Te) SO... WE 4) (9) 


entnehmen. Ersichtlich ist 7=L/(f). 
Ist die Menge {x,} eine Alternante fiir f(x), so ist Q,,(¥) das Minimalpolynom 
fiir f(x). 
§ 2. Holomorphe Funktionen 


Von nun an machen wir die folgende Voraussetzung: Die Funktion f(z), 
z=x-+7y, sei holomorph in einem Gebiet der z-Ebene, welches das Intervall 
[—1, +1] im Innern enthalte. Auf dem Intervall [—1, +1] stimme f(z) mit 
der zu approximierenden Funktion f(x) tiberein. 

Die im §1 betrachteten linearen Funktionale zu einer Menge {x,} kann man 
auf dem Unterraum der obigen Funktionen in bequemer Form durch komplexe 
Integrale darstellen. Wir betrachten ein einfach zusammenhangendes Gebiet G 
der z-Ebene, welches das Intervall [—1, +1] im Innern enthalt, und in welchem 
f(z) holomorph ist. Es sei C eine rektifizierbare geschlossene Jordan-Kurve in G, 
welche das Intervall [—1, +1] umfasse. Sodann bilden wir mit einer noch zu 
bestimmenden Konstanten c das lineare Funktional 


Ves f(z) d 
ier tee Gar ee a 
Cc 


wobei der Integrationsweg C im positiven Sinne durchlaufen werde. Durch An- 
wendung des Residuensatzes erkennt man, daB 


n+1 
=>) Yul %,) 
u=0 
ist, mit 
n+1 
Ve edd (4 — 4) = c7,- (11) 
v=0 


ve 
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Da das Nennerpolynom in (10) den Grad +2 besitzt, ergibt sich durch Ver- 
schiebung des Integrationsweges unmittelbar 

Lee) = 0. fie pO, A, ane 2y 
Setzt man nun 


rl n+1 
=D (= 1) Ye» (12) 
so folgt 4 
L*(f)=L(f), 


wobei L(f) das nach (5) zur Menge {x,,} konstruierte lineare Funktional ist. 


Fiir spatere Zwecke betrachten wir zwei Spezialfalle: 


I. Bezeichnet P, (x) das Minimalpolynom fiir die Funktion /(«)=x"~1, so ist 
bekanntlich 
OES eo Rta (a) * 


Hier ist T,,,,(x) das (n+ 1)-te Tschebyscheffsche Polynom erster Art. Wegen 


Ty +1 (Cos g) = cos (n+-1) p 


besteht die Alternante aus den Punkten 


tse Ss 006 mae (4 =0,4,...,% +4). (13) 
Das zugehorige Funktional hat die Gestalt 
Nay f(z) dz 
Ly (/) 21 f (22—1) U,, (2) ” (14) 
Cc 


wobei U,,(z) das n-te Tschebyscheffsche Polynom zweiter Art ist, welches fiir 
z=cos g durch 


__ sin(a+1) 
U,, (cos p) = sing = (15) 
oder allgemeiner durch 
Ree wert Koc a 
U,(@) = — 4 {(e +81) — (Yea, (16) 
oder endlich rekursiv durch 
Upl2)=1, Ujl2)=22, U,r2(2) =22U, (2) —U,@) (17) 


definiert werden kann. Die Formel (14) ergibt sich aus (11) und (12) wegen 


= 


aire - 
nt a i al 
a2 (—1)4 2" oe 
waPRs ae ee | Hie ewe, 
=i etl al > 
oe — fir w=n+1, 


und somit 
C= Dis 


a 
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II. Es sei f())=——, wobei a eine reelle Zahl gréBer als 1 sei. 


Satz 1. Die Alternante der Funktion (x)= —— besteht aus den Punkten 


+41 und —1 und, fiir n>0, aus den Nullstellen des Polynoms 
n Va® — 1 T,, (x) + (a x —1) T;, (x). (18) 


Beweis. Der Fall n=0 ist trivial. Es sei nun »>0. Ist P,(x) das Minimal- 


polynom fiir die Funktion — 


~» so gilt nach S.N. BernstzIn ([1], S. 120/121) 


ee ange Ay" 
x*—a Es, (x) az—1 cos (n Y aig 0), (19) 
wobel 
%=cosm und Wises BUSS 
xA—a 


gesetzt wurde. Die Ableitung der Funktion 


(ax—1) T, (x) +—-Va®—=1 (x2 —1) Ti (x) 
h(x) =cos(ng + 6) = es 


X—a 


hat die Form 


ete REN 
h'(x) = (n 4 EEL Va? 1) Me 


(xa)? 


Sie verschwindet genau z-mal im Intervall (—1, +1) und einmal fir 
x= a+— Va? — 1. Daher gehéren auch die Punkte +1 und —1 zur Alternante. 
Aus (19) folgt (S. N. BERNSTEIN [/]) 


E,( 1 )= (a—Va?—1)" (20) 


%—a a?—1 


Wir bestimmen nun das zu der Alternante des Satzes 1 gehérende Funktional L®. 
Der Ansatz 


Cy eel eres _ f(z) dz : ae yn 
Eyl) ani (z2—1) (x Va®—1 T, (2) +(az—1) Ty (2) oe 


gestattet die Bestimmung der Zahl c®) aus der Forderung 


Le(—*_) = B,(=+,). (22) 


x—a xX—a 


DaB wir hier das negative Vorzeichen zu wahlen haben, wird sich nachher als 
richtig herausstellen. Ersetzen wir den Integrationsweg durch einen Kreis vom 
Radius R mit 1<R<a um den Nullpunkt, und lassen wir R unbeschrankt 
wachsen, so folgt wegen 


T,, (a) = 4 {(a + Va? = 1)" + (@ — Va? = 1)"} 
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nach dem Residuensatz 


(2) i! S —_ cae a : 
En ima (a?—1) (n Va?—1 T,, (a) + (a2—1) Ty (a) 
(2) (a—Va®—1)" ; 


n(a2—1)8 


Mit (20) und (22) wird somit 
c®) =n Va®—1, 


und das Funktional (21) hat die Gestalt 


on Vat A ee f(z) dz a 


ant J (2—1) (nVa®—1 T(z) +(az—1) Ti (2)) * 
Cc 


(23) 


Fithren wir nun in (23) den Grenziibergang a—>co aus, so ergibt sich wegen 
nT, (2) +2; (2) =n U, (2) 


das Funktional L®. Hierdurch stellt sich die Wahl des negativen Vorzeichens 
in (22) als richtig heraus. 

Zu Abschatzungszwecken benétigen wir noch eine weitere Betrachtung. 
Gegeben seien die Funktionen /(z) und g(z). Es sei P,(*) das Minimalpolynom 
fiir g(x) und L das Funktional zu einer Alternante von g(x). Dann ist 


| L(g)| = En (8). 
Nach (3) ist ferner 
ILOISEn)- 


Wir interessieren uns fiir eine Abschatzung der Differenz 


nach oben. Zu diesem Zweck betrachten wir nach (9) das zu der Alternante von 
g(x) gehdrende Polynom Q,,(x) fiir die Funktion f(x). Dann verschwindet der 
Ausdruck 


L(g) (#(*) — Qn(%)) — LY) (g(x (x)) (24) 
fiir die 7+ 2 Punkte x, der Alternante von g( 
L(g) Qn (x) — L(A) B, () 
das diesen Punkten entsprechende Interpolationspolynom der Funktion 
L(g) f(x) — L(f) g(a). 
Fir den Interpolationsfehler ie gilt (vgl. N. E. Nortunp [2], S. 199) 


= Es ist demnach 


Cass 2)—L o(z 
R(#) = nial) gap | SR Ee ae (25) 
mit 
Wn+2(%) = (% — %) (%¥ — 4%)... (% — X41). (26) 


Hierbei soll x innerhalb des von C umschlossenen Gebietes liegen. Setzen wir 
noch voraus, da8 g(x) kein Polynom n-ten Grades ist, also E,,(g)>>0, so ergibt sich 
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aus (24), (25) und (26) die Abschatzung 
Rn 
E,() Slf— Qs || +. (27) 
§ 3. Asymptotische Approximationen 
Wir betrachten eine Folge von Polynomen Q,,(x) mit Grad Q,,(x)<n fiir 
w— (0, 4,2... Ferner sei 


en (f )= a Mt = Qn |. 


Wir nennen 9g, (f) eine asymptotische Approximation von /, wenn 


on(f) =£,(f) (4+o0(1)) fir n> 00 (28) 


ist. Unser Ziel ist es, derartige asymptotische Approximationen zu finden, wobei 
das Hauptinteresse den Polynomen Q,,(x) gilt. 

Hierzu gehen wir folgendermaBen vor Zu der gegebenen Funktion f(x) suchen 
wir eine geeignete Funktion g(x), die durch die analytischen Eigenschaften von 
f(z) 1m wesentlichen charakterisiert wird. Diese Hilfsfunktion muB so beschaffen 
sein, daB zu ihr fiir jedes m Alternanten bekannt sind. Zu diesen Alternanten 
bilden wir nach (5) die linearen Funktionale L, und nach (9) fiir die Funktion 
f(x) die Polynome Q,,(x). Dann bleibt zu zeigen, daB 


en(f) =|L,(f)| 4+0(1)) fir 200 (29) 


gilt, denn dies ist wegen 
[LO = 2) = On) 
gleichbedeutend mit (28). 

Beim Aufsuchen einer geeigneten Funktion g(x) lassen wir uns von den 
Bernsteinschen Gedankengangen leiten [/]. Es gelten die folgenden Satze 
(S. N. BERNSTEIN [J], S. 112/113): 

I. Ist die Funktion f(z) holomorph im Innern einer Ellipse B, mit den Brenn- 
punkten +1 und —1, wobei gq gleich der Summe der Halbachsen ist, und ist f(z) 
stetig auf B,, so gilt 


Peal spe (30) 
Hier ist 
eee ik.) 
II. Ist 
lim sup VE, () = ge1, (31) 


so ist die Funktion f(z) holomorph im Innern einer Ellipse B,, die wie unter I. 
definiert ist. 

Zu einer Funktion f(z) gibt es eine solche Ellipse B, mit maximalem q>1. 
Zuniichst sei g endlich. Wir betrachten den Fall, daB es Funktionen g(z) gibt, 
zu denen wir die Alternanten einfach bestimmen k6énnen, und fiir welche 

f(z) — g (2) 
in einer Ellipse B,, mit ¢,>q holomorph ist. Dann gilt nach (1) die Beziehung 


|Z, — E£, (@)| Ss £,.(f— 8). (32) 
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Aus (32) folgt aber 


E,(f) =£,(8) (1+0(1)) fir 10. (33) 
Denn es ist s 
PRUE «eps 
Elf —8)< Gay 
fiir jedes g, mit g<q,<q und 
M=Ma VANE 
Max |7(2) — 8 @)| 


Andererseits gibt es aber keine Zahl g, mit g,>9, fiir welche 
E,, (8) =9 (gs ") 
gilt, denn sonst ware 


lim sup VE, (g) _S-< 


Docent 
I3 ee 
im Widerspruch zum klassischen Ergebnis (31) von BERNSTEIN. 


Es gilt noch eine etwas scharfere Aussage als (33). Bilden wir némlich zu 
einer Alternante von g(x) ein lineares Funktional L,, nach (5), so ist 


E, (f) 2|L,()| =|Z4,@) + £2, — | 
= rae (f — g)| 
YS Fs Peng I ter 4 


|Z, (| =E,() (14+ 0(1)) fir oo (34) 


Wworaus 


folgt. 


Wir werden also die Konstruktion von Polynomen Q,, (x), deren asymptotische 
Eigenschaft im Sinne von (28) nachgewiesen werden mu8, so vornehmen, daB 


wir eine Alternante von g(x) bestimmen und das lineare Gleichungssystem (9) 
losen. 


Nunmehr betrachten wir einige Spezialfalle. 


1. Die Funktion /(z) sei holomorph im Innern der Ellipse B, mit 
q=a+|a—1; areell, a> 1. 


Abgesehen vom Punkt z=a sei f(z) auf B, holomorph. Bei z=a liege ein einfacher 
Pol vom Residuum A. Wir setzen 


und 


Die Alternante von g(x) wird durch den Satz 1 gegeben. Hierzu gehort das 
Funktional L® nach (23). Es ist nach (20) und (22) 


und 
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Wir haben nun noch zu zeigen, daB 


lf — Q,]= 2,9 (4+ 0(1)) (36) 


ist fiir moo, Dieser Nachweis geht aus von der Formel (27). Da es bei der 
Funktion yw, ..(z) auf konstante Faktoren nicht ankommt, setzen wir 
Z gl ’ 
vl-a(2) = (@— 1) (T,(2) + “7=4_ F(x). (37) 
n Va?—1 


Aus (35) folgt 


1 f L(g) f2)—LOUA) ele) 5, 
sag f EEO a2 =o(E, (0) 38) 
fiir 7—>oo. Man erkennt dies durch Verschiebung des Integrationsweges iiber den 
Pol bei z=a hinaus. Wir verzichten auf die explizite Durchfiihrung dieser Ab- 
schatzung. Ferner ist fiir —1<*x<+1 offenbar 


Yn+2(x) =O(1). (39) 


Aus den Aussagen (35), (38) und (39) ergibt sich mit Hilfe der Formel (27) die 
Abschatzung (36). 


2. Bei der Betrachtung von Funktionen f(z) mit anderen Singularitaten auf 
B, als nur einem einfachen Pol ignorieren wir die hierbei unwesentliche Funktion 
h(z)=f(z) —g(z). Dies stellt keine Beschrankung der Allgemeinheit dar, erleich- 
tert aber einige Abschatzungen. Es sei jetzt 


A 
f(z) = Goa)?’ 
Fiir diese Funktionen kénnen wir explizit keine Alternante angeben. Wir zeigen 
aber, daB wiederum die Alternante des Satzes1 zu asymptotischen Approxi- 
mationen fiihrt. Zunachst ist L')(f) zu berechnen. Eine Verschiebung des Inte- 
grationsweges liefert auf Grund des Residuensatzes fiir das Integral 


DMG 2 Va?—1 A (z—a)*dz 


2010 . z2—1) (nJa®—1 T, (2) + (az—1) Ty (e)) 
den Wert: 
ING WWE pee 
DO) a ee (40) 
Die Abschatzung (27) mit tos 
Oe Seay: Cyd a 
g(j=_, LP) ues 
und pr .2(z) nach (37) liefert nun 
oe Va?—1)?” [  4\2n 
R, () = vhial4) Gog ae =O (la — Va 1)"") 
(a?—1)? (7—a) (41) 


= 0(E,,(g) |LY (f)]) 
fiir n—>oo, gleichmaBig in x fir —1<*<+1. Aus (27), (40) und (41) erkennt 
man 

lf — Q, |= £,, (A (1+ 0(1)) fir n>o, (42) 


wobei Q,,(x) das nach (9) mit der Alternante des Satzes 1 gebildete Polynom ist. 
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Ahbnlich lit sich die asymptotische Eigenschaft von Q,, (x) beweisen, falls 
die einzige Singularitat von f(z) auf B, ein Pol m-ter Ordnung ist. Man kann 
hierzu auch, wie es im folgenden Beispiel geschieht, einen Kunstgriff mit Hilfe 


einer Integration anwenden. 


3. Es sel 
{ (2) =log(a — 2), 


wobei wir den Zweig des Logarithmus wahlen, der fir reelles a mit z<a reell 
ist. Wir legen nochmals die Alternante von Satz 1 zugrunde. Die Durchfihrung 
der Abschatzung (27) soll hier nur skizziert werden. 


Mit festem b> a schreiben wir 
log (a — z) = | com + log (b — 2). 


a 
& 


i 


b 


Dann wird 
b 
Le) (f) = 1 ( nies (140(1)) fir noo. 


y—x 


Somit ist nach (23) und einer kurzen Zwischenrechnung 


b 


iy ee dy 1 | . 
eal reser oy ere) (43) 


a 


fiir 7>co. Man erkennt, etwa durch partielle Integration, 


(a— Va®—1)” j 
= — Onl 
nVa2—1 ( ( ) 
fiir n> oo (vgl. [1], S. 123). Mit Hilfe des gleichen Kunstgriffes laBt sich R,, (x) 
asymptotisch bestimmen. Auch hier folgt, daB 


On (/) = I ~~ 0, | 


eine asymptotische Approximation ist. 


Den Kunstgriff der Integration kann man nicht immer anwenden. Haben wir 
beispielsweise eine Funktion 


f (2) = (a — 2)* log? (a — 2) 


zu untersuchen, wobei « und £ beliebig reell sind, so kann man folgendermaBen 
vorgehen. Man verforme den Integrationsweg, bis man auf einen Schleifenweg 
kommt, der unterhalb der reellen Achse von + co kommend um den Punkt z=a 
im negativen Sinne herumfiihrt und dann wieder oberhalb der reellen Achse bis 
+ co verlauft. Hier denken wir uns den Verzweigungsschnitt von z=a bis z=-+ co 
vorgenommen. Die Verformung ist stets gerechtfertigt, sofern «<n-+41 ist. So- 
dann werden die Integrale fiir L\)(f) und R,, (x) asymptotisch berechnet. 

Es sei hier bemerkt, daB diese Uberlegungen alle zur Begriindung der im 
nachsten Abschnitt zusammengefaBten Konstruktionsverfahren dienen. Sie 
fiihren methodisch einen Schritt iiber die Bernsteinschen Untersuchungen hinaus, 
denn wahrend BERNSTEIN asymptotische Bestimmungen von E, (f) behandelte, 
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ergeben sich hier auch noch Aussagen iiber die asymptotischen Approximationen 
Qn (f)- 


4. Geben wir zum AbschluB dieses Paragraphen noch einen Hinweis auf die 
Behandlung ganzer Funktionen! Betrachten wir fiir eine ganze Funktion /(z) 
die Entwicklung nach Tschebyscheff-Polynomen 


(@)=Dahe, 
so hat BERNSTEIN [1] gezeigt, daB 
on =f) —La.5 


eine asymptotische Approximation ist, daB 


Any ty a: Ef) (1 si 0(1)) 
und 
n, 
} (x) —) a, T(x) > an, +1 T,,, (%) (1 = 0(1)) 
v=0 
gleichmaBig in x fiir —1<*<-+1 fiir eine Folge 1, ,,...,m,>0o gilt. Die 


Extrema der Tschebyscheffschen Polynome bilden daher fiir diese 1, eine asym- 
ptotische Alternante fiir ganze Funktionen. 

Dies wird auch dadurch plausibel gemacht, daB die Nullstellen des Polynoms 
(18) fiir a—oo gegen die Nullstellen des Polynoms 


wT, (x) + %T, (x) =nU, (x), 


also gegen die Tschebyscheffsche Alternante streben. Man kann die Darstellung 
des Funktionals L® in Verbindung mit (27) dazu benutzen, Abschaétzungen und 
asymptotische Aussagen bei speziellen ganzen Funktionen zu gewinnen. 


§ 4. Konstruktionsverfahren 


Das lineare Gleichungssystem (9) liefert uns die Koeffizienten eines Approxi- 
mationspolynoms Q,,(x) fiir eine gegebene Funktion f(x). Die Uberlegungen 
des §3 ermoéglichen nun eine Bestimmung der in (9) auftretenden Gréfen x,. 

Wir bezeichnen die Bernsteinsche Ellipse B, mit maximalem g auch als Regu- 
laritatsellipse von f(z). 

Hat /(z) auf der Regularitatsellipse genau eine Singularitat bei z=—a, a reell 
und >1, so wahle man die x, nach Satz 1 als Nullstellen von (18). Auf eine 
bequeme Berechnungsart gehen wir noch ein. Liegt die Singularitat bei z= —a, 
a>1, reell, so nehme man natiirlich die am Nullpunkt gespiegelten ¥,,. 


Ist f(z) eine ganze Funktion, so wahlt man die x, nach (13). 


Liegen mehrere Singularitaten auf der Regularitatsellipse, so laBt sich im 
allgemeinen keine asymptotische Alternante bestimmen. Es gibt aber wichtige 
Spezialfalle, z.B. 


+ h(z) 


oder 


f(z) = aor + hie), 
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bei denen es offenbar méglich ist, eine Zuriickfiihrung auf die Alternante von 
Satz 1 zu erreichen. Hierzu betrachtet man die Approximation von 
A A 
Fae 

durch Polynome vom Grade n/2 im Intervall [0,1], denn die Approximation 
von g(x) in [—1, +1] ist ein Polynom in x*. Durch eine lineare Transformation 
erhalt aus Satz 1 die Alternante fiir das Intervall [0,1], woraus dann endlich 
fiir die urspriinglichen Funktionen /(x) asymptotische Alternanten folgen. 

In allen anderen Fallen geht man zweckmaBig von der Tschebyscheff-Alter- 
nante (13) aus und fiihrt noch eine direkte Verbesserung dieser Punkte durch, 
die im folgenden beschrieben werden soll. Das gleiche Verfahren kann auch bei 
ganzen Funktionen Verwendung finden. 

Es sei fiir den Augenblick 


2 
x,—= — COS 


Fi PEON, i= Os aes Ha Al: 


Dann erklart man sog. Grundlésungen, das sind +1 Polynome n-ten Grades 
v,(x), durch die Vorschrift 


U, (%y) =f Uy (%,4+1) 7 Dn (44) 


fiir vy, x=0,1,..., 2. Diese Polynome sind explizit dargestellt in [4]. Macht 
man nun fiir Q,,(%) den Ansatz 


nN 


Qn (%) = 2, Ba Yn (%), (45) 


Oo 


so lautet das lineare Gleichungssystem, wenn wir je zwei aufeinanderfolgende 
Gleichungen addieren, 


=F (oats A eel ge ay PAA (46) 


und es wird 
AA)+2E (—1)" t,)+(— 1)" (9) 
z Bn? 


Man kann nun, wie in [4] naher beschreiben, eine verbesserte Alternante x, durch 
die Vorschrift 


~ Pra Qh (yeah im) a 

Lu le Opla,) —#" (x4) =) [ul = tye Pe. elie) lg nN (47) 
gewinnen, wobei %)= %)=—1, und %,.,=%,4,—-+41 unverdndert bleiben. Zur 
Berechnung der Zahlen %,, brauchen wegen (45) nur die Werte v,(%,) und Vn (%) 
vorzuliegen, abgesehen von den Ableitungen f'(x,) und f’’(%,). 


. Berechnet man nun wiederum Q,,(x) nicht direkt nach (9), sondern aus dem 
linearen Gleichungssystem, welches sich aus dem Ansatz (45) ergibt, so wird 
dieses neue Gleichungssystem nach Addition je zweier aufeinanderfolgender 


Gleichungen im allgemeinen zur Iteration nach dem GauB-Seidel-Verfahren ge- 
eignet sein. 
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Die letzte Bemerkung gilt auch fiir den Fall, daB die GréBe g der Regularitats- 
ellipse nicht zu nahe bei 1 liegt, also etwa bei der Alternante des Satzes 1, falls 
a nicht zu nahe bei 1 liegt. Man vermeidet durch den Ansatz (45), daB das 
Gleichungssystem zur Berechnung der Koeffizienten von Q,,(x) illconditioned 
wird, was bei der Gestalt (9) fiir nicht zu kleine Werte von 1 zwangslaufig der 
Fall ware. : 


Allgemein kann gesagt werden, daB das auf diesem Wege berechnete (Q,, (x) 
in vielen Fallen bereits eine hinreichend gute Approximation fiir f(x) liefert. 
Man kann es stets als eine besonders giinstige Anfangslésung ansehen, wenn 
man etwa ein Iterationsverfahren nach der Austauschmethode von E. REMEz [5] 
bzw. von MURNAGHAN und WRENCH [6] im Auge hat. 

Die Bestimmung der Nullstellen des Polynoms (18) kann fiir groBe Werte 
von ” Schwierigkeiten bereiten. Durch den folgenden Kunstgriff lassen sich 
diese Zahlen jedoch miihelos berechnen. 

Es sei x,, die y-te Nullstelle des Polynoms (18) in der Reihenfolge —1< x, < 
Whe ea rie ot. nd 


Wir wissen, daB x, gegen 0, strebt fiir aco. Man setze a=1/v und fasse die 
Stelle x, als Funktion von v auf: x,—%,,(v). Diese Funktion ist implizit gegeben 
durch die Identitat 


n(1— v2)! T, (x, (v)) + (x, (v) — v) Ty (x, (v)) =0. (48) 


Durch Differentiation nach v ergibt sich 


ax i 
Fk = mo(t— oT, (x,) + 


nm (1 — v?)8 T, (a) a 


Beriicksichtigt man noch 
(1— 24) Ti (x) = x Ti (x) — wT, (2), 


so ergibt sich die zweite lineare Beziehung zwischen T,, und 7, : 


ay {n v(1— v2)—2 + LACS pS AED) Os ape aot bb Cae 
ax (49) 
| { 1+ (n(1— v2) +14 x, (x, —v) (1 — x7)) set iH) eo 


Die Bedingung fiir die nichttriviale Lésbarkeit des aus (48) und (49) gebildeten 
linearen homogenen Gleichungssystems fiir 7, («,) und T;, (x,) liefert die folgende 
Differentialgleichung fiir x,, (v): 


AK y 


| : Laren 50 
dv nV1—v?(1—v %,) +1—v? eo) 
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Wir betrachten nun das Anfangswertproblem dieser Abelschen Differential- 
gleichung zweiter Art mit der Anfangsbedingung x,(0)=g,- Hier kann man 
irgendein Naherungsverfahren verwenden, um x, an der gewiinschten Stelle v 
zu berechnen. Bequemer, wenn auch nicht von hoher Genauigkeit (die meist 
nicht erforderlich ist), ist ein Potenzreihenansatz 


xy (v) = Ou iz C10 


dessen erste Koeffizienten rasch gefunden werden kénnen. Auf diese Weise 
ergibt sich, wenn wir jetzt wieder a an Stelle von 1/v schreiben, fiir die y-te 
Nullstelle des Polynoms (18) die folgende asymptotische Entwicklung 
eS SUE Oa ene een 
u~Cut Gari! a(n+1j22! | 
, Aen) {2 of (n—2) +n? + 5n+2} 
a? (n+ 1)3 3! 


(54) 


peer 


Auf den naheliegenden Beweis, daB es sich bei (51) tatsachlich um eine asym- 
ptotische Entwicklung handelt (sie ist sogar fiir |a|>1 konvergent), soll hier 
verzichtet werden. Liegt a nahe bel 1, a=1-+¢, so lassen sich ebenfalls all- 
gemeine Naherungsformeln fiir die «, gewinnen. 


Tabelle 


der Nullstellen des Polynoms (18) fiir einige Werte von a 


a=1,5 a—2,0 C= Des = one) Cs a=4,0 

m=1 | +0,38197 | +0,26795 | +0,20871 | +0,17157 | +0,14590 | +0,12702 
n=2 | —0,30902 | —0,36603 | —0,39564 | —0,41421 | —0,42705 | —0,43649 
+0,69098 | +0,63397 | +0,60436 | +0,58579 | +0,57295 | +0,56351 

m=3 | —0,61803 | —0,64329 | —0,65685 | —0,66551 | —0,67157 | —0,67607 
+0,19098 | +0,13397 | +0,10436 | +0,08579 | +0,07295 | +0,06351 

__| +0,80902 | +0,77726 | +0,76121 | +0,75130 | +0,74452 | +0,73957 
n=4 | —0,76226 | —0,77503 | —0,78204 | —0,78658 | —0,78978 | —0,79217 
—0,17784 | —0,21545 | —0,23553 | —0,24830 | —0,25721 | —0,26380 
10,45325 | +0,40900 | +0,38640 | +0,37236 | +0,36273 | +0,35568 
+0,86882 | +0,84943 | +0,83989 | +0,83409 | +0,83016 | +0,82731 

m=5 | —0,83883 | —0,84606 | —0,85010 | —0,85273 | —0,85460 | —0,85600 
—0,41355 | —0,43748 | —0,45054 | —0,45894 | —0,46485 | —0,46925 
+0,12595 | +0,08884 | +0,06935 | +0,05707 | +0,04856 | +0,04229 

+ 0,604 53 | +0,57145 | +0,55489 | +0,54473 | +0,53780 | +0,53276 
+0,90387 | +0,89120 | +0,88511 | +0,88145 | +0,87899 | +0,87722 

m=6 | —0,88386 | —0,88833 | —0,89085 | —0,89250 | —0,89367 | —0,894.56 
—0,56526 | —0,58096 | —0,58966 | —0,59530 | —0,59930 | —0,60229 

— 0,124 91“! —0,15262 | =0)16753- | O49 05 —0,18370 | —0,18863 
+0,32990 | +0,29730 | +0,28049 | +0,27001 | +0,26280 | +0,25752 
+0,69976 | +0,67493 | +0,66273 | +0,65533 | +0,65032 | +0,64670 
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Einen groben Uberblick iiber die Nullstellen x, von (18) kann man sich auch 
graphisch verschaffen, indem man 


cies (52) 


schreibt. Da die linke Seite von (52) von a unabhangig ist, braucht man nur die 
Schnittpunkte ihres graphischen Bildes mit einer Geraden zu bestimmen. Die 
Konstruktion wird noch einfacher, wenn man beachtet, daB die Enveloppe der 
Geradenschar 


Va —1y+-ax=1 
die gleichseitige Hyperbel x2— y?=1 ist. Die Zusam- 
menhange sind leicht aus der Abb. 1 zu erkennen (dort 
ist. —4). 

Besonders wichtig, weil zeitsparend, sind die asymp- 
totischen Aussagen iiber die GréBen E,,(f), die man 
haufig fast ohne Rechnung gewinnen kann. Sie geben 
einen guten Anhaltspunkt, wenn man eine Abschatzung 
fiir den Grad des anzusetzenden Polynoms haben 
moéchte, um eine vorgegebene Genauigkeitsschranke — Abb. 1. Graphische Bestimmung 
nicht zu wtiber- und nicht zu stark zu_ unter- Se ASSES Pay) 
schreiten. 

Es ware wiinschenswert, ahnliche Aussagen tiber asymptotische Alternanten 
auch bei rationalen Approximationen analytischer Funktionen zu gewinnen. 


SS 
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§ 5. Beispiele 

1. Die Funktion J"(x) soll im Intervall [1, 2] durch ein Polynom vom Grade 4 
approximiert werden. Die ausschlaggebende Singularitat ist der Pol erster Ord- 
nung bei x=0. Transformieren wir auf das Intervall [—1, +1], so liegt der 
entsprechende Pol bei — 3. Hieraus entnimmt man die asymptotische Alternante 
nach (48), die man wieder auf das obige 
Intervall zu transformieren hat. Die 
Lésung von (9) ergibt das Polynom 


QO, (x) = +3,099 63982 — 
— 4,467 42700 x —3,375 48216 x? 
— 1,182 467 77 x8 -++0,17457650 x4 


und Abb. 2. Q,(x) —L' (x) 


|7| = 1,9629 - 10 


als untere Schranke fiir die Approximation. Man rechnet nach 


Max |I"(x) — Q,(x)| < 1,9680- 1074, 


x €(1, 2] 


so daB der prozentuale Fehler dieser Approximation gegeniiber der bestmég- 
lichen kleiner als 0,26% ist (s. Abb. 2). 
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2. Die Funktion 1 


=| 
soll im Intervall [—1, +4] durch ein Polynom vom 4. Grade approximiert 
werden. Die Funktion besitzt bei +227 und bei — 2z7 je einen einfachen Pol. 
Zur Summe der Hauptteile laBt sich die Alternante bestimmen. Damit ergibt 
sich das Polynom 


OF (x) = +1,000000983 2—0,5000000000 x + 0,083 315.479 7 x®—0,001 3407391 x4 
und 
[7p] 30,8314 A0- 
Es ist nun 
Me ()¥ 
aay Qa (x) 
so daB der prozentuale Fehler gegeniiber 
Abb. 3. e)-=— der besten Approximation kleiner als 
: 0,14% ist (s. Abb. 3). 

3. Wir approximieren die Funktion sinx in dem asymmetrischen Intervall 
[0, 2/2] (was normalerweise nicht sinnvoll ist) durch ein Polynom 3. Grades und 
benutzen als Alternante, da es sich um eine ganze Funktion handelt, die trans- 
formierte Tschebyscheff-Alternante. Es ergibt sich 

Q3 (x) = — 0,001 358 + 1,025 225% — 0,070685 x? — 0,112497 x8 


< 9,845 - 1077, 


und 
|| = 41,358°40°. 
Dardas 


Max |sin x — Q,(x 277 .40-3 
dia, [sins Qs (x)| < 1,377 -10 


ist, ergibt sich ein prozentualer Fehler gegeniiber der besten Approximation von 
héchstens 1,4% (s. Abb. 4). 


-F 


1 m 
0 10% 
if 7 
! ! 
0 eS 7 425 45 1,68 
-~7 ef: 
Abb. 4: Q,(x)—sin x Abb. 5. Q.(x)—log x 


4. Bei der Approximation von log x im Intervall [1, 1,68] legt die kritische 
Singularitat bei »=0. Die Transformation auf das Intervall [45-444 ) and 


die Bestimmung einer asymptotischen Alternante auf Grund von (18) ergibt 
das Polynom (fiir 7 = 2) 


Q, (x) = — 1,244809 + 1,536469 x — 0,290 217 x2 
und 


[7 | = 1,443 - 410-3. 
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Da das 
- : = 
Pies llog i Qs () Z tie 1 


ist, ergibt sich eine prozentuale Abweichung von der besten Approximation um 
hdéchstens 0,14% (s. Abb. 5). 
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Zur Theorie der Approximationen 
vermittels stetiger Operatoren 
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Von P. L. Burzer [1] und E. H1rre [2] wurde das Randverhalten der Poten- 
tialfunktionen am Beispiel des Einheitskreises untersucht. Aus Annahmen iiber 
den Grad der Approximation der Randwerte ergeben sich dabei Folgerungen fiir 
die Existenz der Ableitung in Normalenrichtung. Diese Ergebnisse werden durch 
eine Kombination von Methoden der Theorie der Halbgruppen mit Begriffs- 
bildungen der Approximationstheorie erhalten. Fassen wir die vorgegebenen 
Randwerte als Elemente eines reflexiven Banach-Raumes L? (1< p< oe) auf, so 
gilt mit den Bezeichnungen der klassischen Potentialtheorie: 


4. Notwendig und hinreichend fiir (7 > 1—0) 


27 


i 


0 


27 


5 (1—72)dq@ . P sy, #3 
[ @) parece | =O} 


1 
Dobe 
) 


ist die Existenz einer Funktion /(y) aus L?, so daB 


2a Qn 
tin, [de & fe ge —y — Men] ayo 
gilt. . : 
2. Notwendig und hinreichend fiir (y + 1— 0) 
27 Qa ; ; 
/ Fe | &(9) Seana = Bey = aia —7) 


0 0 
ist g(p) =const. 

Diese Ergebnisse regen die folgenden Uberlegungen an: 

Die Beschrankung der Aussagen auf L’-Normen mit 1<~< oe scheint nicht 
durch die Natur der potentialtheoretischen Fragestellung, sondern durch die 
Technik des Beweises begriindet zu sein. 

Es ist klar, daB Ergebnisse dieser Art nicht auf den Fall des Kreises beschrankt 
sind. Die Methode der konformen Abbildungen lat naémlich unmittelbar er- 
kennen, daB sich analoge Aussagen fiir eine wesentliche gréBere Klasse von 
Gebieten formulieren lassen. Andererseits ist es nicht schwer, die entsprechenden 
Ergebnisse auch fiir den Fall der Einheitskugel zu formulieren. Obwohl das 
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Mittel der konformen Abbildungen im Raume die Ergebnisse nicht unmittelbar 
auf eine gréBere Klasse von Gebieten tibertragen laBt, kann doch kein Zweifel 
daran bestehen, daB entsprechende Satze auch fiir allgemeinere Gebiete gelten 
miissen. 

In der Terminologie der Operatoren in Banach-Raumen besteht das Anliegen 
der genannten Probleme darin, aus Annahmen iiber den Differenzenquotienten 


R(x, h) @= 5 {T(x +h) g — T(x) g} 
fiir x>0 die Existenz des Differentialquotienten im Sinne von 


lim | {T(h) g— L(O)e} 


h—-0 


zu schlieBen. Der Beweis des ersten Ergebnisses benutzt folgende Voraussetzun- 
gen: 
1. Der zugrunde gelegte Banach-Raum ist reflexiv. 


2. Der Operator T(x) erfiillt die Voraussetzungen: 
a) T(x) ist stetig (s— lim T(x) g=T(x))g); F(0) =I. 


b) T(x+y)=T(x) T(y) (Halbgruppeneigenschaft). 
3. Ai) e= 5 


{T(h) g—g} ist beschrankt fiir kleine h. 


Aus den Voraussetzungen 2 und 3 folgt fiir die Elemente /(x)=T(x)g, daB 
der Differenzenquotient 


_ {f(x +h) —f(x)}=R(x,h) g= T(x) A(h) g 


gleichmaBig beschrankt und beziiglich / gleichgradig stetig in x ist. Die aus 
der klassischen Theorie der reellen Funktionen bekannte Technik des Auswahl- 
satzes ist somit anwendbar, wenn bei festem x>0 die Menge R(x, h)g kompakt 
ist. Diese Forderung wird durch die erste Voraussetzung erfiillt, da reflexive 
Raéume im schwachen Sinne kompakt sind. 

Dies zeigt, daB die bisher erzielten Ergebnisse unter wesentlich allgemeineren 
Voraussetzungen zu erhalten sind, wenn die gleichmaBige Beschranktheit und 
gleichgradige Stetigkeit des Differenzenquotienten, sowie die Kompaktheit der 
Menge R(x, h)g bei festem «> 0 garantiert sind. 

Eine erste Verallgemeinerung bietet sich dadurch an, daB wir die Forderung 
der Kompaktheit nicht durch Voraussetzungen iiber den Banach-Raum, sondern 
durch Annahmen iiber den Operator erzielen. Dies erscheint deshalb fiir die 
Fragen der Potentialtheorie besonders naturgemaB, weil jede harmonische Funk- 
tion im Inneren ihres Definitionsbereiches analytisch ist, so daB der Nachweis 
der Kompaktheit der Menge R(x, /)g nicht schwierig ist. 

Die vorliegende Arbeit untersucht die Differenzenquotienten R(x, h)g und 
leitet allgemeine GesetzmaBigkeiten daraus ab, daB diese Folge gleichgradig 
stetig, gleichmaBig beschrankt und fiir x>0 kompakt ist. Die Voraussetzung 
der Halbgruppeneigenschaft ist als Spezialfall darin enthalten. Es enthalt §4 
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die Beweise der verallgemeinerten Fassung der Satze aus der Approximations- 
theorie. Im §2 werden diese Ergebnisse dann auf die Untersuchungen des 
Randverhaltens der Potentialfunktionen im Einheitskreis angewandt. Dabei 
konnen die Ergebnisse von ButzerR und HiLte auf den Raum der stetigen Funk- 
tionen und alle Riume L? (1<p<ov) iibertragen werden. Die Ubertragung der 
potentialtheoretischen Ergebnisse auf allgemeine raumliche Gebiete wird in der 
nachfolgenden Arbeit [4] des einen von uns vorgenommen. 


§ 1. Approximationssatze 


Wir wollen nun zeigen, daB sich die anfangs erwahnten Ergebnisse unter 
wesentlich allgemeineren Voraussetzungen erhalten lassen, und gehen dazu 
folgendermafen vor: 

Es sei 8 ein Banach-Raum und 7(x) ein Operator, der von dem reellen 
Parameter x (0<.x< co) abhangt und der 8% in sich abbildet. Wir setzen 


ey! 
(1.1) R(x,h) g= 5 


{T(x + h) g — T(x) g} 
und beweisen 
Satz 1. Es geniige der Operator T(x) den Voraussetzungen 
a) T(x) g ist beziighich x stetig fir alle g. 
b) R(x, h)g ist gleichmaBig beschrankt und beziiglich h gleichgradig stetig in x. 


c) Fir jedes x>0 ist die Folge R(x, h)g fiirh—>0 kompakt (im starken oder 
schwachen Sinne). 


Dann gibt es eine stetige Funktion q(t) mit 


Zum Beweis wahlen wir auf der positiven reellen Achse eine abzahlbare, 
iiberall dichte Menge ® und nehmen zunichst an, R(x, 2) g sei im starken Sinne 
kompakt. Dann kénnen wir eine Folge h, so auswahlen, daB R(x, h,)g in allen 
Punkten von & konvergiert. Es gilt also etwa 


(4:2) R(x, h,) g=9,(x) +(x) fir «ER. 


Ist x9 ein beliebiger Punkt ¢8, so wahlen wir zur Definition von (xj) eine 
Folge x,C8 mit x,—2%). Dann ergibt sich aus Voraussetzung b) und (1.2) 
|p (x,) Sale (x,,)|| = |p (%,) geal ts (x,)|| oi 


(1.3) 
hs lp (%,) — Pn (x,)|| ale Il, (x,) —_— Pn (x,)| = € 


’ 


fiir alle |x, —x)|<6(e) 


%,—% |Sd(e) und alle n=N(e). Wir setzen daher 


(1.4) P(Xo)\= lim @(x,). 


Xy— Xy 
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Die so gebildete Funktion g(x) ist eindeutig bestimmt und stetig. Bilden wir 
namlich 
+ |p (%2) aa p (x,)|| + le) — 9(x,)| Eee eis 


und wahlen |x,—¥,|, |%.—~,| geniigend klein, so folgt aus (1.3) 
(1.6) |p (%1) — p(*2)| Se. 


Nun denken wir uns nach dem Diagonalverfahren eine Nullfolge h,,>0 so 
ausgewahlt, daB R(x, h,,) g=q, (x) in den Punkten 


Bes, (F=0,1,..., acl} 


aus ® konvergiert. Dann gilt fiir #,,>0 auf Grund der Definition des Integrals: 


x 


Jf pl) dt= hy 2s p (En) +O (Fr) 


0 


(1.7) = DT (En + hy) = T(&)} e+ 9(4) 
=T (hy [7] +i] 8 — T(0) ¢ + 0(4). 
Fiir /,, 0 folgt aus (1.7) und der Stetigkeit von T(x) 


(1.8) J p(t) dt= T(x) g—T(0)¢g; 


damit ist Satz1 im Falle starker Kompaktheit bewiesen. Gilt diese Voraus- 
setzung nur schwach, so folgt entsprechend fiir alle 6 aus dem konjugierten 
Banach-Raum %8* 


x 


(1.9) (0, fo at) = (0, T(x) g— T(0) 8). 


0 
Nach dem Satz von HAHuN-BANACH ist dies gleichbedeutend mit Gl. (1.8), 
und Satz 1 ist auch im Falle schwacher Kompaktheit bewiesen. 


Wir beweisen ferner 

Satz 2. Der Operator T(x) geniige den Voraussetzungen: 

a) T(x)g set beziiglich x stetig fiir alle g. 

b) Fiir jedes x>0O existiere eine Folge h,,, so dak R(x,h,)g gegen Null konver- 
giert (im starken oder schwachen Sinne). 
Dann gilt T(x) g=T(0)g. 

Zum Beweis gehen wir wie bei Satz 1 vor. Die Funktion @ (?) verschwindet 
in diesem Falle, woraus die Behauptung folgt. 
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Zur Verallgemeinerung der eingangs genannten Approximationssatze fiihren 
wir nun eine Klasse K von Operatoren T(x) ein, die durch die folgenden Eigen- 


schaften definiert wird. Es sei 


(1.10) 5 {T(x +h) g — T(x) g} = R(x, h) = T(x) A(h) + B(x, h) 
mit 

(1.11) A(h) = ; {T(h) —DB 

und 

(1.12) lim A(h)g=Ag; gE D(A) 


falls dieser Limes existiert. Der Operator T(x) sei stetig; es gelte 7(0)=—J, und 
B(x, h)g sei gleichmaBig beschrankt und beziiglich / gleichgradig stetig in x. 


Der Fall, da8 7(x) Halbgruppeneigenschaft besitzt, ist hierin enthalten. 
Im Falle einer Halbgruppe gilt namlich R(x, h)= T(x) A(h). 

Wir beweisen nun 

Satz 3. Der Operator T(x) gehore zur Klasse K. Der zugrunde gelegte 
Banach-Raum 8 set reflexiv. Dann ist die Bedingung g© D(A) notwendig und 


hunretchend fiir 
[7 (x) ¢—gl]=0() 
Liegt g in D(A), so folgt die Behauptung von Satz 3 unmittelbar. Gilt um- 
gekehrt 
(1.43) 1A(2) SC, 


so erhalten wir 


|R (41,2) g — R(x2,h) g|| <|(T(4) — Tee h) g|+ 


(1.14) 
+ || (B(x, h) — in eine 

fiir alle # und |x,—x,|<6(e). R(x, h)g ist also gleichmaBig beschrankt und 

beziiglich / gleichgradig stetig in x. Da % als reflexiv vorausgesetzt wurde, ist 

R(x, h)g nach (1.10, 13) beschrankt, also schwach kompakt. Die Voraussetzungen 

von Satz 1 sind somit erfiillt, es gibt daher eine stetige Funktion gy (¢) mit 


(1.15) A(x) g=— [ pat 
und ; 
(1.16) lim A(x) g=Ag=9(0). 


Damit ist Satz 3 bewiesen. 


Als nachstes verzichten wir auf die vorausgesetzte Reflexivitat des Banach- 
Raumes und ersetzen sie durch Vollstetigkeit des Operators. Es gilt 
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Satz 4. Der Operator T(x) gehére zur Klasse K und sei fiir jedes x>0 voll- 
stetig. Die Folge B(x, h)g (h->0) set fiir jedes x > 0 kompakt. Dann ist die Bedin- 
gung gC D(A) notwendig und hinreichend fiir 


T(x) g—g|=O(»). 
Gehoért g zu D(A), so folgt Satz 4 wieder unmittelbar. Es sei also 
(1.17) lA slSC. 


Dann finden wir analog zu (1.14), daB R(x, h)g gleichma&Big beschrankt und 
beziiglich / gleichgradig stetig in x ist. Wegen 


(1.18) R(x, h) = T(x) A(h) + B(x, h) 


erhalten wir aus (1.17) und den Voraussetzungen tiber T(x) und B(x, h), daB 

R(x, h) g fiir jedes x >0 kompakt ist. Die Voraussetzungen von Satz 1 sind also 

wiederum erfiillt, Gl. (1.15, 16) gelten daher. Satz 4 ist damit ebenfalls bewiesen. 
SchheBlich zeigen wir 


Satz 5. Der Operator T(x) gehore zur Klasse K. Die Folge B(x, h)g (h->0) 
set fiir jedes x > 0 kompakt. Dann folgt aus 


lim ~ T(x) g—g||=0 
bereits Ag=(. 

Zum Beweis von Satz 5 benutzen wir die Ergebnisse von Satz 1 und Satz 2. 
Das Restglied R(x, h) setzt sich nach (1.10) aus den Anteilen F(x) A(A) und 
B(x, h) zasammen. Der erste Anteil geniigt den Voraussetzungen zu Satz 2, der 
zweite denen von Satz1. Wahlen wir daher wie beim Beweis zu Satz1 die 
Teilfolge h,, zunachst so, daB A(h,)g gegen Null konvergiert, dann erhalten wir 
wieder 


(1.19) T(x) e—e=S (0) at, 


wobei die stetige Funktion @ (¢) nur durch den Anteil_B (x, h) bestimmt ist. Nach 
(1.10) gilt aber fiir alle 


(1.20) B(0, h) =0. 
Aus (1.20) folgt y(0)=0 und daraus 
(1.21) lim A(x) g=Ag=0. 


Damit ist auch Satz 5 bewiesen. 


§ 2. Das Randverhalten der Potentialfunktionen im Einheitskreis 


Mit Hilfe der Poissonschen Formel 
2% 
jedd 
(2.1) Ury)=s- | 89) 
0 


eg A ane 
1+7?—2rcos(p—y) 
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lésen wir die Dirichletsche Randwertaufgabe fiir Funktionen g (y) aus dem Raum C 
der stetigen Funktionen oder aus L? (1 pce). Der Operator T(x) wird durch 


(2.2) COG v) —Le)2 x=—lIgr 
definiert. T(x) ist beziiglich x stetig und geniigt der Beziehung 


T(x+y)=T(x)T(y); T(0)=!, 
(2.3) 


auf deren Nachweis wir im einzelnen hier verzichten. Setzen wir weiterhin die 
Beschranktheit von 


22 
(2.4) {L-sup|U(r, 9) —e(p)| baw. SE [10 #) —8l@)lPag 
BY Gy eu 
0 


fiir y>1—O voraus, so bedeutet das, daB der Operator 
(2.5) A(h) ¢=5 {T(h) — Tg 


im Sinne der C bzw. L’-Norm beschrankt ist. 

Die Voraussetzung der Kompaktheit erhalten wir nach dem Vorgang von 
Butzer, HILLeE und DE LEEuw [3] dadurch, da wir uns auf reflexive bzw. 
konjugierte Banach-Raume, d.h. 1<p<co beschranken. Nach Satz 3 folgt 
dann g€ D(A), es existiert also im Sinne der zugrunde gelegten Norm 


: 1 : 0 
(2.6) lim — (T(x) g—g) =lim U(r, 9). 


Die Beschrankung auf reflexive Banach-Raéume entfallt, wenn wir die Voll- 
stetigkeit der Operatoren T(x) fiir «>0 benutzen. Die Vollstetigkeit ist leicht 
einzusehen, da vermittels des Poissonintegrals 

27 
2. ya [s eberw Ce Pts, 
( 7) fn(7, P) 2n : 8n(@) 1+72—2rcos(p—y) 

0 
jede gleichmaBig beschrankte Folge g,(y) aus C oder L? (1< pov) fiir r<1 
in eine gleichgradig stetige transformiert wird. Aus Satz 4 folgt daher wieder 
gCD(A), der Differentialquotient 


2.8 im © ) 

8) lim 5, UU. 9) 

existiert also im Sinne der zugrunde gelegten Norm. 
Wegen 

(2.9) Igy=O(1—7) fiir 11 


haben wit damit fiir den Raum C der stetigen Funktionen und alle L?-Raume 
gezeigt, daB die Existenz des Grenzwertes (2.8) notwendig und hinreichend fiir 


(2.10) IU, 9) — ¢(|=00—7) 
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ist. Aus diesem Ergebnis ersehen wir speziell, da8 notwendig und hinreichend fiir 


(2.11) IU, &) —g@|=ot—7) 
die Bedingung 

| mca = 
(2:12) ; lim ae p= 0 


ist. In der Potentialtheorie folgt aber aus dem Verschwinden der Normalableitung 
daB U(r, y) konstant ist. Dann ist auch g konstant, und diese Bedingung ist 
notwendig und hinreichend fiir (2.11). 
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Uber das Randverhalten harmonischer Funktionen 


Rour LEIS 


Vorgelegt von C. MULLER 


§ 1. Einfiihrung 
Es sei g(p) eine auf dem Einheitskreis definierte Funktion. Dann erhalt 
man durch das Poisson-Integral 


(1 =P) ap 


——. — - —- VA 
1+72—2y7 cos(py—y) 


(1.4) Ur.v = | 8) 


die Potentialfunktion fiir das Innere des Einheitskreises, die die Randwerte g (¢) 
annimmt (Dirichletsche Randwertaufgabe). Das heift, es gilt 


(1.2) lim U(r, y) = g(y)- 


Wir wollen das Verhalten von U(r, y) am Rande in Abhangigkeit von g(y) 
genauer diskutieren, insbesondere die Frage nach der bestméglichen Approxi- 
mation der Randwerte stellen und notwendige und hinreichende Bedingungen 
dafiir finden. 


Dazu wahlen wir statt 7 den Parameter x= —lgv (0S x%< ce) und schreiben 
fiir (1.1, 2) symbolisch 
(1.3) U(r,y)=T(x) g(y); lim T(x) g=¢. 


Legen wir fiir g den Banach-Raum der in p-ter Potenz integrierbaren Funktionen 
zugrunde (gE L?, 1<p<oo), dann gilt nach BurzeEr [J] 

Satz 1. Notwendig und hinretchend fiir 

Z(*) ¢— gl]=o(x) 

ist g=const. 
und 

Satz 2. Notwendig und hinreichend fiir 
| IZ ¢—sl=0() 
isto Ce 

Fiir eme Approximation der Randwerte von der Ordnung O(x) ist also die 
Differenzierbarkeit von g notwendig und hinreichend. Werden die Randwerte 
von der Ordnung o(x) approximiert, dann ist g bereits konstant. 


Der Beweis dieser Satze erfolgt durch Anwendung von Ergebnissen aus der 
Halbgruppentheorie. Man kann namlich leicht nachrechnen, daB fiir den durch 
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das Poisson-Integral definierten Operator 


bi) I(x + 9) g=T(x) Ty) g 
Soin (d)es= Tape dalle..)=>0; 


(1.5) A(x) g=— (T(x) —1)g 
und g€D(A) falls der Limes 
(1.6) Ag=s—lim A(x) g 


existiert. Dann gilt (vgl. Butzer [1], Hitre & Puixtirs [4]) 


Lemma 1. Es sei S$ ein Banach-Raum. Dann verschwindet 


lim + T(x) ¢ — ¢| 


x0 *¥ 
genau dann, wenn g invariant unter T(x) ist 
und 
Lemma 2. a) Es set g€D(A). Dann folgt 


|T(*) ¢ — g|=O(2). 
b) Es set 8 reflexiv, und es gelte 


I7() ¢—g|=O(%), 
dann folgt g€D(A). 

Wir ersehen also aus Lemma 2, daB die Bedingung g€D(A) fir reflexive 
Raume notwendig und hinreichend fiir ein Approximationsverhalten von der 
Ordnung O(x) ist. Als nachstes taucht daher die Aufgabe auf, den Definitions- 
bereich von A festzulegen. Fiir das Poisson-Integral folgt durch Anwendung 
von Ergebnissen aus der Theorie der Fourier-Reihen [7], daB sich in den Banach- 
Raumen 8=L? (1<p< ov) die Eigenschaften g€D(A) und g’€L? entsprechen. 

Es ist nun unser Ziel, zu Satz 1 und 2 analoge Aussagen fiir die Operatoren 
zu beweisen, die durch die Dirichletsche Randwertaufgabe der Potentialtheorie 
fiir beliebige glatte Flachen definiert werden. Im allgemeinen (von der Kugel 
abgesehen) ist jedoch die Halbgruppeneigenschaft nicht mehr erfiillt. Wir werden 
dadurch dazu gefiihrt, etwa zum Beweis von Lemma 2 nicht vom Additions- 
theorem auszugehen, sondern dafiir Voraussetzungen iiber Kompaktheit und 
Stetigkeit des Operators in den Vordergrund zu stellen. Solche Uberlegungen 
sind in der Arbeit [8] dargestellt. Dabei gelingt es, auf die in Lemma 2b voraus- 
gesetzte Reflexivitat des zugrunde gelegten Banach-Raumes zu verzichten, wenn 
man dafiir die Vollstetigkeit des Operators T(x) ausnutzt. Das ist fiir die An- 
wendungen auf die Potentialtheorie besonders wichtig. 

Wir werden daher als erstes (§ 2) unter Benutzung der Ergebnisse von [8] 
fiir eine Klasse von Operatoren zu Lemma 2 analoge Aussagen fiir nicht not- 
wendig reflexive Banach-Réume beweisen. Sodann zeigen wir im dritten Para- 
graphen, daB die durch die Dirichletsche Randwertaufgabe der Potentialtheorie 
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definierten Operatoren zu dieser Klasse gehéren. Im letzten Paragraphen schlieB- 
lich wollen wir den Bereich D(A) charakterisieren. Da uns jetzt Ergebnisse aus 
der Theorie der Fourier-Reihen nicht mehr zur Verfiigung stehen, geben wir 
eine Darstellung aller Funktionen g€(A) durch potentialtheoretische Uber- 
legungen und gelangen so zu Differenzierbarkeitsaussagen tiber g. Auf diese 
Weise gewinnen wir zu Satz 1 und 2 analoge Ergebnisse. 


§ 2. Ein Approximationssatz 
Es sei T(x) [0< x< co] ein Operator, der den Banach-Raum 8 in sich ab- 
bildet. 7(x) sei im Sinne der zugrunde gelegten Norm stetig. Das heiBt, es gilt 
(2.1) s— lim T(x) g= T(x) g 
fiir alle g¢¥. Es sei ferner 
(2.2) T(0) g=8. 


Wir setzen 
(2.3) 


und bezeichnen 


(2.4) Ag=lmA(h)g, g€D(A), 
falls dieser Limes existiert. Weiterhin sei 
(2.5) R(x, h) = T(x) A(h) + B(x, h), 


wobei B(x, h) gleichmaBig beschrankt und beziiglich / gleichgradig stetig in x 
ist. Besitzt T(x) die Halbgruppeneigenschaft (1.4), dann verschwindet B (x, h) 
identisch. Wir werden sagen, T(x) gehére zur Klasse 8, wenn T(x) den Eigen- 
schaften (2.1, 2,5) geniigt und R(x, h) fiir jedes x>0 kompakt ist. 


Dann beweisen wir 


Lemma 3. Der Operator T(x) gehore zur Klasse 8. Dann ist g€D(A) not- 
wendig und hinrerchend fiir 


|T7(2)g—g|=O(h) fir hoo. 
Zum. Beweis von Lemma 3 verwenden wir folgendes Ergebnis [8]: 


Satz 3. Es set T(x) ein stetiger Operator und 


R(x, h) g= 7 {T(x +h) — T(x)} g 


gleichmipig beschrankt, beztiglich h gleichgradig stetig in x und fiir jedes x>0O 
kompakt. Dann gibt es eine stetige Funktion p(t) mit 


T(x) g—@ =f el) dt. 
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Wir beweisen nun Lemma 3. Wenn g zu 2(A) gehort, folgt die Behauptung 
aus der Definition (2.4). Umgekehrt, es sei A(h) g beschrankt. Dann ist 


(2.6) T(x) A(h) g 

beziiglich / gleichgradig stetig in x. Auf Grund der Voraussetzung tiber B(x, i) 
folgt damit aus (2.5); daB auch R(x, 2) beziiglich h gleichgradig stetig in x ist. 
Somit sind die Voraussetzungen von Satz 3 erfiillt, und es gilt 


x 


(2.7) T(x)g—g=J p(t) dt. 


0 
Da (é) stetig ist, kann der Grenziibergang 
(2.8) lim + (T(x) g—g) 


x0 ¥ 


durchgefiihrt werden. Ag existiert also. Damit ist Lemma 3 bewiesen. 


Im nachsten Paragraphen wollen wir zeigen, da die in der Potentialtheorie 
auftretenden Operatoren zur Klasse ® gehéren. 


§ 3. Diskussion der in der Potentialtheorie auftretenden Operatoren 

Die geschlossene, glatte, orientierbare Flache 2 zerlege den dreidimensionalen 
Raum in das Innere G; und das AuBere G,. Es sei &(r, y) die Greensche Funktion 
fiir das Innengebiet G; und g (3) eine auf 2 definierte Funktion. g(3) gehére zum 
Banach-Raum %; wir wollen $ hier auf den Raum der stetigen Funktionen C 
oder den Raum der in #-ter Potenz absolut integrierbaren Funktionen L? 
(1 <p) beschranken. Dann erhalt man die Potentialfunktion U(x), die auf 
Q die Randwerte g(3) annimmt, durch 


(3.1) Ue) = [ ¢6)z, OWA ak. 


Zur Diskussion des Randverhaltens von U(r) fiihren wir in einer Umgebung von 
Q2 ein Normalensystem ein und stellen einen Punkt y dar als 


(3-2) r=3—on) 3€2, 


wobei der Normalenvektor nach auBen gerichtet ist. In einer Umgebung von 2 
setzen wir daher 


(3.3) U(r) = To) g. 
Der Operator T(o) ist beziiglich der Norm stetig in g, und es gilt 
(3.4) T'(0)g=¢. 


Da die Greensche Funktion © (r, 3) fiir alle 3¢€2 und r¢Q stetig differenzierbar 
ist, 1st 

al 
(3.5) R(x, h) = 4 {T(x +h) — T(a)} 
fiir jedes x>0 kompakt. Um zu zeigen, da T(x) zur Klasse & gehért, muB 
daher noch nachgewiesen werden, dab 


(3.6) R(x, h) = T(x) A(h) + B(x, h) 
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gilt, wobei B(x, h) gleichmaBig beschrankt und beziiglich / gleichgradig stetig 


in x ist. Das soll nun geschehen. 
Es sei Q, eine Parallelflache zu Q im Abstand y, die in G; liegt, und V(z) eine 
beliebige in G; erklarte Potentialfunktion (vgl. Abb. 1). Es sei ferner 


se (3.7) W(2, 75%) baw. W(Q, f; 2) 


we Gs die Potentialfunktion, die durch die Randwerte / (3) 
he auf Q bzw. auf Q, festgelegt ist. Dann folgt aus dem 
Eindeutigkeitssatz der Potentialtheorie im Inneren 

von 22, 


(3-8) V(x) = W(2, V(Q); x) = W(Q, V(Q,); x) 
% P mit 
oan (3.9) V(Q,) = W(Q, V(Q); Q). 
Der Aufpunkt x habe die Koordinate ge=A+ im Normalensystem (wir be- 


schranken uns auf eine Umgebung von Q, in der diese Darstellung méglich ist). 
Setzen wir dann 


(3.10) U2Q,5)0=TP (Af, 

so folgt aus (3.8, 9) 

(3.41) A+ u)g=T(A) T(u) g 

oder 

(3.42) TA+ ph) g=T(A)g+uT(d) A(u) + {P(A) — T(A)} T(u)- 
Nach (3.6) gilt also 

(3.13) B(A,w) =" {T*(a) — TA} T(W). 


Wir wollen zeigen, daB B(A, uw) gleichmaBig beschrankt und beziiglich « gleich- 
gradig stetig in A ist. 

Dazu miissen wir etwas weiter ausholen und eine Darstellung der Greenschen 
Funktion © (x, y) geben. Bekanntlich gilt 


1 


ee) waa epee de Ut) 
mit 
(3.15) Una = ap Ue) =U(.2). 


Zur Konstruktion von U(r, y) machen wir den Ansatz 


A6 UG wea i 12 yaa 
(3.16) (x, ») an | 708 tens 


und erhalten die Integralgleichung 


(3.17) 
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die eindeutig lésbar ist [5]. Um einen Einblick in das Verhalten von »(y, 3) zu 
bekommen, iterieren wir Gl. (3.17) und erhalten mit 


(3.18) (9,3) =K"y(y,8)3 Yn (0,3) =A" y(y, 8), 


=—y,+Ky,, 
(3.19) V4 4 4 


Ferner gilt mit 


1 1 a 1 0 1 
p: ( 5) —— 7 WT eer dk, 
£1 (65) aa ony ss ony |e 3) * 


(3.20) ; 3 ; 
ins) 55 | &ar6 SEN rs pee rg te 
Q 
(a) 1 
a= oe ne = at 
V1 (0, 3) 2m [n—3'| Ony [3—3'| 
(3.21) 
13) dF. 
Yn(¥),3) = Batt 1 (3, 3’) 


Aus (3.20) folgen die Abschatzungen 
$1(3, 3’) =O (Ig|3 —3'l) 
£2(3,3') =O (1) 


Die Funktion g,(3, 3’) ist mithin hdlderstetig beziiglich 3 [7]. Aus (3.21) erhalten 
wir fiir pees, 


(3.22) fiir |, —3'| 0. 


1 


an Grig, 12603) = Fn Tega + £1 & do) 
(3523) - 72 (803) =O (lg |4— ol) 
30 
a - 1 é 
Ging, V4 o> 8) = 8a (3, 40) real ans, | : a oo 3) dy. 
2 


Daraus folgt, daB ses V4 (so, %) existiert und beziiglich 3 hdlderstetig ist. Dann 
30 


gilt das auch fur cee 4 (30> 3). 
Ong, 
Es sei nun 
r=3—An(), 31 = § — (3) 
m=3—(A+tu)ni), @5,=V1—2u 46) +? K (3) dk. 


Dann vergleichen wir die Integrale 


f 7) D 
i= {it sai") Ony S(t, 5) dy 


(3.24) 


(3.25) 
P= [100 3 G1 (ry, 31) dF, 


wobei 1 (r,, 9,;) die Greensche Funktion fiir das Innere von 2, ist. 
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Als erstes ersehen wir aus (3.24) und (3.25) 
i onWay: (2) ’ 
(3.26) PHA) t= [ £6’ — wn) go OMe, a1) AE +o By (41) F- 
a ; 
Zur Diskussion von 1 (x1, 31) setzen wir 
eh ag aN r=|t| 


(3.27) w=u—H=s—s— (Ate) ng) +e), n=(4I 


und fiithren im Punkte 3 ein Tangenten-Normalensystem ein. Dann gilt 


g=0; 3 =(%,y,2); m3) =(0,0, 1) 
ue) ae Zy,1+0(%)), a = Vx? + v2 


2=O(e*); 22,2 44=O0(a). 


(3.28) 


Aus (3.28) folgt fiir wu, yO 


1 ee es ee uo 2) __ is 
7 Fel = lenleer| rere oe mall 
(3.29) : 
ieee .0() me =n0(4 
Be DE ! (aula Ps Ci i e) 
Weiterhin gilt mit 
(3.30) n(1) =n’) tuo, v= (O(a), O(a), 1) 
mt mh] _ (ne omy] bn 
ie Y pe "3 ti 
alla 1 meen F ‘Y 1 
(3.31) a pre = 8) Fa{t— 4) 5] +MO(, + 
Le er 
="0(- ae ca 


8t (Gy 5:31) = ey) we O ela 
Less iC i) 1 (5 s) 0 (Ig |3 —3']) 
83 (31, 31) = 8a (3, 3’) + uO (1) 


und fiir ) 4, 


7) ) 
an VY (bi) = 37 (0.8) + ol | 


(a) Z 
(3.33) Bae vi (D1, 1) = = V1 (9, 3) +0(23) 


0 4) 
oq V4 1 di =z val.) +m O(1). 
SchlieBlich erhalten wir aus den Integralgleichungen 


Ww=—HytKy,= Py, 


1 teeta 1 
Y¥4=— yy Ky Ty, 


(3.34) 
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wegen der Beschranktheit der Inversen (T)+ und 
(3.35) Kg=Kg+O(u), 


(3.36) yi Goan) == "4 (ood) LOU: 


ON3,, 
Wir erhalten daher aus (3.29, 31) sowie (3.14) 
0 


ce, ae 


und ersehen aus diesen Abschatzungen, daB der Operator B(A, ) gleichmafig 
beschrankt und beziiglich wv gleichgradig stetig in 4 ist. 


Damit ist gezeigt, daB die in der Potentialtheorie auftretenden Operatoren 
zur Klasse ® gehoren. Wir werden daher im nachsten Paragraphen die Ergebnisse 
des § 2 auf Potentialfunktionen anwenden und Aussagen iiber das Randverhalten 
dieser Funktionen gewinnen. 


§ 4. Das Randverhalten der Potentialfunktionen 
Es sei 2 wieder eine geschlossene, glatte Flache und U(r) eine im Innern von 
Q erklarte Potentialfunktion, die auf 2 vorgeschriebene Randwerte g (3) (g (3) € 8, 
B=C oder L? 1<p<oco) annimmt. In der Umgebung der Flache Q fithren wir 
wieder ein Normalensystem ein 


(4.1) £=3 — en(s) 
und setzen 
(4.2) U(r) = T(e) g (3) = T(e) g. 


Wir beweisen dann 
Satz 4. Notwendig und hinreichend fiir 
|T(e)g —gl=e(@) fir e>o 
ist g=const, d.h. T(o) g=const. 
sowie 
Satz 5. Notwendig und hinreichend fiir 
IT(e)g—gl=O(@) fire>o 
ast 


g@)= [76 74% mit v(3)CB. 


Als erstes beweisen wir Satz 4. Es sei g konstant. Dann ist auch T()g 
konstant fiir alle 9 (Eindeutigkeitssatz). Muithin gilt 


(4.3) IZ(e) s —gl=0. 
Umgekehrt, es sei 


(4.4) |Z (e) g — gl =9(e) 
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Dann verschwindet A g nach Definition. Es ist aber 


Ag=lim + (T(o) g— 8g) 


—0 0 
(4.5) F 1 é 


d.h., die Normalableitung von U(r) verschwindet auf dem Rande. Daraus folgt 
aber wieder U(r) =const, d.h. 


(4.6) T (0) g=g= const. 


Damit ist Satz 4 bewiesen. 
Zum Beweis von Satz 5 wenden wir Lemma 3 an und ersehen, daB die Be- 
dingung g€ D(A) notwendig und hinreichend fiir 


(4.7) IZ(e) g — gl =O) 
ist. Satz 5 ist daher bewiesen, wenn wir gezeigt haben, daB sich die Aussagen 
(4.8) geD(4) und ¢G)= fog) ah 

Q 


entsprechen. Wir wollen daher D(A) festlegen. 
Nach (4.5) gilt fiir alle g€¢ D(A) 
= id = me 
(4.9) Ag=zUl)| = re. 


Nach dem Verfahren der Integralgleichungsmethode laBt sich U(r) aus y ein- 
deutig konstruieren. 


Aus dem Ansatz 


(4.10) Ua)=3. fr) + 4k 
Q 


(4.11) Q 
y=v+Ky »,yES 


fiir die die Fredholmsche Alternative giiltig ist. Im Falle 8=C oder L® ist das 
bekannt. Fir G=L? gilt 


|K»|? =|fv(3') KG’, 3) dB |? 
(4.12) <J|rG') KG’, al 4B {Sf |K@’,a)| dee" 


er SAJ|r@)?|KG.a|4B, pg=pt+a, 


(4.13) Ky 


| Sel, 
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der Operator K ist also beschrankt. K ist auch vollstetig. Um das zu zeigen, 
zerlegen wir 2 in N Teilbereiche Q, und erklaren 


(4.14) joe in Q, 
| QO”  -sonst 
Wir erhalten also 
N 
(4.15) Ky=) Ky. 


Zum Beweis der Vollstetigkeit von K,;» wenden wir als Kriterium [3] an, daB 
es zu jedem ¢ ein 6(e)>0 so gibt, daB 


(4.16) [|Kir+t) — Kyo (@)PaR<e 
2j 
ist fiir alle |t]=¢<6(e) (t€Q) und alle y mit |r| <C. 
In unserem Falle erhalten wir, wenn wir um den Rand von Q, einen Streifen 
der Breite t abschneiden (2¢< 7) und den Rest mit Q* bezeichnen 


[|Kj»@ +t) — Kjr@)|?aB 


2, 
(4.17) S em |»)? |G, 3 + 1) — KG’, 8)| 4B aE + O(c) 
= [pl {sup B(’, t) + 0(75)}. 
Nun gilt, wenn K, im Kreis mit dem Radius t um 3’ ist 
BG.) =f|K3 +1) — KGa) 4% 
j 


(4.18) = i renee f oe 


QRAKr OF —O8NKr 


t 
=0(7+~,) =O) 
fiir alle 3’, wobei wir t=t* gewahlt haben (¢< 274). Mithin erhalten wir 
(4.19) Kj» G+1) — Ky @,=blOW, 


womit die Vollstetigkeit der Operatoren K; und damit auch die Vollstetigkeit 
von K gezeigt ist. Die Fredholmsche Alternative gilt also fiir Integralgleichung 
(4.11) und die dazu konjugierte. Es sei K’ der zu K adjungierte Operator. Dann 
bezeichnen wir die homogenen Gleichungen mit 

pel?; HEL? 

PIP a9: 


Die Funktionen g und H sind hélderstetig differenzierbar, und es gilt 


(4.20) O==geai-om 60H ke 


(4.21) H (3) = const. 
Aus ; 
1 ; 1 
(4.22) ee f Hise yt 
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folgt namlich U(z)|,=0 und daraus U=0 auben, Mithin gilt auch dU/o6n =O 


auf Q und daraus folgt U(r)|;= — 2H=const. 
Aus (4.9) erhalten wir ‘durch Anwendung der Greenschen Formel mit U(r) 


und 1 auf das Innere von 2 
(4.23) J jal =o, 


Daraus folgt aber die Lésbarkeit unserer Integralgleichung (4.11); wir erhalten 


1 1 
ae DG) cece ob Le « 
(4.24) Ue)=s, [6 path 
2 
Aus (4.24) ergibt sich schlieBlich 
1 1 
ee Vet Oly — # GWE LE, 
(4.25) SW=s5 [rele Ole 
92 


Insbesondere folgt aus (4.25) nach einem Ergebnis von CALDERON und ZyYGMUND [2] 
(4.26) Vi geL’, 


wobei wir mit VY, den Flachengradienten von g(3) auf {2 bezeichnen. 


Wir haben damit gezeigt, daB aus der Existenz von A g die Darstellung (4.12) 
fiir g(z) folget. Umegekehrt, es gelte (4.12). Dann erhalten wir in 


1 1 
4.27 U(r) = — ') —_—, dF, 
(4.27) @) = fro) anak 
Q 
eine Potentialfunktion U(r) mit den Randwerten g(3). Die Darstellung (4.27) 
ist eindeutig. Aus (4.27) folgt die Existenz von 


0U 


(4.28) FE Bea 


=7 3) oS 
und nach (4.5) also die Existenz von A g. 


Im Falle B= L?(2<p< oo) laBt sich Satz 5 folgendermaSen umformulieren: 


Satz 5*. Notwendig und hinreichend fiir 
IT(¢)¢—el= 


ist Vy gEL? (2< p< co) und g(3) totalstetig. 
Wir sahen bereits, dafB aus (4.25) die Existenz von VY, g folgt. Wir bilden weiter 


(4.29) 8 (i) — 8 (2)| S | |» | : 
1 (3 )| fl (3 )| eal lae= a 
Q 
Es sei nun |, dann folgt 
(4.30) a 
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oder 


(4.34) Is) el =0(, +e + )=0(4 +2) 


mit «>0 (wegen g<2). Setzen wir nun etwa t=73, so folgt aus (4.31) fiir alle 
i<2-3 


(4.32) |g (x) — §G2)| S 4a, — 30" 
und damit die Totalstetigkeit von g (3). 


Umgekehrt, es sei g (3) totalstetig und VY) g€L?. Wir wollen zeigen, daB dann 
die Normalableitung von 


(4.33) T(e) = UG) = [ 66’) 80s) 4h 


existiert. Nach (3.141) haben wir dazu zu diskutieren 
on 8) = a ta Bri oS ed 


ag [ay Rnd a ole wg Tt gi 


Im Integranden von — U, (x) 1a4Bt sich die Differentiation nach dem GauBschen 


(4.34) 


Satz umwailzen (vgl. [9] oder [6]). Man erhalt im wesentlichen ein Integral von 
der Form 


(4.35) J (Fos) Vora ah 
F L 
welches wegen Vy g€L? nach [2] fiir r +2 existiert. 


Zur Diskussion von ~~ U,(t) bilden wir 


ani, ee! IP) Oy (a9 re 
5 Ua t= a ff e6 ) saree a) An, ony =e dk, 
ae é a a a 
eens / ee eas on hanes ; 
=< Sg | LW) see ap ea et oe 
2 
mit 
(4.37) HG) = f 86") 5.90.8) ae 
a” 
Q 
4 oN ees p ates He AL : sO 1 : ee 
und BO Sie . Da OR y(3’",3') nach (3.19) wie eae singular 


wird, folgt analog [2] aus g€L? bereits V) fEL’. Die Differentiation kann daher 
auch in (4.36) nach dem Gaufschen Satz umgewalzt werden, der Grenzwert fiir 
y—>Q existiert. Wir haben damit aus VY, g¢L? und g totalstetig die Existenz von 
Ag gefolgert. Satz 5* ist damit bewiesen. 

Im Falle 1<f/<2 laBt sich die Totalstetigkeit von g aus der Darstellung 
(4.25) nicht mehr beweisen. Die Bedingungen in Satz 5* sind daher fiir die 
Existenz von A g hinreichend, jedoch nicht notwendig. 
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Wir bemerken noch, daB im zweidimensionalen Fall Satz 5* fiir alle # 
(1<p<oco) giiltig ist, da Bedingung (4.25) hier 


(4.38) e@=— frW) epi. 
€ 
lautet. Die Funktion g(3) ist daher fiir alle >1 totalstetig. 
Das Randverhalten der Lésungen des Dirichletschen Problems fiir den Halb- 
raum wurde von P. L. ButzEr mit Hilfe von Fourier-Transformationen disku- 
tiers | fa), 
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